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XCVIL. Deformation Textures in Face-Centred Cubic Metals. 


By E. A. Caunan, B.Sc., and C. J. B. CLzws, Ph.D. 
National Physical Laboratory, Metallurgy Division*. 


[Received July 11, Revised August 8, 1950.] 


INTRODUCTION. 


‘THERE have been several previous attempts to explain the deformation 
behaviour of polycrystalline aggregates in terms of the slip mechanisms 
operative when single crystals are deformed, but no wholly satisfactory 
treatment has been devised. In discussing these, reference must first 
be made to the work of Boas and Schmid (1931), which assumes the 
operation of the three most favourable slip systems for any particular 
crystal orientation and derives the end position which is stable for these 
three systems. This predicts tension and compression textures, for both 
_face-centred and body-centred cubic metals, which are in good agreement 
with the main features of the observed textures. Some details of the 
observed textures are, however, not explained. The treatment of 
Taylor (1938), mathematically more rigorous, has as its basic assumption 
the perfectly homogeneous deformation of the aggregate. Following the 
analysis of Mises (1928) that this demands the simultaneous operation of 
five special glide systems, Taylor has derived the systems and 
corresponding shears which conform to the principle of virtual work. 
The magnitudes and directions of the crystal rotations may be deduced 
from the shears selected, while the critical shear stress for a particular 
elongation is calculated from the sum of the shears and the single crystal 
shear-hardening curve. Good agreement with the experimentally 
determined tensile stress-elongation curve for a polycrystalline aggregate, 
and the main features of the tension and compression textures for 
face-centred cubic metals are obtained. The disadvantages of the 
treatment are first that homogeneous deformation rarely occurs in 
practice, and secondly that the heavy computation makes its application 
to metals of other crystal structures virtually impossible. 

Pickus and Mathewson (1939) have put forward a theory of the origin 
of rolling textures in which they assume that unequal participation of 
the three most favourable glide systems brings the deformation axes to 
symmetrical end positions, where the resolved shear stresses are equal 
and the rotations cancel out. The most favourable systems are those for 
which the product of the resolved shear stress and the cosine of the angle 
between the glide direction and the direction in which free flow may 
occur are highest. The relative probability of the end positions are 
considered to depend on the values of the product function at the end 
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positions. It is shown that all the observed rolling textures in face- 
centred metals fit in reasonably well with these ideas. The main criticism 
of this treatment and that of Boas and Schmid appears to be that, 
although the crystals oriented to correspond to the end positions will 
no doubt remain there and will slip without rotation, there is no complete 
explanation of how other crystals reach these positions. 

One further criticism of all the previous treatments may be noted. 
Although the simultaneous operation of several slip systems is demanded, ~ 
there is no explanation of the means by which the most favourable system 
or systems can apparently withstand resolved shear stresses greater than 
the critical value until the resolved shear stresses on less favourable 
systems have reached the critical value and deformation occurs. 


Main CoNcEPTS OF THE PRESENT TREATMENT. 


It is an experimental fact that considerable deformation can occur 
without, cracks appearing at the grain boundaries. In the ideal case of 
homogeneous deformation this demands the operation of five special 
slip systems, and even allowing for some inhomogeneous deformation, 
particularly near the grain boundaries, the critical shear stresses must be 
attained on more than one system at the same time. Moreover, apart 
from minor differences, for example, the observations on duplex slip. 
in «-brass crystals, Géler and Sachs (1929), the critical shear stresses. 
must be equal. Therefore the constraints imposed by the surrounding 
grains of the aggregate together with the applied stress must be such 
as to produce equal resolved shear stresses on the glide systems which 
operate. That is to say, if a grain deforms by the simultaneous operation 
of two systems then the resolved shear stress on these systems must be 
equal, and consequently it will be assumed that the effective stress may 
be represented by a tensile stress, T,, so oriented as to give duplex slip, 
v.€. its representative point on the stereographic projection lies on one 
of the boundaries of the unit triangle. If, however, the grain deforms by 
the operation of more than two slip systems then T, must lie at one of the 
corner points of the triangle, [100], [110], and [111], where there are 
respectively 8, 4, and 6 equally favourable slip systems. 

This concept, although probably an over-simplification of what must: 
be a complex system of stresses, does give a physical picture of the 
mechanism. Further consideration of the movement of T, leads to an 
explanation of the means by which T, reaches positions favourable for 
duplex or multiple slip without the resolved shear stress on the most 
favourable system exceeding the critical value. 

The most favourable system for a grain of applied tension axis in 
unit triangle IT of fig. 1 is the direction A on the slip plane of normal B. 
As the applied stress, T,, rises the resolved shear stress rises pro- 
portionately until it reaches the critical value at which slip would occur 
if the grain were an isolated single crystal. In general this single slip 
mechanism does not operate for all the grains in the polycrystalline 
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aggregate, and therefore the constraints due to the surrounding grains 
must prevent the resolved shear stress from reaching the critical value. 
Due to these constraints the direction of the effective stress, T,, moves 
away from the direction of the applied stress in such a way ae: ce 
resolved on the most favourable system is always below the critical 
value. 

In fig. 1 are shown the contours of equal resolved shear stress for the 
most favourable system, the numerical values being cosy cosA, where 
x and A are respectively the angles which the slip direction and slip plane 
normal make with the particular T, orientation. For orientation P (fig. 1), 
the critical resolved shear stress “C is reached when T,=C/0-46, and in 


order that the critical value shall not be reached as T, increases, T, 


Layee 


B 


Stereographic projection showing shear stress contours for a face-centred 
cubic metal. 


will move away to a point of lower cos y cosA value, 7.e., move down the 
contour gradient. When T, has increased such that T,—C/0-40, and slip 
has not occurred, then T, will lie on the 0-40 contour. The position on 
the contour will be determined by the stresses exerted by the neighbouring 
grains, and with a fine grained polycrystalline aggregate the general 
direction of movement may be taken as the path of quickest descent, 
i.e. perpendicular to the contours. On reaching a boundary of the unit: 
triangle, T, will travel along it in the direction of decreasing cosy cos 
until it reaches one of the corner points. These are minimum positions, 
and any further increase in T, must now result in slip, for in whatever 
direction T, attempts to move from the corner point some system will 
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attain the critical shear stress. Here, then, with the multiple systems 
available, deformation can occur without immediate misfit at the grain 
boundaries, and with the symmetrically disposed systems all exactly 
equivalent there will be no resultant rotation of the crystal. 

In some grains T, may not reach the corner point before slip occurs 
owing to insufficient lateral stresses, as perhaps in the case of a grain at 
a free surface. Consequently, rotation will occur in the usual directions 
of single or duplex slip, the latter taking place if T, has reached the unit 
triangle boundary. In addition to those at the surface, the grains most 
likely to slip before T, reaches a corner point are those for which the 
greatest lateral stresses will be required. These, which are the grains 
most likely to rotate, will thus be those furthest removed from the 
appropriate corner point, and conversely rotation will be least probable 
near the corners. In the instant at which slip occurs the lateral stresses 
are released and T, returns to T, or, if rotation takes place, to the new 
position of T,. With the cessation of this slip due to strain hardening, 
the movement of T, may begin again under the influence of a new 
environment which may or may not lead to rotation. Since inhomo- 
geneous deformation is implied in this treatment, T, need not have a 
unique value throughout the whole of a single grain, and the deformation 
corresponding to a single movement of T, may be as small as that 
associated with the passage of a single dislocation. Thus the deformation 
process is envisaged as step-like in character rather than as smooth and 
continuous. 


TENSILE STRESS-ELONGATION CURVE. 


From the concept of the deformation of a polycrystalline aggregate 
occurring when T’, in each grain has reached a corner point, the relation- 
ship between tensile stress and elongation may be derived from the 
corresponding single crystal shear hardening curve. First, however, the 
validity of considering T, in all the grains as having reached the corner 
points must be examined, since in this treatment the rotations which 
lead to deformation textures occur in precisely those grains, or parts of 
grains, for which this assumption is not true. The fact that the grains 
do not pull apart at the boundaries must imply that in the boundary 
regions at least. multiple slip occurs. Recent X-ray work on creep 
deformation in aluminium, Calnan and Burns (1950) has shown the greater 
complexity of the deformation processes in the regions adjacent to the 
grain boundaries. In view of the comparatively large size of these 
boundary regions, in which the lateral stresses must be sufficient to cause 
considerable displacement of T,, the lateral stresses elsewhere in the grain 
are probably not very much less, with the result that T, will be somewhere 
near a corner position. It can be shown that the tensile stress for a given 
elongation is not greatly different for multiple slip deformation (T,. at 
a corner point) and simple slip deformation with T, about 5° from the 
same corner point. Thus an approximation to the polycrystalline curve, 
within the limits of experimental variations, may be made on the 
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assumption that in all the grains T, has reached a corner point without 
vitiating the arguments employed in deriving the deformation textures, 
namely that in the interior of the grains at least T , has not reached this 
point. 

On the shear stress contour diagram, fig. 1, the peak position F is joined 
to the points of maximum shear stress on the boundaries of the unit 
triangle, C, D and E, by the lines of quickest descent. Since the mean 
path of T, is along the line of quickest descent, T, for all the grains of 
orientations contained in the area a,,,; between CFD and the [111] 
moves towards this last direction. For example, T, for the grain of 
orientation represented by P moves along its appropriate line for quickest 
descent until it reaches the [111]-{110] boundary at Q. From the shear 
stress contours in the adjacent triangle IV, corresponding to a different 
most favourable slip system, it will be seen that the boundary itself 
constitutes a line of quickest descent, and thus T, moves from Q towards 
[111]. Similarly, T, for the orientations represented by the area a@y19 
between DFE and the [110] moves to the [110], and for those in area a, 
between CFE and the [100] to the [100]. With an initially random 
orientation in the polycrystalline aggregate the proportions of grains with 
these end positions for T, are given by the ratios of the solid angles to 
which these areas correspond. It is however more convenient to use the 
ratios of the areas on the unit triangle, and since this involves an error of 
only about 1 per cent this method has been adopted. 

Anelongation e in the direction of the applied stress requires an elongation 
e/[cos* ¢(1+-c)—o] in the direction of a given corner point, where ¢ is the 
angle between the direction of the applied stress and the corner point, and 
o is the value of Poisson’s ratio chosen to give approximately constant 
volume over the range of deformation considered. A weighted mean 
value of cos: ¢ was calculated for each of the three areas referred to in the 
previous paragraph, viz. (cos? ¢),,;, (cos? ),49, and cos? $),99. For 
symmetrical slip it can be shown that the elongation in the direction of the 
corner point, i.e. in the direction of the effective stress for multiple slip, 
is S$ COS x9 COS Ay, Where s is the sum of the moduli of the shears on the 
individual operative systems, all of which have the same values of x) and 
Xo, the initial angles between the effective stress direction and the slip 
direction and ‘the slip plane normal respectively. Substituting the 
appropriate values of x), A, and (cos? ¢) in the expression 

€=8 COS Xo COS Ay[(cos? $)(1-++c)—o], 
gives the three values of shear s corresponding to the selected value of e. 
The critical shear stress for each of these s values may be determined from 
the experimental shear stress—shear strain curve for a single crystal. The 
effective tensile stress in the direction of the corner point to produce this 
shear stress is T,—critical shear stress/cos yp) cos Ag, which in turn implies 
a tensile stress in the applied stress direction 


T,,—(cos ¢) . critical shear stress/cos Xp COS Ag, 


where (cos ¢) is the weighted mean value of cos ¢ over the area considered. 
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The final tensile stress is obtained by adding these three values of T, 
weighted in proportion to.the areas defined above. The tensile stress has 
been determined by this method for aluminium up to extensions of about 
25 per cent, using the experimental single crystal results of Elam referred 
to by Taylor (1938). Details of the calculation are listed in Table I., and 
the resulting polycrystalline stress-elongation curve is shown in fig. 2 
together with the experimental points determined by Elam using the same 
material as that used for the single crystal work. The curve predicted by 
the Taylor theory is also shown. It will be seen that the curves derived 
from the Taylor theory and from the present treatment give equally 
satisfactory agreement with the experimental observations. Although 


Fig. 2. 
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Comparison of predicted tensile stress-elongation curve with Taylor’s 
predicted and experimental curves. 


the derivation of a stress-elongation curve must be regarded as an essential 
part of any deformation theory too much significance should not be 
attached to the exact agreement of prediction and experiment in view of 
the various assumptions made. The lack of agreement at higher elonga- 
tions calls for some comment however. Since the Elam curve is for load- 
elongation it may be expected to fall below the calculated stress-elongation 
curves in this region where reduction of cross-sectional area is now no 
longer negligible. Furthermore, in the calculations random orientation 
has been assumed throughout, whereas with increasing elongation some 
slight measure of preferential orientation must occur. 
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DEFORMATION TEXTURE IN TENSION. 


The development of deformation textures by marked grain rotations: 
implies deformations much more severe than those considered in the stress— 
strain relationship, where the departure from random orientation due to 
rotations was not taken into account. It has been proposed above that 
rotation will occur when the deformation takes place before T, has reached. 
a corner position. It will be that corresponding to duplex slip if T, has. 
reached a triangle boundary, and otherwise that corresponding to single 
crystal slip. In fig. 3 are shown the lines of quickest descent to the corner 
points, to the [211] direction, and to the maxima on the triangle sides. In 
addition the rotations for duplex slip are marked with arrowheads along 
the triangle sides. For an applied stress direction represented by p, the 
mean path of the effective stress, T,, is down the contour gradient to reach 
the [110-111] boundary between [110] and the maximum D. If deform- 
ation occurs before T, reaches this boundary, there is slip on the single 
most favourable system with consequent rotation along the great circle 
towards the operative slip direction A. This is shown by the arrow I. 
If T, reaches the boundary however, but deformation occurs before it has 
reached the [110] corner point, rotation corresponding to duplex slip occurs 
in the direction shown by the arrowhead on the boundary, producing the 
rotation of p shown by the arrow II. The single and duplex rotations for 
other typical points q, 7, t, wu, v and w, are deduced from similar reasoning. 
As the deformation proceeds, the orientation representing a particular 
grain (or part of a grain) moves across the unit triangle, and its probable 
rotations vary according to the particular region in which it is temporarily 
situated. Thus, with the treatment as oultined so far, it is possible for an 
orientation such as w to move into the region EF[100], typified by v, then. 
into the region CF[100], and finally into the .CF[112] region. One addi- 
tional feature mentioned in an earlier section will however modify this 
effect. The grains least likely to rotate are those which require the smallest 
lateral stresses for T, to reach the appropriate corner point, i.e. those 
orientations near the corner point. This will apply to the grains near the 
[111] direction in the DF[111] region, and those near the [100] direction in 
the EF[100] region. It will not apply to those near [110] in the DF[110] 
and EF[110] regions or those near [100] in the CF[100] region for the 
following reason. Consider an orientation similar to p near [110]. 
Although the effectiveness must travel only a short distance to the [110], 
any single crystal or duplex rotation which may occur increases the distance 
from the [110] ,and thus increases the probability of rotation, whereas in the- 
DF{[111] region near the [111] any rotation decreases the distance and the 
probability of rotation. The region between the [211] and the [111] 
directions is a special case in which, although the rotations lead away from 
the [111] corner point, the end point of the rotations, the [211], is still 
sufficiently near this corner point for the probability of rotation to be 
comparatively low. This should lead to a spread of orientations between 
[111] and [211], with a tendency for concentration near each of these: 
directions. 
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Fig. 4 shows the general trends of the rotations which give rise to the 
tension texture in fig. 4.__There appears to be no experimental work on 
pure tension textures with which these orientation spreads may be com- 
pared, and it has been the custom to consider the textures developed in 
cold-drawn wires as tension textures. While it is reasonable to suppose that 
in the wire drawing process longitudinal tension may be the predominating 
factor, it is probable that the ideal tension texture will be modified by 
radial compressive stresses, and discussion of this effect will be left until 
after some consideration of compression textures. 


Fig. 3. Fig. 4. Bien 
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from the shear stress tations in a_ face- tations after tension as: 
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tions of rotation due in tension. applied stress directions. 
to single crystal slip (tension texture). 


(+I) and duplex slip 
(+II)in aface-centred 
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COMPRESSION TEXTURE. 


The compression texture may be derived in an exactly similar manner to 
that described above, bearing in mind the fact that the rotation directions. 
for single and duplex slip are not simply reversed. For single slip the 
applied stress direction moves now towards the pole of the operative slip 
plane, while under duplex slip it moves towards the great circle joining the 
two operative poles. Thus for an orientation represented by the point r 
in fig. 3 and fig. 6 the single crystal rotation is in the tension case towards 
the operative [110] slip direction A, and in the compression case towards 
the pole of the operative (111) slip plane B. The duplex rotations are 
however merely reversed except in the small region between the [111] and 
[211] directions. 

Considering now the rotations in more detail, it will be seen that the 
orientations contained in the area F[111][110] will move towards the line 
F[110] and, in addition, an appreciable proportion of the orientations in 
CF[111] will also move towards this line. This will clearly denude the 
regions surounding [111]. Some of those in the CF[111] region, namely, 
those near to CF, will move across this line. In the region CF[100] there 
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is some tendency for rotation towards the [100] so that a small residual 
concentration of orientations may be expected near to [100]. Most of the 
orientations here however, together with those in the EF[100] region, will 
move in the EF[110] region. In this latter region the path of T, is primarily 
along the line E[110], and consequently duplex rotation towards [110] 
is the more probable. This will also have the effect of turning those 
orientations which have reached the F[110] boundary from the area above 
the line sharply towards the [110] direction. Fig. 7 shows these general 
rotations which would lead to a distribution of orientations shown in fig. 6. 


Fig. 6. 


B 


Lines of quickest descent from the shear stress peak (F), and directions of 
rotation due to single crystal slip (1) and duplex slip (+ ID) in a 
face-centred cubic metal in compression. 


The prediction of a [110] texture is of course amply confirmed by 
experiments on face-centred metals and the detailed agreement between 
the observed distribution of orientations, fig. 8 (Barrett and Levenson 
1940), and that predicted on the present treatment, is satisfactory. 


Rouuine TEXTURE. 


The simplest combination of stresses for the approximate simulation of 
the rolling process is that of a compression in the direction of the rolling 
plane normal together with a tension in the rolling direction (Wever 1931) 
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This has been adopted in the present treatment, which examines the ways 
in which both tension and compression end conditions may be satisfied 
simultaneously. The resulting rolling texture is therefore not that which 
would result from a mere super-position of the tension and compression 
textures. 


Fig. 7. Fig. 8. Fig. 9. 
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Stereographic projection showing the compression texture and the resulting 
spread in the rolling direction. 


The compression normal to the rolling plane will produce a texture such 
that, for the majority of the crystals, the direction of compression will lie 
within the shaded area shown in the unit triangle of fig. 10. To simplify 
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the ensuing discussion this texture may be considered as a spread of orien- 
tation along the great circle from the [110] direction towards the [311], 
shown dotted in the figure. Since the rolling direction is perpendicular 
to the rolling plane normal, the rolling direction must lie on the great circle 
for which the rolling plane normal is the pole. Thus the rolling direction 
lies on the circle ©, for the rolling plane normal in the [110] direction, and 
on the circle C, for the rolling plane normal at the other limit, the [311] 
direction. The great circles corresponding to the intermediate positions 
of the compression spread are obtained by oscillation, about the axis [112], 


Fig. 11. 
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of a great circle between the limits C, and Cy, and the distribution of rolling 
directions will be as indicated by the cross-hatching between these two 
great circles. At the same time, the tension in the rolling direction will 
tend to produce a texture similar to that illustrated in fig. 5, in which most 
of the crystals have the direction of tension at end-points of the form [111] 
[211] and [100], the tendency for a spread of orientation about these ena 
positions being ignored at present. Such directions are shown by the 
points A, B, D, E and F in fig. 10. The only one of these points which 
completely satisfies both a tension end condition and the compression 
spread is E, the [112] direction. Thus there will be a strong tendency for 
rotation to occur such that this condition is fulfilled, giving a strong 
(110)[112] texture. The [111] and [100] end positions of the tension 
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texture satisfy only the [110] end of the compression spread, but since this 
is the most densely populated region in the compression texture, there 
will be some tendency towards (110)[111] and (110)[100] components in 


the final rolling texture, which will be made up of the following preferred 
orientations. 


Direction of tension Direction of compression . 
(rolling direction) (normal to rolling plane) NOEEgS 
[112] [110] — [811] (110) [112] to 
(311) [112] 
[111] [110] (110) [111] 
[100] [110] (110) [100] 
Big. 12. 


OBSERVED 

Cu 97-7% REDUCTION 
(IWERONOWA AND 
SC HDANOoW) 


> T.D 
(111) pole figure showing predicted and observed rolling textures. 


The ideal (111) pole figure corresponding to these four textures is shown 
in fig. 11. The lines joining the (111) poles of the (110)[112] and (311)[112] 
textures are obtained by rotation of the (110)[112] texture about the [112] 
direction by an amount appropriate to the compression spread from [110] 
to [311]. In view of the fact that only the poles corresponding to (110) 
[112] and (311)[112] orientations and the spread between them satisfy 
both the compression spread and a tension end-position these will pre- 
dominate. Thus a semi-quantitative pole figure corresponding to these 
textures with some contributions from the (100)[111] and (110)[100] 
textures may be obtained as in fig. 12. 
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On the same figure is shown the (111) pole figure for copper after 
97-7 per cent reduction, obtained by Tweronowa and Schdanow (1934). 
The predicted pole figure shows a much higher degree of preferred orien- 
tation, but is otherwise in good agreement with the experimental obser- 
vations. It is perhaps worth noting that Barrett and Steadman (1942), 
in an investigation on the rolling of copper using both large grain poly- 
crystalline material and single crystals mounted in polycrystalline blocks, 
found (a) multiple orientations developed in nearly every grain, with the 
exception of those with initial orientation (1 10)[112] which were un- 
changed even after 97-6 per cent reduction, and (b) the entire polycrystalline 
texture developed in a crystal with the initial orientation (100)[001]. 
This supports the view that the (110)[112] orientation will predominate. 


Fig. 13. 


Stereographic projection showing ideal tension texture and the corresponding 
range of compression directions. 


DRAWING TEXTURE. 


The drawing process will be considered as equivalent to a tension exerted 
on the metal in the direction of drawing together with uniformly distri- 
buted radial stresses in the plane perpendicular to the drawing direction 
The longitudinal tension, considered to be the predominating influence, 
will tend to produce the tension texture already described, namel 
concentrations of orientations near [111], [211] and [100] directions The 
compression directions, being perpendicular to the tension direction must 
be at 90° from these end points, and therefore lie on the great ent C 
C, and C; respectively, fig. 13. A compression at any point on these eid 
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attempts to bring a [110] direction, the compression texture being regarded 
for simplicity as predominantly at [110] texture, into alignment with itself 
whilst retaining the tension textures. For any point on C, this requires a 
rotation of not more than 30° to the points A, B, and D, and on C, of not 
more than 45° to Aand E. On C©,, however, a rotation of as uel as 90° . 
to A would apparently be demanded, so it is far more likely that the 
orientations corresponding to C, will, except for those near to A, rotate to 
join the [110] end points Band D. This will reduce the density of orienta- 
tions around [211] in the tension texture, and correspondingly augment the 
[111] texture, thus producing the double fibre texture having [111] and 
[100] as the fibre axes. 

The experimental results of Schmid and Wassermann (1927) show the 
existence of this double texture in copper, silver and gold, although the 
proportions of the two textures vary from metal to metal. Aluminium 
however shows only the [111] component. 

It is not possible to predict from purely geometrical considerations these 
variations in texture, but in the present treatment the development of 
textures depends on the extent to which the effective stress, T,, can move 
before deformation occurs, that is to say, on the ability of the boundary 
regions of individual grains or sub-units of structure to withstand large 
stress gradients. This property may well be dependent on the temper- 
ature of the experiment in relation to the melting point of the metal and its 
purity ; the nearer to the melting point and the higher the purity, the less 
likely is the metal to withstand these stresses. Ifthe movement of T, is so 
limited that mainly single crystal rotation occurs the texture will be 
different from that developed in another face-centred cubic metal in which 
duplex rotation is the predominating feature. When mainly single crystal 
rotations occur the [100] component of both the compression and tension 
textures will be absent. Thus considering the development of the drawing 
texture, fig. 13, C, will be absent and a single fibre axis, the [111], will 
result. It may be significant that aluminium, which at room temperature 
is nearer its melting point than the other face-centred metals, copper, silver 
and gold, does indeed show this single texture. 


CONCLUSIONS. 


From this qualitative treatment of the deformation of polycrystalline 
agregates it has been possible to derive for face-centred cubic metals :— 

(1) a tensile stress-elongation curve which is in satisfactory agreement 
with experiment and with that previously derived by Taylor, 

(2) tension, compression, rolling, and wire-drawing textures which are 
in reasonable agreement with observation. 

An attempt has been made to give a clearer physical basis of the deform- 
ation process and since the above results can he obtained without laborious 
calculations it may be possible to extend the treatment to body-centred 
cubic metals. 

A very tentative suggestion is put forward to explain the variation of 
deformation textures among metals of the same crystal structure. 


1100 On the Deformation Teatures in Face-Centred Cubic Metals 


ACKNOWLEDGMENTS. 


The authors wish to acknowledge the encouragement and _ helpful 
comments received from Dr. N. P. Allen, and the assistance given by the 
Mathematics Division of the Laboratory in the calculation of the Resolved 
Shear Stress Tables. 

The work described above has been carried out as part of the research 
programme of the National Physical Laboratory, and this paper is pub- 
lished by permission of the Director of the Laboratory. 


REFERENCES. 


Barrett, C. 8., and Levenson, L. H., 1940, Trans. A.I.M.EH., 127, 112. 

Barrett, C. S., and StEapman, F. W., 1942, Trans. A.I.M.E., 147, 57. 

Boas, W., and Scumip, E., 1931, 1, Z. tech. Phystk, 12, 71. 

Caunan, HE. A., and Burns, B. D., 1950, J. Inst. Metals, 77, 445. 

GOLER, V., and Sacus, G., 1929, Z. Physik, 55, 581. 

Iweronowa, W., and Scupanow, G., 1934, Tech. Phys. U.S.S.R., 1, 64. 

Miss, R. v., 1928, Z. angew. Math. Mech., 8, 161. 

Pickus, M. R., and MatuEwson, C. H., 1939, J. Inst. Metals, 64, YRW 

Scumrp, E., and WASSERMANN, G., 1927, Z. Physik, 42, 779. 

‘TAYLOR, G. loss ae Inst. Metals, 62, 307; Stephen Timoshenko, 60th 
Anniversary Volume, p. 218, Macmillan, New York. 

WEvER, F., 1931, Trans. A.I.M.E., 98, 51. 


jer Todor | 


XOVIIT. The Energy Distribution behind Decaying Shocks.—l. Plane Waves. 


By M. J. Ligurutiu*. 


[Received July 11, 1950.] 


SUMMARY. 


A hypothesis is made concerning the degree of accuracy of a theory of 
decaying plane shocks due to Friedrichs. This hypothesis is carefully 
examined in the light of energy balance ; it is substantiated by the checks 
made on it, and in the course of the investigation the picture of the dis- 
turbance to the fluid as a whole is considerably improved. The work will 
be used in Part II. to clarify other phenomena involving decaying shocks, 
including some in steady supersonic flow. 

Most of the analysis refers to the following problem. A plane wave, 
headed by a shock, is produced by a piston, which is moved impulsively 
into homogeneous stationary fluid and then retarded arbitrarily until it 
comes to rest. One may usefully ask what is the whole history of the 
disturbance to the fluid if viscosity and heat conduction be neglected, since 
the modifying effects of these can be sketched in afterwards. The 
Friedrichs theory replaces the disturbance by a “ simple wave ”’ (with the 
fluid behind it in its original undisturbed state) and satisfies conditions at 
the shock neglecting the cube of its strength (see §§2-4 below). A precise 
hypothesis about the degree of accuracy of this theory is made in §1, but 
there is no experimental or logical certainty that this is correct at large 
distances from the piston. A comprehensive check is made as follows, 
arbitrary thermodynamic properties being assumed for the fluid. 

The energy of the simple wave at any time ¢ (§6) is compared with the 
work done by the piston (§5). The difference (§7) is a quantity of the 
second order in the initial shock strength, obtained with an error which (on 
the accuracy hypothesis) is of the third order. This residual energy R(t) 
can only represent energy of the fluid behind the simple wave (where 
disturbances are of the third order but their integrated effect may be of the 
second order). It is comparable with the total entropy gain S(t) of the 
fluid behind the simple wave, which is known from the theoretical variation 
of shock strength. In fact the relation 


RO=T(1-5 mp. — 3 pe) 8 5 APNE TS 


is found to hold, neglecting third order terms, where the suifix derivatives. 
are taken with respect to minus the volume at constant entropy (see §2. 
for notation). 


* Communicated by the Author. 
SER. 7, VOL. 4I, NO. 322.—NOV. 1950 4@ 
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This rather surprising relation is checked as follows. The residual 
energy is divided (§8) into that due to changes in entropy and pressure 
separately, the fotmer changes being known. Hence the integrated 
pressure over the fluid behind the simple wave is determined at all times. 
Since also these small variations of pressure satisfy the equation of sound 
(even though the entropy is non-uniform) and the pressure gradient at the 
stationary piston is zero, this determines completely the third order 
pressure waves in this part of the fluid. They are such as would be pro- 
luced if the shock emitted a pressure pulse of magnitude 


Las WPS MeV Pee cds 
‘eeeePs 


ak ap? i pe eae (2) 
backwards through the fluid at all times ¢, which was then reflected at the 
piston. The bracket in (2) is negative for most fluids but zero for mon- 
atomic gases. 

The wave (2) is next shown (§ 9) to consist of the reflection of the simple 
wave at the shock. The reflection coefficient in this process may be 
calculated exactly by considering the shock conditions with the local 
equations of motion, but only the third order terms in the reflected wave 
can be accurate after it has interacted with the simple wave. These are 
precisely (2), which completes the check by showing that the conclusions 
of the theory in the large agree with those in the very small. 

To demonstrate that the said third-order terms do not change by inter- 
action with the simple wave, the important, but apparently unrecognized, 
fact is used that, in fairly weak waves, a small variation of entropy is 
taken account of to a first approximation if the ordinary theory, at constant 
entropy, is applied with the velocity and pressure as independent variables ; 
in particular the flow behind a shock can then be obtained correctly to the 
third order in the strength. This is explained in §9. 

The extension to an “ N-wave ’’, produced by the piston being retarded 
through zero to a negative velocity, and then stopping abruptly to produce 
a second shock, is made in §10. The accuracy hypothesis, which is again 
substantiated by the checks made, now implies, however, that the motion 
of the rear shock after a long time is given less accurately by the theory 
than that of the front shock. Until the shocks have decayed to a small 
fraction of their initial strength, the accuracy is unimpaired; also 
equations (1) and (2) still hold. The pressure waves (2) are now aaa the 
combined effects of reflection of the simple wave at the front shock and 
transmission of it through the rear shock. But the reflections of these 
pressure waves at the piston gradually overtake the rear shock as time 
goes on and changes its path from that given by the Friedrichs theor 
(equation (79) shows by how much). Thus, when they are Taeisoune 
waves, they communicate their negative energy to the N-wave by slowin 
down the rear shock and so widening the N-wave on its low-pressare ie 
After a very long time, all the residual energy will be due to the increase si 
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- total entropy of the fluid with, ultimately, no change in pressure ; 7. e. 
R(co)=(1—TyP/TPy)S(co). The remaining energy is propagated to 
infinity in the N-wave (whose width increases as its amplitude decreases). 

It is shown in §11 that, if the piston has stopped on returning to its 
original position, then after an intermediate length of time the ratios of the 
energy in the N-wave, in the increased entropy, and in the rear pressure 
waves, are as 5:10: (—1) for air. The third form of energy is sub- 
sequently transferred to the first. On the other hand, if the piston stops 
much ahead of the original position, the first-named form of energy always 
predominates. 


§1. INTRODUCTION. 


DECAYING shocks appear, for example, whenever there is an explosion. 
Also, when any object moves through fluid faster than sound, the bow 
shock is observed to become weaker as we follow its course away from and 
behind the object ; this may also be called decay. A number of papers 
have been written on the mechanics of the process (see References). 
Some give detailed theories of its progress in particular cases. These 
mostly apply only when the shock is fairly weak, and hence only to the 
later stages of decay for a strong shock. Taylor (1950) treats the early 
stages of decay for a very strong explosion. 

A careful study of the published theories reveals certain doubts, 
difficulties, and even divergences of opinion. These are connected with 
considerations of energy and entropy. It is the purpose of this paper, 
with Part II., to clarify these points, and to answer such questions as 
(i) When a wave, with a decaying shock in front, and possibly one behind 
also, passes through fluid, what relationship, if any, is there between the 
rate of gain of energy, and of entropy, of the fluid behind the wave? 
(ii) If either is substantial, can a theory of the propagation of the wave, 
which makes use of assumptions valid only under adiabatic conditions, 
have any value? 

No reference to any but a single type of decaying shock will now be 
made until §10. The interval will be filled with a very detailed analysis, 
in this one representative case, of the questions broached above. 

Consider a very long tube indeed, of unit cross-sectional area, filled with 
homogeneous fluid at rest. Let a piston at one end (say x=0) begin to 
move into it impulsively, with a certain finite velocity, at time t—0, and 
be subsequently retarded in an arbitrary manner until after a certain time 
it is again at rest. Then at any instant, if the tube is long enough, the 
fluid in motion will be confined between the piston and a plane shock. 
The shock will decay in strength as it moves along the tube faster than 
sound. 

Viscous stress and heat flow by conduction will be very small, except inside 
the shock, and where otherwise they are largest (near the front of the wave) 
will not act on any element of fluid for very long. A theoretical treatment 


4AG2 
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which neglects them, and replaces the very thin shock by a plane surface of - 
discontinuity, should therefore be quite satisfactory. Then the motion is 
governed by Huler’s equations of continuity and momentum, the Rankine- 
Hugoniot relations at the shock, and constancy of entropy for each 
element of fluid after the shock has passed it by ; the boundary condition 
is that the velocity of the fluid in contact with the piston is that of the 
piston. 

No exact mathematical solution is known, even to a particular case of 
this problem. An approximate theory of Friedrichs * should be good, 
according to. plausible hypotheses, when the velocity of the piston is always 
rather less than the velocity of sound in the undisturbed fluid, so that the 
shock is always weak. But there has been no experimental verification, 
or mathematical proof, of the validity of the approximate theory. It 
would not perhaps be difficult to produce either if the approximations were 
required to be good only within a limited distance of the piston. But 
there are serious difficulties, even in the mathematical discussion, at 
arbitrarily large distances from the piston, where any terms neglected may 
have had an arbitrarily long time to influence the propagation of the wave. 
(It may be noticed here that, the weaker be the shock, the longer does it 
take to decay to, say, one half of its strength.) 

The Friedrichs method is to determine the “simple wave’’, at the 
entropy of the undisturbed fluid, which is compatible with the motion of 
the piston. This is taken as the flow behind the shock. (In a simple 
wave the velocity and the thermodynamic variables are all functionally 
related. Any value of one of them is propagated through the fluid at a 
constant speed, namely at the local speed of sound relative to the local 
fluid velocity. The simple wave can exist only when the specific entropy 
is uniform. In the present problem it will terminate with the sound wave 
emitted by the piston as it comes to rest. Behind this, on simple wave 
theory, the fluid is once more homogeneous and at rest.) The shock is 
calculated approximately as a plane surface moving in such a way that the 
physical quantities in the simple wave thereat take those values which the 
motion of the shock into fluid at rest would produce, with the cube of the 
shock strength neglected. The calculations were given by Friedrichs only 
for the case when the fluid is a perfect gas with constant specific heats, but. 
the extension to a fluid with general thermodynamic properties is a fairly 
simple matter. 

It would be impossible to satisfy the conditions at the shock more 
exactly ; for example, the specific entropy, which has to be taken uniform 
in the simple wave, must actually increase at the shock by an amount of 
order the cube of the strength. 


—_———. ES 

* Friedrichs, 1948, pp. 220-231. This theory of the phenomenon, which is the 
subject of the present §§ 1-9, will be constantly referred to as “the Friedrichs 
theory”, even though in the paper cited many other phenomena were also 
investigated. 
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Friedrichs (1948) is cautious about claiming any precise degree of 
accuracy for this theory, but an attempt at estimating the errors may be 
made. These errors arise from using the incorrect value of the entropy, 
and from the fact that only a second order boundary condition can be 
applied at the shock. These considerations suggest the following 

Accuracy Hypothesis.—The errors in the Friedrichs theory are compar- 
able with the changes produced by altering the simple wave in two 
respects : (i) the specific entropy of the whole wave is altered to its known 
maximum, 2. e. by an amount of order the cube of the maximum shock 
strength ; (ii) the velocity distribution is altered, say by changing the 
velocity of the piston (continuously, and even after it has actually come 
to rest) by an amount whose ratio to the velocity of sound is of order the 
cube of the maximum shock strength. 

If this hypothesis is correct, the Friedrichs theory is of considerable 
value at all distances from the piston, as will be seen in §4. Even after 
many multiples of the duration of motion of the piston, the position of the 
shock (relative to its position if it has been an acoustic wave) is in error 
only by a factor one plus a term of order the square of the maximum shock 
strength. Even though this factor multiplies a quantity which grows 
indefinitely, the approximation may be regarded as highly satisfactory. 

The accuracy hypothesis might be questioned on two grounds, both 
based on the fact that the total entropy gain of the fluid as a whole, as a 
result of the passage of the shock, is of order the square of the initial shock 
strength *. (i) The effect on the wave of the terms neglected, if it accumu- 
lates as time goes on, might be held as a consequence to be of order greater 
than the cube of the maxiumm shock strength. (i) It might be expected 
that the total energy of the fluid is defective (on the Friedrichs theory) 
by an amount of order the square of the maximum strength (and this is 
true), and that this indicates errors of such an order in the wave. 

The author rejects objection (i) on the grounds that the form of the 
accuracy hypothesis already permits propagation of errors along the 
diverging characteristics, and hence their progressive accentuation ; 
thus, in the model simple wave used for comparison, the effects of the new 
value for the specific entropy accumulate as time goes on. The defective 
total energy mentioned in (ii) is the subject of this paper, which will 
analyse it in detail, assuming the accuracy hypothesis. It is sufficient 
here to remark that, though errors within the simple wave itself could not 
account for this deficiency in total energy, the region behind it (which the 
Friedrichs theory gives as undisturbed) is permitted, on the accuracy 
hypothesis, to include disturbances of order the cube of the maximum 
shock strength. The integrated effect of these over the extensive region 
behind the simple wave is sufficient to account for the discrepancy in 


energy. 


* That this is so, although the gain in specific entropy of any element is of 
order the cube of the strength, is because the weaker shocks take longer to decay. 
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But the hypothesis remains unproved, and it was to make a further test 
of its correctness, and in particular to make a detailed critique of obj ection 
(ii), that the further deductions from the hypothesis, given below, were 
made. If they had led to a contradiction, or to physically unplausible 
results, the hypothesis would have had to be modified. — In fact, the results. 
are neither impossible or unlikely, and the considerably fuller picture of the 
flow as a whole, which this investigation in the large provides, is confirmed 
in quantitative detail by an investigation of the local reflection of the 
simple wave at the shock. 

The deductions which are to be made have already been sketched in the 
Summary. At any time t the total increase in energy, kinetic and internal, 
of the fluid in the simple wave over its internal energy in the undisturbed 
state—briefly, ‘‘ the energy in the simple wave ’’—will be calculated on the 
accuracy hypothesis. Now the energy in the whole fluid must equal the 
work done by the piston. Hence, if the accuracy hypothesis is true, the 
difference must represent energy of the fluid behind the simple wave. 
From the expression for this residual energy, the pressure distribution in 
this region is deduced, to a first approximation, as a third order function 
varying only relatively gradually. The same pressure distribution is. 
then deduced in a totally different manner, as a result of reflection of the 
simple wave at the shock followed by reflection at the piston. 

The author has thought it worth while to set out the work for a fluid with 
perfectly general thermodynamic properties, for three reasons:—(i) The 
proportionality of residual energy and entropy shown in equation (1) 
(in the absence of which the theory might find itself in difficulties), and the 
agreement in expression (2) as obtained by two methods, then appear as. 
the result of the cancellation of a large number of terms, in a manner 
impossible to believe fortuitous. (ii) The work shows up explicitly what 
terms arising in it are due to departures from adiabatic conditions * 
(roughly those bringing in the temperature). (iii) It is desirable to place on 
record a statement of the Friedrichs theory for fluids other than perfect 
gases. 

The “ tail” behind the simple wave is found to be a compression wave, 
in which the slightly reduced pressure directly behind the simple wave is. 
gradually made up to its undisturbed value. Technically this compression 
wave must ultimately form a shock, but it is so weak that this could not 
happen for a very long time, and need hardly be considered in detail ah 
One may also notice that the increased entropy in the general neighbour- 
hood of the piston, which the theory leaves as permanent, would gradually 
be flattened out, for a real fluid, by conduction of heat, causing a further 
slight increase in the total entropy, though not in the total energy. 


* See § 3 for a discussion of this point, and of the errors in a paper of Martin 
due to his ignoring it. 


} The indication that an N-wave must ultimately form, even in the present 
problem, is, however, interesting. 
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§2. THe Stmpte Wave. 


The notation to be used for the various physical quantities, at points in 
the undisturbed fluid and elsewhere, is indicated in the table below. 


Quantity Undisturbed Disturbed 
Pressure P P+ 
Specific volume Vv V-—v 
Velocity 0 u 
Velocity of sound A=V\V/Py a=(V—v),/(ep/dv), 
Temperature Laer ee iY Rae Ee 
Internal energy E E+e 
Specific entropy 0 28 


Derivatives, in states not that of the undisturbed fluid, will be indicated 
in the manner usual in thermodynamics ; thus (@p/dv), indicates a deri- 
vative at constant entropy. But the values in the undisturbed state of 
derivatives along an adiabatic (7. e. at constant entropy) will be denoted 
by suffixes and by the use of capital letters ; Py is the value of the deri- 
vative just mentioned in the undisturbed fluid, Pyy the corresponding 
second derivative, and so on*. Suffix S will also be used to denote 
differentiation with respect to entropy at constant volume in the undis- 
turbed state. 

The Euler equations of continuity and momentum are 


Ov ‘dv Ou ou ou OD ee 
a T% a, ap Mo Mepees OE ns el Oi aE . (3) 


Hence at constant entropy 
a (dv dv du ow du a ov 
mente +uz) = G +035) an Cees ) me 
a (ov Ov Ou Ou Ou a Ov 
Vo (5 +e 7) aE & +05) tt (a Vp 7) ee 


A solution is therefore given by putting 


eae 
= i, Orpen. ru dar Maal «trey (0) 
to satisfy the first equation, so that the second becomes 
Ou Ou 
a = ee Paine ce oO 
Ot +(u-+-a) as 0, (6) 


showing that wu is constant along straight characteristics dz/di=u—a. 
This solution is the “‘ simple wave ”’. 
* Notice that Py is therefore minus the derivative of pressure with respect to 


volume (and Pyyy minus the third derivative, but Pyy plus the second 
derivative). This is convenient as making the quantities all positive for most 


fluids. 
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Friedrichs showed that it was convenient to regard the simple wave as 
specified by a single function o(€), which he defined as follows. € isa 
parameter constant on each straight characteristic, namely the value of x 
where the characteristic produced would meet the axis t=0. It is always 
negative in the problem of this paper. Then, writing wt+-a—=A(1 +o), it 
follows, since wu and a are constants on any characteristic, that o and € are 
functionally related. Actually, in the present problem, it is more con~ 
venient to regard é as a function of o than vice versa, for it is a problem in 
which this function is one-valued, unlike some considered by Friedrichs. 
Then the equation of a characteristic is 


v= &(o)+At(l+o). woke) he Se ema) 


Since o is a non-dimensional measure of the disturbance to the fluid, its 
cube may always be neglected in the Friedrichs theory. 

All the physical quantities are related in the simple wave ; all will now 
be expressed in terms of o, neglecting cubes of disturbances. Firstly, 
by (5), 


. °° A+ Avy N. A i} A Ay 2 >. 

w= | ae d= Fot5 (i re (8) 

whence ; 
ar uta—A Pe ] Ay 1 i} Ay Ayy 2 

a = (pt ets (tap te)” = Sag) 


Replacing A by VP¥, and its derivatives accordingly, (9) may be rewritten 


Evy aN (ae Paw Peay ie 
eA Os (sve 7 ae an) be (19) 
oP 1 Pp PoP 
Genes ey vies 4 v- vvv \ 2 
whence v (eh, E (5 : VP,. PEsnt al guider ore CLE 
2AP Ps pap a 
and u= Ww | (: see MSs AERA |S P 
VPyy al VPyy Poy ) : | Papa 


The coetiicient of o* in u vanishes in the case considered by Friedrichs, of a 
perfect gas with constant specific heats. 

According to the accuracy hypothesis the errors in (7), (11) and (12) are 
such as would be produced by adding terms of order o? (or iather the 
maximum of this quantity) to those in o and o-. 


§ 3. CONDITIONS aT THE SHOCK. 


The values of disturbed quantities behind the shock are related by 
equations expressing conservation of mass momentum and energy for an 


element of fluid as the shock passes through it, say with velocity U 
These are 


. 


Ug: U Wl 
Veen ae (P+ p)u= = (3u?-+e). ie 1 3) 
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Elimination of U and w gives the Hugoniot relation between the thermo- 
dynamic variables _ 


Cee 2) USSCay alga Meee each ks) 


It is convenient to expand (14) in series, using that P+p=(de/dv),. The 
best form is as 


[Ey+}(Eyyyo?-+2Eyygvs-+ Eyes?) Ju= Kyo+ Egs 


mail 
+ 54 (Eyyyu?+ 3Eyygv7s-+ 3Eygcvs? + Eggos*) + O(v' +85), . (15) 
where (in this one equation only in the whole paper) suffix derivatives 
refer not to the undisturbed state but to a state with specific volume ,and 
entropy halfway between those ahead and of behind the shock. (In this 
way quadratic and quartic terms on both sides of the equation are avoided.) 


This gives the entropy as 


Hiyv¥in.g toh vives 
cok Set eee an ae UG) 


the error being O(v*) if Pyy and T are taken in the said intermediate state, 
but O(v*) if (as will be sufficient in what follows) they are taken in the 
undisturbed state. 

Equation (16) is a convenient form of the relation between the thermo- 
dynamic variables behind the shock ; it states how near to the adiabatic 
this relation is. The velocities w and U are deducible, by (13), as 


u=(pv)}, = VCore se ew Fates: © (LT) 


It should be noticed that it follows from these equations, with (16), that 
in the expansion of uw in powers of p or v alone, neglecting cubes of the 
disturbances, the adiabatic relation may be used, but that this is not so in 
the expansion of U, which is 


U=V[Py+4PyyotsPyyye?+ Aor +P,(s/v)+ a) Ket ts iE 


oi Pyy Pyyy Ty Pyy 2 3 


since P,=Ey,=Ty. 

Martin (1947) falls into the error of neglecting the effect on U of the 
entropy change. Hence second order terms in his calculations are wrong. 
One result of this mistake is that he obtains a reflection of a simple wave 
at a shock, whereas in fact, for the special case of a perfect gas with adiabatic 
index 5/3 which he considers, the reflection coefficient is at most of the 
fourth order, as will be seen in $9. Friedrichs avoids the mistake by using 
the explicit solution for a perfect gas of the Rankine-Hugoniot relation. 
The T,/T in (18), which will recur in many subsequent formulae, is a 
constant reminder that in one place, at least, the Friedrichs theory does 
have to take account of the entropy change at the shock. 
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§ 4, DETERMINATION OF THE SHOCK. 
i 


The shock condition must now be expressed, using the “ simple wave 
values (§2) for the various physical quantities behind it. Then v is 
related to the fundamental variable o by (11), and the expansion (18) 
becomes 


U=A(1+}o-+2o?), 


h fee Py 3 Py Ty 
eee P4ESERGP \ Velo Donel 


For a perfect gas with constant specific heats \ is simply $. 

At time ¢ let the variable o take the value o, at the shock, so that o, is a 
meagure of the strength of the shock at time ¢. Then the condition 
dz/dt—U along the shock, with (7) and (19), gives 


(19) 


d 
(E(o,)+At) 7 +A(L+0,)=A(L+4o, +09), . - - (20) 
dt 
whence ($0, —Aoy)A F— + At=— €"(04). ee eee eee 
O71 


The initial condition for the shock is x=0 at time t=0, and so by (7) 
€=0 whent=0. If oo, initially, so that o, is a measure of the maximum 
shock strength, then (21) may be solved exactly as 


Wen (la 2AGT\? oe ala) 
At—a( = ) ienae ee or es Sty" 


but, since cubes of o have been neglected in the formulze on which this 
equation is based (without which it would be impossible to satisfy the 
shock conditions completely) the only significant terms in (22) are 


2s 8A pe 12 (% 

AS é - =) Joe dé(o)+ + of dec). . .-. (28) 

This tells us the time ¢ at which the non-dimensional measure of distur- 

bance o takes the value o, just behind thé shock. It is easy to verify that 

the condition (17) on wu is also satisfied to the required order, by comparing 

it with (5). 

Here the significance of the errors in (23) which the accuracy hypothesis 

permits should be noticed. They are terms like 


% % rOo 
or? | a dé, o;1| o* dé, | o dé, ete. 
Their ratio to the terms included is always of order on, although the first: 
two become indefinitely large as time goes on and o, decreases towards 
zero. In consequence, the value of oc, for given t is also in error only by a 
factor 1-+O(o2), as stated in §1. 
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§ 5. Tue Piston. 


It is necessary to have an expression for the motion of the piston in 
terms of the fundamental function &(). The former should logically come 
first, as being the cause of the whole motion, but it is much more convenient 
to work from o throughout. It will be sufficient for our purposes to neglect. 
the square of o in determining the time ¢ at which the piston emits the 
characteristic specified by o. 

At the piston dv=w dt; hence, by (7) and (12), 


, 2Py 
do(é'(o)+ At)+A dt (: +o(1— vEx)) A) Yin pk (4). 


whose solution, still neglecting o?, is 
At=—&(o)+o&(o)+ + au Pp 
taking x=0, and so €=0, initially ae fuea 

Where the piston comes to rest w=0, and hence, by (12), c=0 (neglecting, 
on the accuracy hypothesis, o3). The “‘ simple wave ”’ will be taken as the 
region satisfying 0<o<o, exactly, 7. e. lying between the shock and the 
characteristic whose speed is that of sound in the undisturbed fluid. This 
characteristic is not emitted precisely as the piston comes to rest, but (as. 
will be seen in §9) very slightly earlier. 

The work done by the piston against the pressure of the fluid, neglecting 
a3, is 


{ Bie vutdeidey de 


| ode, eee (26) 


2PPy acd 3Py a) (ox 2 | 
= dé |, 
veel |. ade (5 ae Ss oak Joe é 
Te. (6) 
where (12) and (25) have been used and, in integrating of(c) by parts, 
the fact that it vanishes at both limits. 


To the order of approximation employed, (26) must represent the total 
gain in energy of the fluid at any time after the piston has come to rest. 


§6. ENERGY OF THE SimpLE WAVE AT TIME ¢. 

The energy of the simple wave will now be calculated, on the accuracy 
hypothesis, at any time ¢ after the piston has come to rest. Under 
adiabatic conditions the increased kinetic and internal energy per unit 
mass of a disturbed element of fluid is 


pe feted 
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which in the simple wave is 


; 2PPS 1 Py  PyPyvv , 2P% ) | 27 
du?+Pov+43Pyv?= po o4 Gs ate VPyy rel p24 EPee o ’ ( ‘ 


neglecting o°. 
The total energy in the simple wave 0 <o<o, obtained by integrating 
(27) multiplied by da/(V—v), with x given by (7), 1s 


2PP, (2: l° Pare Pabeacihe nee ) | ( Je ueeby ) 
oe Gi (5 + Yi) me Paqe | cPPee) ai as 
roy i Ope 
, Lots y 12 as eae 
x (E(a) + At) do= Fp | f. o dé+ Athe? (5 + PPS 


aie eS Ath : ete 
ae VP ay Pz, j, o d€+ 3°71 ( ) 
With At expressed in terms of o, by (23), this energy becomes 


2PP, 


"Oo 


Ee Epa ee is odé+6r) ode 


O71 


1 2P% 3Py AES] ( i : ec ts aie 

= = : o* dé+20 a dé ake t) 
a (5 shy ae) gram we oe Easy An in| (29) 
Since in (28) Atis multiplied entirely by terms o7, o?, etc. *, the error in 
applying (23) to obtain (29) is at most of order o? (arising from integrals 
like | a oF dé). Thus, on the accuracy hypothesis, (29) is the correct first 


and second order expression for the energy of the simple wave. 


§7. THe RestpuaL ENERGY. 


The difference between (26) and (29) must therefore represent energy in 
the region behind the simple wave. Subtracting, one obtains for this 
residual energy 


Pha /leeoP2 ewer. un Pope 
MAT Aare (5 PP avPe Bae —6r) 


x (| o dé— jo, | ode). 2 Vala. ERS 
Equation (30) has been proved only at times t after the piston has come to 
rest. But, before this, At=O(1), and study of (23) shows that then 
| (a1) |=O(o), so that (30) is O(o3). The residual energy at this early 
stage may also be shown to be of the third order, and so (30) remains true 
in the sense that the difference between its left- and right-hand sides is 
always a third order quantity. 
eee 
* This is because the simple wave has been considered as terminating with the 
characteristic on which c=0 exactly. This point will be examined further in § 9. 
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It will be observed that all integrals in (26) and (29) involving the range 
O<o<c, have disappeared by the cancellation of eight distinct terms. 
This suggests that R(t) depends only on the properties of the shock previous 
to time ¢. The total entropy of the fluid at time #, say S(t), also has this 
property, and a comparison of the two is clearly desirable. 

Unit mass of fluid, on passing through the shock, gains an entropy s 
given by (16), and subsequently retains it. To a first approximation mass 
passes through the shock at a rate A/V ; also vis given by (11) as 2Pyo/Pyy. 
Hence the rate of increase of total entropy of the fluid is 


dS__ A Pyy /2Py\3 , 
G7 ¥ it pe) Peter ee (31). 


Hence, using the first approximation to (23) to integrate (31), 


S(t) = Li i= 7 ist (Po) [ts(-= —3| a dé do, 


to} 1/4 oy 


eA So Go 
a 2 2 dé— : ba a 2), 
3 VIP2, (s ie o* d&—4o, |. o az) (32), , 


Comparison of (32) with (30) shows that R(t) and S(¢) are proportional 
at all times (neglecting third order terms). Substitution for » in (30) 
from (19), and division, gives their relationship in the form already stated 
in the summary, namely 


SEE age Te) S(t) (1 bis), 

The “‘ coincidences ”’ necessary to obtain this relatively simple relation- 
ship already add to one’s confidence in the Friedrichs theory. On the 
other hand, the relationship is not so simple that any rough-and-ready 
arguments, based on familiarity with the thermodynamics of systems in 
equilibrium, would have been able to predict it. 


§8. StaTE oF THE FLUID BEHIND THE SIMPLE WAVE. 


According to the accuracy hypothesis, the velocities, and the distur- 
bances to the thermodynamic variables, are of the third order in the region 
behind the simple wave. The total energy is R(t). Since kinetic energy 
terms are of the sixth order, they are negligible in building up the total 
energy, and so 

] At 
R= 7 |. Epa mee oN os Fie eae (33) 
Here the limits are not exactly correct for describing the region behind the 
simple wave, but the error involved in using them in (33) is still of the third 
order, since the true boundaries are exactly parallel to x=0 and x—At 


respectively. 
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Now, to a first approximation, 


[8 elt 
e=TetPypSTs-eP ("st apes ol ee 
Pog, 
“A. PT 
et | pde= = {VR()— (t > *) vsio} , 
0 “2 
VPpaNaw 
=r (44 — 4) 80 7 Bae aa 


by (34), (33) and (1). Thus the fluid behind the simple wave cannot be 
entirely at its original undisturbed pressure unless the bracket in (35) is 
zero (so that the bracket in (1) becomes 1—TyP/TPy) *. 

For a perfect gas with constant adiabatic index y the bracket in (35) is 
3y—5. Thus it vanishes for monatomic gases, but is negative for other 
gases, even if they are only approximately perfect, since for these y is 
always considerably less than 3. The bracket is probably not positive 
for any fluid. Thus for liquids the second term, corresponding to the 
relative decrease of compressibility with decreasing specific volume, is 
always considerable, exceeding about 3 even at high pressures, and at least 
6 at normal pressures. On the other hand, the first term may be written, 
in terms of the coefficient of expansion « and the compressibility x, as 
4xVi/c,«. Here 4V/c, is limited by the law of Dulong and Petit to less than 
about 3x 10-7 cm.?/dyne.deg. But «a/x, the pressure increase due to 
raising a rigidly enclosed gram of liquid through one degree, is not as 
much as twenty atmospheres at normal pressures, as (it is seen) would be 
needed to make the first term in the bracket attain to the lower limit 6 of 
the second term. On the other hand, at high pressures the coefficient of 
expansion falls off rapidly, while the compressibility increases, indicating 
that here too the bracket in (35) is negative. It will be assumed, through- 
out what follows, that the bracket is non-positive +, as has here been 
indicated may be true for all, or nearly all, fluids. Then, by (35), the 
average pressure behind the simple wave is below the undisturbed value 
(or equal to it, to the third order, for a monatomic gas, a gloss which will not 
be repeated henceforth in similar contexts). 

Now in this region of small variations of pressure, entropy and velocity, 
those of pressure at least are governed by the equation of sound. For by 
(3), with small disturbances, ; 


dv Ou Ou Op 
ae ae == (); Ot ot Vee Ue ° . . . . (36) 


* Nor can the fluid be at its original undisturbed specific volume, which would 
mean that (1) is replaced by R(t)=TS(f). 

t In words, the conjecture is that the slope of an adiabatic, as graph of pres- 
sure against volume, increases at least as rapidly as the fourth power of the 
temperature. 
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But entropy is effectively independent of t, because constant for a given 
fluid element. Hence dv/dt=(V2/A?)dp/dt, and, therefore, 


OP _ yx OP 
BP an Ont 


(37) 


Also, by the second of equations (3), at the piston (where x 0, and u=0 
at all times ¢ after the piston has come to rest), 


(Gr0x) ep 0 ee i ai Ke *(88) 


Equations (35), (37) and (38) determine p completely (to the present 
third order approximation). By (37) and (38) 


p=f(At+a)+f(At—z) MA ae eesti (39) 
for some f. Equation (35) then gives 


= Wty _ WPor) sy 


fv) dy=¥ (4 (40) 


Py 


Hence, behind the simple wave, 


EEL Gld) a) 


where S(é) is given by (31) *. 

The first term in square brackets in (41) represents. a pressure wave 
propagated in the direction x decreasing, opposite to that of the shock. 
‘The second term is the reflection of the first at the piston. At v=Af the 
first is simply S’(¢). Since the shock satisfies x=At with a proportional 
error which is uniformly small, we may say that the pressure fluctuations 
behind the simple wave are such as would be produced if at each time ¢ the 
shock sent backwards a pressure pulse 


T /,Vy_VPw\d8__ PY (,VTv_VPwv\ 2 (4p) 
16A dime 2AY PERE Tn Dp ort 


p Py 


as stated in the Summary, which is subsequently reflected at the piston. 

Since 8’(t) decreases monotonically to zero, it follows that after a long 
time the pressure (41) is only considerable within a fringe or * tail”, say 
A(t—t,)<a<At, behind the simple wave, wherein 8’(}(t—a/A)) is still 
significant. Thus only the reflected wave is of importance after a long 
time. The main simple wave “ over-expands ”’ slightly to a negative 


pressure 


¢ 


Ay ( VTy 1 VPyy 


miltan pe) 8), a eed etic Co) 


: ih 
* It may be remarked that, by (36), the velocity takes a form similar to (41), 
‘ut with an additional factor V/A and with the sign of the first term in square 


brackets changed. 
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for large ¢ (and to twice this value for small ¢), and this pressure is recovered 
during the passage of a very much wider compression wave. At every 


point fae 


t>0 
thus the initial undisturbed pressure is finally recovered. The entropy, 
however, remains permanently altered, unless, as pointed out at the end of 
§1, heat conduction is taken into account. 


$9. CAUSE OF THE BACKWARD-EMITTED PRESSURE WAVE. 


It has been inferred, by a somewhat roundabout investigation of the 
disturbance to the fluid as a whole, that the shock emits a pressure wave 
backwards, the pulse emitted at time ¢ being given by (2), or (42). This 
will now be verified directly by a study of conditions just behind the shock. 
This will show that the pressure wave may be regarded as a reflection of 
the simple wave at the shock. 

Immediately behind a plane shock, which moves into undisturbed fluid, 
the pressure may locally be divided into two waves each moving with the 
local speed of sound relative to the local fluid velocity, one in the direction 
of motion of the shock and the other in the opposite direction. The ratio 
of the pressure gradient across the latter to that across the former may be 
called the reflection coefficient of waves at the shock. It may be calculated 
exactly for shocks of arbitrary strength, and is proportional to the cube 
of the strength for weak shocks. It does not seem to have been given 
explicitly in the literature, though an essentially similar problem, of two- 
dimensional steady supersonic flow, was solved by the author (Lighthill 
1949 a). 

The reflected pressure wave is of the same order of magnitude as the 
entropy variations behind the shock. However, it remains unchanged to 
this order as it passes through the simple wave into the fluid behind. 
This is because pressure and velocity, unlike the other variables involved, 
remain insensitive to variations of entropy to a first approximation—i.e. 
ignoring products of variations of entropy with those of other quantities. 

For sound waves this fact is known, and has been used in §8 ; but this 
section is concerned with waves of finite, though not large, amplitude. 
Special cases of the result have been used by the author in his work on the 
flow behind a stationary shock (Lighthill 1949 a) and on the diffraction of 
blast (Lighthill 1949 b, 1950), but he has realized its truth in the generality 
given above only more recently. 

The proof fills only one paragraph. In the full equations of continuity 
and momentum for an inviscid fluid, derivatives of the pressure with 
respect to all three space variables occur. But the only derivative of the 
density which occurs is the “ total *’ derivative following a particle of fluid. 
Now, neglecting conduction and radiation of heat, the entropy is constant 
for a given particle of fluid, so that this total derivative of density may be 
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replaced by that of the pressure, divided by the square of the velocity of 
sound. The four equations are then homogenous in the first derivatives of 
the pressure and the velocities, the coefficients of these derivatives being 
the velocities and various thermodynamic variables. Hence, in replacing 
the latter by their values as functions of the pressure on a particular 
adiabatic, only variations in entropy multiplied by variations in velocity or 
pressure are neglected in the equations. So, if such products are negligible, 
it is sufficient to solve the equations of continuity and momentum for the 
velocity and pressure as if the entropy were constant. 

Thus in the case of plane waves the Riemann theory can be used to . 
determine the velocity and pressure. It will perhaps add to the reader’s 
conviction if this fact is deduced independently from the exact character- 
istic equations for plane waves (Courant and Friedrichs 1948, p. 82): 


ae V—v 


a 


au dp on a au+a, du=++— dp on Lara 


dt 


d 5 
ds=0 on en 


dt 
(44) 


If in these equations (V—v)/a be replaced by the value as a function of p 
which it takes for one particular value of the entropy, say (V—v)/a=G'(p), 
changes in wu and p will be given errors of order the entropy multiplied by 
those changes themselves. So if variations in wu and » are first order and 
those of entropy are third order, the errors in taking 


dx } 
u+G(p)=constant on Fete: | 
(45) 
u—G(p)=constant on —=u—a, 
dt j 


(to replace the ordinary Riemann equations) will be of the fourth order. 
(But if G(p) were replaced by a function of v which is identical along the 
adiabatic s—0—in fact by the integral (5)—the error would be of the 
third order.) 

Now consider equations (44) at a point directly behind the shock. 
Here the shock conditions (14) and (17) hold. Given any one of wu, p, v, 
these determine the others. Also at adjacent points on the shock the 
conditions hold. Hence changes of wu along the shock are related in a 
specific way to changes of p; that is du/dp along da/dt=U is a known 
function deducible from conditions (14) and (17). To be precise, it is 


7) __ (P—(de/dv),,)v+ (0e/Op) ,p 
Ndp/) sno  — 24(P+2p—(0e/0r),)  ” 
SER. 7, VOL. 41, NO. 322.—NOV. 1950 4H 


(46) 
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But if dp/dt=« on dx/dt=u+a, so that du/dt=—«(V—v)/a, and if 
dpjdt=B on da/dt=u—a, so that du/dt=B(V—v)/a, then 


(=z) _ (dufdt)arja—0 
shock 


dp ~ (dp|dt)axja=v 
U—uta/«(V—v)\ | —Utura (“—) 
2a (eerese i ae : 
“= Ueeaca —U+u+a 
eile tag OE 
S Veo(eresU) pe ee ee 


a (a+u—U)B+(a—u+U)z° 


This gives the reflection coefficient exactly as 


ae atu—U a [du a du 
B - a—u-+U (hs V—» (Foca) ft Ved ae a seis 


For a weak shock the numerators in (48) are respectively of the first 
and second orders in the strength, so that the whole is of the third order. 
For since, by (16), behind the shock 


while by (17) w=(pv)?, it is easy to deduce that , 
dp Pyy Pyvy , PyyTy  5P% 3 
(Fe) sox Vv 7 (1+ 5p apy ee T STP, Sr)? ea )), shoe 
whereas, by expanding a in powers of v, 
@ A Pyv Pyvv  Piv\ 
eon a (i+ ap" IP, ot & a) #+0(0%)) «ae oe OL 


Hence the reflection coefficient (48), to the third order, using (17) and (11), 
is 


Oo ae 2 (Gaps a) 2Py At Py 4 Vy VPvv 3 
Bae sTP 32P2 (Gen ni / 5 wae Ae ape ) 
| (52) 

The bracket in (52) is recognized as that of (2). Notice, however, that 
the term in Ty/T due to entropy change has arisen in (52) in a manner 
different from that in which it appeared in calculating the residual energy. 
It is due to the second order term in dp/du, not to that in U. This fact 
seems to the author to reduce still further the likelihood that the agree- 
ment of the check on Friedrichs’s theory is fortuitous. 

Now f, the value of dp/dt on dx/dt=u—a, is known approximately from 
the simple wave theory (on which, however, «=0). To a first approxi- 
mation p=(2P%/Pyy)o; also, by (7), if da/dt=u—a, 


—Adt=(é(c)+At)dot+Adt, .. . . . . (53) 
d fe —2A 
so that B= a = CPP wl aan: aay 
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Hence by (52), and since for c=o,, to a first approximation, by: (20), the 
denominator of (54) is —40, dt/do,, 


a 


PB ( pees) »do, oss 


eco eae erat Dias dia 


To see that this value of (dp/dt)q» 4-440 18 propagated through the simple 
wave unchanged to a first approximation, consider two adjacent character- 
istics dx/dt—u—a. The time interval At between them remains unchanged 
to a first approximation. At the shock the pressure difference between 
them, along a characteristic of the other family, is «4t, and the velocity 
difference is —G’(p)xAt by (45). As p varies on each characteristic, (45) 
shows that wu varies by the change in G(p), so behind the simple wave (for 
example) 


Au—(—G'(p)aAt)=G'(p)(Ap—a« At). ee ae (0) 


But if the differences between the two characteristics are again taken 
along a characteristic of the other family, then 4u=—G'(p)4p. Hence, 
by (56), dp—a«At as before; that is, behind the simple wave, dp/dt on a 
characteristic dx/di=u-a is still given by (55). 

Since the shock is, to a first approximation, a characteristic of the family 
dx/dt=u-—+a, it has been shown that the flow behind the simple wave is 
the same as if the shock sent out a pulse at time ¢ whose time derivative 
is (55) and which can therefore differ by at most a constant from expression 
(42). (Strictly it has only been shown that the waves travelling towards 
the piston are of this form ; but waves travelling in the opposite direction 
must represent the reflection of the former waves at the piston if the 
boundary condition there is to be satisfied.) 

This checks (42) to within a constant. The constant if present would, 
by the theory of §8, make the total energy of the fluid tend to infinity 
with time. But we can also see that it must be absent from the point of 
view of the present section. For a complete third order solution of the 
equations of motion and boundary conditions could now be given by 
specifying the pressure and velocity fields behind the shock as the combin- 
ation of the simple wave, determined (to the third order) by the motion of 
the piston, and its third order reflection in the shock (with zts reflection at 
the piston). Now the value of the pressure at the piston immediately 
after it first moves may be calculated independently, (i) from the shock 
conditions, and (ii) from the simple wave theory. The difference must be 
twice the pressure in the reflected wave at this initial stage (twice because 
the reflected wave is itself reflected at the piston). The two quantities 
are in fact the integrals with respect to u, from 0 to w, of (50) and (51) 
respectively. Their difference is easily verified to be twice (42). This 
check of (42) at t=0, with the previous check of its time derivative, 
completes the investigation of the accuracy hypothesis concerning the 
Friedrichs theory. It has been found watertight at every point. 


4H2 
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The complete third order solution, which could be obtained in the 
manner just described, would add a negative pressure distribution to that 
given by simple wave theory. Now, at large distances from the piston, 
the reflection of the simple wave at the shock becomes negligible (by (42)), 
and only the reflection of this reflected wave at the piston modifies the 
original simple wave. ‘This, travelling in the same direction, turns it into 
another simple wave (the entropy being effectively uniform at these large 
distances from the piston), with reduced pressure distribution. This 
‘means that the characteristics are slowed down and, therefore, that which 
now travels with velocity A must be one emitted slightly before the piston 
comes to rest. This is the one which was taken as the boundary of the 
simple wave in §5, for the purpose of calculating its energy in §6. 

It is now clear that the form of the accuracy hypothesis adopted in §1 
was necessary as well as sufficient. For a modification of the simple wave 
of the kind described really takes place at large distances from the piston, 
while near it the entropy is modified as stated. It is not surprising that 
the errors inthe Friedrichs theory should be comparable with those obtained 
by introducing both these modifications, and no others, uniformly. 


§10. THe Decayine N-WAVE. 


After the careful examination of this single class of problems, each 
involving only one decaying shock, it is important to discover whether the 
conclusions for phenomena involving more than one shock are substantially 
identical or not. 

The problem which will be investigated is like that treated in §§1-9, 
except that the piston is continuously retarded (from its initial impulsively 
attained positive velocity) through zero speed to a negative velocity (so 
that it is moving away from the fluid) and is then instantaneously stopped 
(not necessarily in its original position), and subsequently held at rest. 
The piston, on being instantaneously stopped, emits a second shock, which 
moves along the tube behind the first while both gradually decay in 
strength. The motion of the front shock is exactly as in the problem 
previously studied. 

Friedrichs (1948) only treats a special case, but his method is quite 
general. The flow between the shocks is taken to be the simple wave, at 
the entropy of the undisturbed fluid, which is compatible with the motion 
of the piston. The fluid behind the second shock is in its original undis- 
turbed state, and the shock conditions thereat, as at the first shock. are 
to be satisfied to the second order in the strength. The pressure, after 
rising discontinuously at the front shock, falls in the simple wave mono- 
tonically to below its undisturbed value, to which it then rises discontin- 
uously at the rear shock. This behaviour may be represented graphically 
as something like an N, whence the wave is called an N-wave. 

Considerations as in §1 lead to the formulation of an accuracy hypo- 
thesis for the Friedrichs theory of the N-wave, with identical wording to 
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that used in §1. Conclusions as in §§2-9 can then be drawn, and will be 
displayed in brief, which check the hypothesis and improve the picture of 
the motion. 

But it must be remarked at once that this hypothesis, which is about to 
be checked, does not forecast that the motion of the rear shock will be 
given to a similar degree of accuracy by the simple wave theory, at 
distances from the piston so large that it has decayed to a small fraction 
of its initial strength, as was indicated for the front shock by the theory of 
§4. For the front shock, conditions on one side are known exactly. Asa . 
consequence, conditions on the other side can be specified, as in § 4, with 
an error of order the cube of the local strength. But matters at the rear 
shock are different. The conditions behind it are known, on the accuracy 
hypothesis, with a possible error of order the cube of the maximum shock 
strength. This error communicates itself to the conditions ahead of it, 
which determine the position of the rear shock in the simple wave; and 

~when the shock has become very weak its motion may therefore be given 
quite inaccurately by simple wave theory. It will be seen later, however, 
by use of the improved physical picture, how to find more closely what will 
happen, and hence, for the time being, the consequences merely of the 
Friedrichs theory of the N-wave will be set out. 

Taking the state behind the rear shock as undisturbed, it is easy to show 
that the shock conditions in the simple wave thereat have exactly the 
same form as for the front shock. For example, (19) holds (but notice that 
o is negative at the second shock). At time ¢ let the variable o take the 
value co, at the second shock. Then the differential equation (21) holds 
equally with o, for c,. But the initial condition is different. If c=, at 
the piston when it comes to rest, say at t=t, = —&(c3)/A, then these are 
the initial values for the second shock, at which therefore (instead of the 
solution (22)) ’ 


ee (l= 20a fo. dé(c) 1 om 2 , 
eae ( om ) | «, (1 —2Ac)* a 2 (<) At | Beer eae 
or, rejecting unreliable terms, as in (23), 
— 2 ; 8A i 2 a 2 1 2 
Apo i Bs =) | ° dé(o)— z (or | x dé(o)+ $02 (s)) . (58) 


(But since the expression (19) for U, and hence also the coefficient of 
dt/dc, in the differential equation, may contain errors O(0}), the equation 
(58) is not reliable for indefinitely small o,.) 

The work done by the piston before coming to rest is the same as in 
(26), but with o, replacing the limit of integration 0, and hence with an 
additional term 402£(c3) inside the square brackets, since in integrating 
oé by parts the integrated term does not vanish at the limit oy. 

At time ¢ the simple wave is now specified by the inequality o,<0o<o,, 
i.e. it lies between the two shocks. Hence in the expression (28) for its 
energy the range of integration must be replaced by this one ; and the 
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integrated terms Atso? and At}o} must be replaced by At(4o0?— 403) and 
At(4o3—403) respectively. Then, replacing At by its expression (23) or 
(58), according as it is multiplied by powers of o, or 9», the energy in the 
simple wave becomes 
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and the residual energy, obtained by subtracting (59) from the work done 
by the piston, becomes 


PEEP y lie 2Pe te aby eee ry 
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x I | ) o8 dé 90, | « dé— $e, | car (60) 


(Again, (59) and (60) may be inaccurate, since (58) has been used to obtain 
them, when oc, has become too small.) 

On the other hand the rate of increase of the total entropy S(t) of the 
fluid is given by the term (31), due to the gain at the front shock, plus a 
similar term, obtained by writing o, for o, and changing the sign, due to 
the gain at the rear shock ; the latter is also positive since o,<0. Inte- 
grating with respect to the time t from the initial value S(0)=0, and using 
the first approximation to (23) or (58) to express dt in terms of do, or dog, 
the expression 
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results. Equation (61) is correct to the second order, on the accuracy 
hypothesis, at all times, since after o, has decreased to a small fraction of 
its initial value, changes in S(t) due to the rear shock cannot be significant. 
_ Comparison of (61) with (60) shows that when both are true the equation 
(1) relating the entropy S(¢) and the residual energy R(é) still holds for an 
N-wave. For this period of time the theory of §8, which is based simply 
on (1) and the accuracy hypothesis, therefore holds without change for 
the N-wave. Conditions behind the rear shock are, to the third order, as 
if at time ¢ a pressure pulse (2) was emitted, in the direction x decreasing, 
from approximately x= At, 2.e. (to within the degree of approximation 
permitted) from either of the two shocks. In terms of oc, this pulse (2) 
becomes 


2) 
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(Compare (42).) 
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The first term in (62) corresponds to the reflection of the simple wave at 
the front shock, as worked out in §9. This reflection is transmitted 
through the rear shock together with an additional wave (corresponding 
to the second term in (62)) which may be regarded as the transmission of 
the simple wave through the rear shock. (The characteristics in the simple 
wave are gradually caught up by the shock, which, as it does so, touches off 
a small part of their energy as a transmitted wave travelling in the 
opposite direction.) 

The analysis of the rear shock etyangt the last paragraph 1 is as follows. 
Ahead of it, by §9, 


dp du V—v dx F 
de and jean along aii =uta. - (63) 
But, by (54), immediately ahead of the rear shock, 
dp es 4A du V—v : dx 
adi =6,—=(2P5/Pyy) o, dildo, 5 and dE — ee Be, along at =U ae 


(64) 
Behind the rear shock quantities are of the third order and only their first 
approximation is of interest; this is given by the theory of sound, on 
which the shock is a characteristic and, along it, say 
dp du Vv 
== Spe a Ona week eck i ees OO 
| ria eet eo Co) 
Now the pressure difference and velocity difference across the rear shock 
are so related that changes in them satisfy (50), since the coefficients in (50) 
represent the state behind the shock correctly to the second order. Taken 
with (65) this implies that, ahead of the rear shock, 
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with third order terms neglected on the right and fourth order terms 
neglected in both numerator and denominator on the left. Substituting 
in the left-hand side by use of the information in (63) and (64), and 
maintaining the same degree of approximation, it becomes 
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by (55), (64), (19), (51) and (50). 


so that 
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Integrating (68), it may be deduced that the pressure distribution along: 
the rear side of the rear shock is as in (62). These pressures (along this. 
approximate characteristic) may be taken as propagated acoustically 
along characteristics of the other family, ¢.e. towards the piston—as has. 
been deduced previously by energy considerations alone. 

It may now be seen what are the third order disturbances behind the: 
rear shock which make its equation of motion (58) unreliable after it has. 
greatly weakened. By this time the actual pressure pulses (62) have. 
themselves become weak ; it is the reflection at the piston of pulses (62), 
emitted at a much earlier time which are significant. (This reflected wave 
is given by (41) without the first term in square brackets.) Notice that 
ultimately, though very gradually, the whole reflected wave catches up: 
the rear shock. For, by (58), the distance of the rear shock behind an 
imaginary wave emitted by the piston with velocity A on coming to rest, 
which is practically —Atc,, tends to infinity like a small multiple of #}. 
So all acoustic waves travelling behind it ultimately, though slowly, catch it 
up. Their reflection from the shock is negligible, since the reflection 
coefficient is of the third order. They cannot be transmitted through it, 
since the time decreases on both characteristics ahead of the shock. Their 
energy goes into an alteration of the shock path. Thus, since they are 
rarefaction waves, their incidence on the shock makes it weaken more 
slowly, and makes the N-wave become wider (on the low pressure side) than 
on the theory which neglects them, It is only in this indirect way that 
they transmit their negative energy to the N-wave. 


§11. Tot ENnERGy DiIsTRIBUTION AFTER A Lone TIME. 


This section will treat of the approximate energy distribution, as given 
by the Friedrichs theory, after the shocks have weakened to a small 
fraction of their initial strengths, and also of the subsequent modifications. 
due to the gradual overtaking of the rear shock by the pressure wave 
behind it. 

The energy distribution on the Friedrichs theory, as t> 00, will be 
divided into that remaining permanently in the N-wave (whose width 
increases as its amplitude decreases), that due to increase in entropy of the 
fluid through which the wave has passed (with, ultimately, no change in 
pressure), and that in the pressure wave behind the rear shock. It is only 
after these have become almost stationary at their limiting values that the 
gradual process of transfer of energy from the third form to the first 
becomes important. For this reason the said limiting values will first be 
determined. 

The expressions which have been obtained for all these energies with 
t arbitrary are simplified by letting t+ oo. For then, by (23) and (58), 
o, and o, tend to 0 (naturally, because the shock strengths tend to 0) 
The energy (59) in the simple wave becomes 

2PPy fi - 2 2 
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The energy due to increased entropy, after the original pressure has been 
recovered, by (34) and (61), is 


TP 4P3 TYP\ 2% 
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and the remainder of the residual energy, corresponding to the pressure 
wave behind the rear shock, is 
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by (1). For the problem of §§1—9, without a rear shock, (69), (70) and (71) 
are correct at all sufficiently large times. For that of §10 the gradual 
transfer of the energy (71) to the N-wave will make the energy of the 
N-wave, after a very long time indeed, equal to the sum of (69) and (71), 
namely to 


2PPy [| (% Ee poeee T 
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The energy (70) is positive for all fluids ; the energy (71) is probably 
negative for all fluids (§8), though zero to the third order for monatomic 
gases. The energy (69) may be positive or negative, and normally the 
first term in the square brackets therein will be of a higher order of 
magnitude than the others, and will dominate this part of the energy and 
make it more important than the residual energy (70) and (71). We have 
then an N-wave with one “ upright ” of the N always taller than the other, 
and the internal energy on that taller side of the wave due to greater 
volume changes on that side predominates in importance at all times. 

But there is one problem, of special interest, in which the first term in 
(69) does not dominate the others. This is the case when the piston is 
stopped precisely on reaching its initial position. It is interesting as 
being directly analogous to the problem of the steady two-dimensional 
supersonic flow past an aerofoil (neglecting boundary layer effects). The 
condition is that over the path of the piston 


ie u(dt/dc) do=0, 


whence, by formule of §5, 
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Then the energy (69) in the N-wave *, on substituting for A from (19), 
becomes 
2PPy 1 PyTy % 2 ‘ 
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* This energy is still due to an inequality in the “ uprights” of the N, though 
a smaller one than when the first term in (68) dominates the others. 
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‘Thus, if the piston comes to rest in its initial position, the approximate 
distribution of energy, when the shocks have decayed to a small fraction 
cof their initial strengths, as between the N-wave, the increased entropy 
of the fluid as a whole, and the pressure waves behind the N-wave, is in the 


ratios (by (74), (70) and (71)) 
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‘These depend only on the physical properties of the fluid. 


For a perfect gas with constant adiabatic index these ratios become 
Sy—3 1 3y—5 (76) 
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and for air they are 
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ive ae ae" 


Notice that in this one case the total energy of the fluid is TS(0o). 
‘Thus, after S(t) has reached its limiting value, the energy is that 
associated with the total entropy change when volume changes are small 
and integrate out to zero. That the integrated specific volume changes 
should be zero is inevitable, to a first approximation, since the space 
occupied by the fluid is unchanged if the piston returns to its initial 
position. Notice however that the effect is achieved by cancellation of the 
-decreased volume in the N-wave with the increased volume (due to 
increased entropy without change of pressure) near the piston. Observe 
also that, before S(t)==S(0o), the volume changes in the N-wave are too 
large, compared with the entropy changes, to be treated linearly (and in 
fact the energy exceeds TS(t), being TS(co) at all times after the piston 
has come to rest). 

As time goes on the energy in the third form is transferred to that in the 
first. Ultimately (one may suppose) in all cases, and not merely approxi- 
mately, all the energy which has not gone into increasing the entropy of 
the fluid without change of pressure is in progress to infinity with the 
N-wave, after that portion which has been released and reflected at the 
piston has once more overtaken it. (It seems unlikely that, even theoreti- 
cally and after vast times, this plane N-wave would be damped by viscosity 
‘or heat conduction. For it follows from (7) that do/da~(At)-! as t +00, 
so that the gradients of all physical quantities have an approximately 
uniform value in the N-wave, decreasing as the reciprocal of the time. 
The energy dissipation in all time might then be expected to be bounded by 


a negligible integral of the form O ( | pt? . it) , where is the viscosity). 


* For a liquid at ordinary pressures they are practically 0: 1: 0, owing to the 
low value of xP=P/VPv. (But if the piston does not return to its original 
position, the energy in the N-wave still predominates.) 
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The path of the rear shock, as it gradually moves behind the path given 
by the Friedrichs theory; thus increasing the width of the N-wave on the 
low pressure side to permit the negative energy (71) to pass into the N-wave, 
would be calculable to some extent, if necessary. The errors in (18), as a 
formula for giving the conditions at the rear shock, if the flow is disturbed 
behind it, are due, to a first approximation, to taking A incorrectly, and to 
taking U and v not relative to the velocity and volume behind the shock. 
Thus, if suffix 1 refers to conditions behind the rear shock, the error in (18) 
is u,+a,—A(1+Pyyv,/2Py). When the shocks have weakened consider- 
ably, so that s is practically zero behind them, this expression becomes 
u,+Vp,/A approximately, a quantity constant along characteristics 
dx—A dt, and equal by (41) to 
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Adding this to (19), we see that the differential equation (21) (with o, 
for o,) holds with (77) subtracted from the coefficient of A dt/do,. This 
equation, in which t—2/A is practically —to,, could be treated approxi- 
mately, but this will not be attempted here since the physical nature of 
the process is already clear, and the purposes of the paper have all been 
accomplished. 

Some idea of the time scale of the final energy transfer is however 
desirable. The pulse (62) will still be considerable when o, or oy is, say, 
half its initial value, so that by (23) or (58) ¢ is of order max | €|/ Aoy t/a. 
Its reflection in the piston catches up with the rear shock when, at it, 
4(t—a/A), 1. e. —4tog, is of this order, so that by (58) 


eee i | ee. (78) 
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Thus, when the transfer of energy from the rear pressure wave to the 
N-wave is in full swing, the rear shock strength is of order the square of 
its initial value, and the time ¢ is of order the duration of motion of the 
piston f;, divided by the cube of the initial shock strength. 

The ultimate asymptotic behaviour of o, (which measures the rear 
shock strength, and also its position, at time t) may also be obtained 
without difficulty. After a long time, when disturbances have become 
very small, only the Pv term in the specific energy affects the integrated 
energy in the N-wave, and this is therefore (PPy/VPyy)(o{—o3)At. 
Equating this expression—with of At replaced by its limiting form 


| @o+12d0%) dé 
0 
as given by (23)—to expression (72) gives that (neglecting 3) 
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The Friedrichs theory gives the limit as (79) without its last term. Thus. 
that theory describes the final stage of decay of the rear shock correctly 
only to a first approximation, although it describes that of the front: 
shock correctly to a second approximation. 
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XCIX. Energy Losses of Ionizing Particles at Relativistic 
Velocities in the Photographic Plate. 


By H. Mgsset and D. M. Ritson, 
Dublin Institute for Advanced Studies*. 


[Received June, 26, 1950.] 


‘THE theory of ionization losses at relativistic velocities in condensed media 
was first treated quantitatively by Fermi (1940). Fermi showed the 
importance of considering the dielectric properties of the medium, and his 
calculations were later extended by Halpern and Hall (1948) to cover a 
rather general model for the dielectric properties. This theory is based on 
well established laws of molecular physics and classical electrodynamics, 
and the presence of any large deviations from it would be most surprising. 
In recent investigations Halpern and Hall (1948), Bradt, Kaplan and Peters 
(1950,) Occhialini (1950), with the photographic plate on the grain density 
of relativistic tracks, no increases of grain density greater than 10 per cent 
have been observed up to momenta of the order of 104 Mc. These obser- 
vations have been interpreted as due either to a failure of the Halpern—Hall 
theory or to the hypothesis that ionization caused by distant collisions will 
not contribute to the development of Silver Bromide grains Occhialini 
(1950), Powell (1950), Wilson (1950). Both of these assumptions are 
difficult to accept. 

However, detailed examination of the problem shows that these inter- 
pretations are based on an incorrect application of the results given by 
Halpern and Hall. Halpern and Hall were interested in the energy losses 
of particles traversing large thicknesses of matter and therefore calculated 
the total loss of energy in all collision processes. 

The quantity required for calculations on the photographic plate as also 
for certain applications in cloud chamber measurements is the “ probable 
ionization’. The development of a Silver Bromide grain will be a function 
of the energy actually converted into ionization losses in the grain, and will 
not be a function of the total energy loss. 45-rays with ranges very much 
greater than the radius of the grains will lose only a small fraction of 
their energy in the grain of Silver Bromide, and thus the energy losses, in 
the form of such 5-rays, of an ionizing particle will not contribute to the 
development of the Silver Bromide grain. The quantity required is the 
energy lost in all collision processes involving energies less than some 
maximum energy value E,,,,. This quantity is defined as the “ probable 
ionization ”’. 

The Bethe—Bloch formula for ‘‘ probable ionization ”’ losses is as follows : 

dW . 2anet Mv" nas 5 
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where all the symbols have their usual meaning. 


* Communicated by Professor L. Janossy. 
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For large values of the momentum the theory of Halpern.and Hall shows 
the Bethe—Bloch theory overestimates the energy loss by an amount 4 
given by 

2ane* ne*h* 
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thus giving an asymptotic form for the “ probable “ energy loss for high 
momenta of 
_ dW iss 2anet { THC as +1}, 
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This result is identical with that given by Fermi(1940) and is quite indepen- 
dent of any of the assumptions as to the molecular structure of the medium. 

For low momenta in the absence of conduction electrons the Halpern and 
Hall correction is small and thus the ratio of the * probable ionization ” . 
losses is insensitive to-the minimum value for the “‘ probable ionization ” 
losses to any assumptions made about the molecular structure. 

On detailed examination of the Halpern and Hall theory we find that the 
relativistic increase is entirely associated with the term related to “Cerenkov 
Radiation ’’. We have investigated the possibility that energy loss in the 
form of Cerenkov Radiation will be able to escape from the crystals. 
This contribution is given by the second term in the formula (28) of the 
Halpern and Hall paper. The term is 
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6 is a function which is 0 for (1—f%«)>-0 and z for (l1—f*e)<0. Ona 
multifrequency theory for a dielectric « will be large in the neighbourhood 
of a resonance frequency and less than unity or small elsewhere. Thus the 
contributions to the Cerenkov Radiation are mainly from frequencies in the 
neighbourhood of resonance lines. There will be some loss of ‘‘ Cerenkov 
Radiation ” particularly for frequencies associated with the K and L 
levels for which the mass absorption coefficient will be less than the order of 
2000 g./em.? necessary for stopping the radiation from leaving the crystals 
of Silver Bromide. The theory is thus only approximate. 

Fig. 1 shows a graph of the ratio of “* probable ionization ’ to minimum 
‘“ probable ”’ ionization. For comparison purposes the dotted line shows 
the ratio of total ionization to minimum total ionization. The difference 
is due entirely to the. contribution from high energy d-rays. E,,,, Was 
taken as 5 kV. corresponding to a range of 0-2u in Silver Bromide. We 
used the results of Halpern and Hall, calculated for iron. This is a fair 
approximation as we were concerned only with the ratios of the ionization. 
We note that the total relativistic increase is only 15 per cent. The reason 
for the difference between this result and that actually observed in a cloud 
chamber gas is due quite simply to the fact that the onset of the polarization. 
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effects occur at much lower momenta due to the dielectric properties of the 
medium. In air we should expect the polarization effects to only assume 
importance at momenta of the order of one hundred times greater than those 
for solid Silver Bromide. Thus we should expect an increase of probable 
ionization logarithmically up to the order of 1-7, and beyond that the 
jonization should become constant. 


Fig. 1: 
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Fig. 2 shows the histogram of the ionization versus numbers calculated 
for the sea-level spectrum of Cosmic Rays given by Wilson (1946). We 
note that it is a broad distribution from | to 1-15 I,,;, and the mid point is at 
1:07 Lam. The usual method of determining minimum ionization is to plot 
such a histogram and determine the mid point, and take this as corre- 
sponding to minimum ionization (Brown, Camerini, Fowler, Muirhead, 
Powell and Ritson 1949). Such a method, as can be seen by reference to 
fig. 11, will however determine a value of “minimum ionization” of the 
order of 1.07 1, The results for high altitudes are not likely to be very 
different. The energy spectra will contain more low energy mesons near 
minimum ionization, but this will be counterbalanced by more electrons 
with ionizations near 1-15 I, in- 


1132 On the Energy Losses of Ionizing Particles at Relativistic Velocities 


With measurements as made at present the maximum possible rise in 
ionization thus comes out at about 7 per cent above the experimentally 
determined minimum. 

Present available data gives no decisive evidence either for or against the 
above relativistic increase. It would however be important if it could be 
established with certainty. It would then provide discrimination between 
mesons and protons in the 10%V. range in measurements with the 
photographic plate. 


We are indebted to Professor L. Janossy and Dr. J. G. Wilson for 
discussions on this problem. 


Number of Mesons 
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C. A Theory of the Plastic Bulging of a Metal Diaphragm by Lateral 
Pressure. 


By R. Hitz, M.A., Ph.D., 
Metal Flow Research Laboratory, Sheffield*. 


[Received July 27, 1950.] 


SUMMARY. 


Explicit formule are obtained for the stresses in a metal diaphragm 
which is bulged plastically by lateral pressure. The predicted influence 
of work-hardening on the shape of the profile, and on the relation between 
polar strain and curvature, agrees well with experimental data. A simple 
expression is developed for the instability strain. 


§1. [yTRoDUCTION. 


THe ductility of sheet metal under biaxial stress is often examined by 
means of the so-called “‘ bulge test’. A uniform plane sheet is placed 
over a die with a circular or elliptical aperture and is firmly clamped 
around the perimeter. An increasing hydrostatic pressure is applied 
to one side of the diaphragm, causing it to bulge through the aperture. 
From the measured profile and thickness of the plastically-deformed 
diaphragm near the pole it is possible to calculate the local state of stress 
in terms of the applied pressure. If, in addition, the state of strain is 
measured by means of a scribed grid, the stress-strain characteristics of 
the metal under biaxial tension are obtained. An especial advantage of 
this test over any other simple one is that a greater range of pre-instability 
strain can be attained. 

To understand the significance of this test it is desirable to have a 
theory which explains how the varying shape of the bulge, and the 
distribution of strain, depend on the properties of the metal. The theory 
should also be able to predict the moment when the process becomes 
unstable. The only published work of importance appears to be that 
of Gleyzal (1948) whose analysis is, however, restricted to small strains. 
Gleyzal formulated the necessary differential equations for a circular 
aperture and integrated them numerically in a particular case, but did 
not obtain a general solution. Moreover, he adopted physically 
unacceptable stress-strain relations’ in which the components of total 
strain are made to dependsuniquely on the current stress. The present: 
investigation is an attempt to obtain an explicit solution, valid for large 
deformations and based on the Lévy—Mises stress-strain equations. 


* Communicated by the Author. 
SER. 7, VOL. 41, NO. 322.—NOV. 1950 AI 
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§2. Bastc EQuaTIONS. 


Consider a circular diaphragm of radius a and thickness ¢), where t/a 
is so small that bending and shearing stresses can be disregarded. Let h be 
the polar height of the bulge when the applied pressure is p. The bulge 
being a surface of revolution, we can take cylindrical coordinates (r, 2) 
for a particle, referred to axes in and perpendicular to the original plane 
(fig. 1). The initial radius to a particle is denoted by s, the current local 
thickness by ¢, and the inclination of the surface to the plane of the 
edge by @. 

By well-known geometrical formule the principal radii of curvature are 
dr 

de ? (1 ) 
where p, refers to a section such as fig. 1, and p, to a perpendicular section 
containing the surface normal. The radii are not independent since the 
elimination of @ from (1) leads to 


Pe=r cosec 6,  p,=sec 0 


Tipe 1 Fa ae 


The radii are equal at the pole, but generally nowhere else unless the bulge 
happens to be a spherical cap. 


Let o, and o, be the principal components of stress in the surface (the 
suffix ¢ denotes the circumferential or hoop direction, and the suffix ¢ 
denotes the tangential direction in the plane of fig. 1). The equations of 
equilibrium for a thin membrane in directions normal and tangential to 
the surface are respectively 


(ric, ta.) ie I Pee 


po P 
= 5B, (Pe) = 5 (2 ee Pes 


ae A iS af ma) 
oy 5) Pe = 57 cosec 0. 

The formula for o, may also be obtained by considering the equilibrium 

of any cap of the bulge. The stresses are expressed by (4) directly in 


terms of the pressure and the buige geometry. Notice that if p,2 
then o,So). pew 
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The yield criterion is taken to be 


: (Gee Os =ro(e), ns. peaeeeen wen sgt Paice 5) 
where o(c) is the stress-strain relation in uniaxial tension or compression, 
« is the equivalent logarithmic strain (to be defined later), and fis a 
homogeneous function of the stresses which does not vary with the 
amount of deformation. 

We neglect elastic strains and adopt the Lévy—Mises relations 
€ 2a 


a ee 


oes ee 6 
€é, 20;—0,’ ie) 
where ec, and ¢, are the principal strain-rates in the surface. The condition 
for incompressibility is 


eS ck fas | 
Sec ia cv: : . . . . . . . (7) 


If u=r—s is the radial displacement of a particle, and v=a-is its radial 
speed, then 
ee em OU : 
pare arr @ tan 6. © Waal toe’ 2) eh (8) 
For convenience, rates of change are measured with respect to the 
height h, which is a suitable monotonically varying quantity. The 
derivative following a particle is then 

7] ] 
oh 80) dr ’ 


so that, for example, 


dene _ ou //,_ au 
v= (5,405) ¥ or = 5 ( eg 


The boundary conditions on the clamped edge are 


U, Vv, 2=0 when r=a, for all h. tee er) 


In practice the die is provided with a rounded edge to avoid premature 
fracture there, and the sheet is clamped around a circle of radius somewhat 
greater than a. The assumed boundary condition is therefore not strictly 
satisfied, but, when the radius of the rounded die is less than about a/10, 
the discrepancy between theory and experiment due to this cause appears 
to be negligible. 

The condition (9) implies that the hoop strain vanishes at the edge. 
We shall see later that the tangential strain may or may not be zero at 
the edge, depending on the particular material. When the metal is 
thoroughly annealed, with a vanishingly small yield-stress, it is evident 
that the edge must strain and work-harden in order that a sufficient 
tangential stress is available to support the thrust of the pressure over 
the whole bulge. Hence, according to (6), o,=}0, at the edge. : The 
edge may, however, cease to deform later when it has become sufficiently 


ALZ 
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hardened in relation to the rest of the diaphragm. On the other hand, 
when the metal is cold-worked and has a sharp yield-point, o, is not 
necessarily equal to 4c, at the edge, and deformation there is then 
impossible (elastic strains being disregarded). In particular, when plastic 
bulging is about to begin, both principal stresses are everywhere equal to. 
the initial yield stress. The absence of peripheral strain has been observed. 
in hard copper by Brown and Sachs (1948). 


§3. A SPECIAL SOLUTION. 


Before trying the general problem it is worthwhile examining whether 
there exists a particularly simple solution for some special stress-strain 
curve. Let us see what follows from the trial hypothesis that the particles. 
describe circular paths which are, moreover, orthogonal to the momentary 


Fig. 2. 


0! 


profile (fig. 2). Let R be the radius of the path of the particle initially 
at radius s. The conventional circumferential and tangential strains are 


€=—(1— cos @) (: z) 


€.—=(1— cos @) Cron 
s 


‘s 2 
Assuming further that e,=e,, or 


dR R 
re aes 0; 
a2 
we find that R=} (+ -s) eM an T irri ih La 


since R=0 at the edge for zero hoop strain. This equation, being equiva- 
lent to s(2R+s)=a?, is the condition that the particle paths are arcs of 
coaxal circles having the edges as limiting points. The successive profiles 
of the diaphragm are therefore the orthogonal set of circular ares passing 


through the limiting points. It follows that the diaphragm is deformed 
into a spherical cap of radius 


ar 
p=3(F +h). 5. a) mie a eh Ute ae (11) 
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Since the principal strains in the surface are equal, the distortion of 
every element of the sheet is a uniaxial compression normal to the surface. 
The equivalent logarithmic strain ¢ is therefore In (t/t). Moreover, 
whatever the stress-strain relations, and provided only that the material 
is isotropic, o,=o,. That is, the state of stress in every element is a 
balanced biaxial tension. This is equivalent to a uniaxial compressive 


Stress if the hydrostatic component of stress is without influence on yielding. 
Hence (5) must reduce to 


o,=9,=—a(e), where «=In (t,/t). 


From the first equation of (3), p=2ot/p and hence the product has the 
same value throughout the diaphragm at each stage. This demands that 
the stress-strain curve is 


CET ONC, ae A ry (ban Metsties. nis te a(t) 


Fig. 3. 


where a, is the initial yield stress. The second equation of (3) is now 
satisfied identically, and the formula for the pressure becomes 


P= 2a lole=—40 ol h(a) Lo ee (18) 


The pressure reaches its maximum when the bulge is a hemisphere. Our 
initial hypotheses have thus been proved consistent and we have obtained 
the solution for a special material. Its stress-strain curve is, however, 
unlike that of any actual metal in that the rate of hardening rises with 
increasing strain. In fact, the material extends under a constant load. 
None the less, the relation between the polar strain and curvature is in 
good agreement with that measured for various metals, as will now be seen. 

We first obtain an expression for the strain at any point in terms of the 
polar height. Consider any two circles with centres O and Q (fig. 3) and 
any line through one of their points of intersection P meeting the circles 
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in Z and X. Let ZOY and XQP, be the respective diameters through 
Zand X. Then the points P, Py, and Y are collinear since angles ZPY: 
and XPP, are both right-angles. Let Z, be the meet of the diameters 
XP, and ZY. Then, if the circles cut orthogonally, so do these diameters. 
For angle XZ,Y is the sum of angles OZP and QXP, and therefore of 
angles OPZ and QPX ; the sum of the latter is a right-angle since the 
radii OP and QP cut orthogonally, by hypothesis. Translating these 
results in terms of the geometry of the diaphragm, we see that if Z is the 
pole of the bulge and P any particle (Z, and P, being their initial positions), 
then PP, produced passes through Y, the other end of the diameter 
through Z. To find the final position of a particle initially at Po, therefore, 
it is only necessary to produce YP, to cut the profile of the diaphragm. 

If P’ is the foot of the perpendicular from P to OZ, the conventional 
hoop strain is 


ted a NAL pate Lea 
TiS SE BG Ae NA: 
by similar triangles. Hence, since Z,P’=z and YZ,)=a’/h, 
h 
a= = (14) 


The relation between the curvature and the through-thickness logarithmic 
strain at the pole can therefore be expressed, parametrically in terms of 
the polar height, as 
h? p h?\ /2h 

e=2in(1 +5), £=(1+5) = lace 15) 
Brown and Thompson (1949) found that the relations between polar 
curvature and strain for annealed stainless steel, annealed low-carbon 
steel, and certain aluminium alloys (both hard and soft), differ only 
slightly. The relation (15) is in excellent agreement with their mean 
curve up to a strain of order 0-4. It is only at greater strains that a 
discrepancy occurs, due partly to the occurrence of an instability and 
partly to the increasing dissimilarity of the stress-strain curve (12) from 
any actual one. The reason for the insensitivity of the (p/a, ¢) relation to 
variations in the properties of the metal will become clear when the 
general solution is considered. 


§4. THe INsTapiniry STRAIN. 


Brown and Sachs found that when the pressure reaches its maximum 
a secondary bulge begins to form over about the top half of the diaphragm, 
and that this continues under a falling pressure. In this sense the pressure 
maximum can be regarded as the moment when the process becomes 
unstable. We can calculate the corresponding instability strain by means 
of (15), which we have seen to be a good approximation for a wide variety 
of metals. The relation between the pressure and the polar curvature is 


‘p=2ot/p where ¢«=In(t,/t), 
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since the local stress-state is a balanced biaxial tension equal to the 
yield stress c. The pressure is a maximum when 


ldo 1 dp 
dat ae sys de? (16) 
or, from (15), when 
ida 3. p 
ode 2 2h 
On expanding in powers of ¢, this reduces to 
1 do al 1 1 
= de —_— iy a De > . . . . . . . ( 1) 


with an error of less than 1 per cent in the relevant range. If the stress— 
strain curve can be fitted by a power law c=Ae”, the polar strain at. 
maximum pressure is 


em = = (2n-+1). x teen erect 3) 
For greatest accuracy the exponent » should be chosen to give best fit 
to the actual stress-strain curve in the neighbourhood of the instability 
strain. 

In uniaxial tension the strain at maximum load is equal to n; in a 
thin spherical shell expanded by internal pressure, to 27/3 ; and in a plane 
sheet under balanced biaxial tension, to 2n. In all these tests the 
instability strain is vanishingly small for a heavily cold-worked metal, 
but is equal to 4/11 in the bulge test. This strikingly demonstrates the 
superiority of the latter as a means of obtaining stress—strain characteristics. 

Lankford and Saibel (1947) have attempted to calculate the instability 
strain in the bulge test by assuming (without previous analysis) that the 
bulge is spherical and that the strain is given by (14), but with logarithmic 
strain instead of conventional strain. This therefore agrees with the 
present analysis only when the deformation is small, which is certainly 
not so at maximum pressure. Lankford and Saibel’s final result is still 
further vitiated by a gross error committed in taking partial derivatives. 


§ 5. GENERAL SOLUTION FOR COLD-WORKED METALS. 


The stress-strain curve of a metal in a fairly severely cold-worked 
state can be approximated by a linear relation 


o/Y=1+He, Si es ghee ere LL) 


where Y is the initial yield stress and H is a non-dimensional constant 
of order unity. We look for a solution giving the various unknown 
quantities as power series in h/a. We shall take the solution only as far 
as the second approximation, and it will become apparent that the yield 
criterion is then required to be satisfied only in first approximation. 
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Any yield criterion (5), which reduces to von Mises’ criterion for stress- 
states near to a balanced biaxial tension, can then be written as 


(opto, )/ 2 Ve Ia Hey 3. ee ae ee (20) 


where the equivalent strain ¢ is equal to lu (t,/t) to first order. 

Elastic strains being neglected, bulging does not begin until the whole 
sheet is stressed to the yield point. Initially, therefore, o,+-0,=—2Y 
everywhere. The only possible stress components satisfying this 
condition and the second of equations (3) (with t=f)) are 


C= 0;==Y. 
This is the first approximation to the stresses. Substitution in the 
second of equations (4) gives 6=pr/2t,)Y to first order. Now 


it tan @. 


dr 


Hence, after substituting for 0 and integrating 
z= (1 3). ee te 


ar az 


2 


on introducing the necessary conditions that z=h when r=0, and z=0 
when r=a. From (6) and (8) after substituting 


(bet 


mee 
we find that (to first approximation) 


2hr ( 5) 
== Por 1 St Si ) 
a a 


since v=0 when r=a. Finally, from (7) 


) 


ie aed os li 
Ines aes 
n ; Z (1 a): 
‘This completes the solution to first approximation. To this order the 
bulge is spherical (or rather a paraboloid), irrespective of the rate of 


work-hardening. 


The same sequence of solution is used for the second approximation. 
The yield criterion (20) becomes 


oto, |  2Hh Pr 
BN: ait ae (: ‘ 


Combining this with the second of (3) we find that the second approxima- 
tion to the stresses is 
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We see that o,20, according as H21. Continuing as before, we find 
successively that . 


2hr ee h? | 
d=, 1+ {ay 5— sy} S], PTO CL)) 
i-(1- a) +t SE], Bails AE ous etlgn 4. (23) 
p  A4ht : h? 
ee 14-43H—5) 5 I, ere (24) 


bo 


bp 


= (1 —B)[1+4{H—s) 5 ey Ti “8. (25) 


wu hr 1 re t h 

aaa) EE ORG], OH 

Fon he ; : 
a ee ee 


‘The polar radius of curvature is given by 
a 2h h 
i 4. be 
: =F). HH+1) 5 |. ea (28) 
The equation of a circle passing through the edges and the pole is 
(to this approximation) 


, 2) 2 fy 
3- (1-3) (1+ G)- Sabla Moccia AOL) 


On comparing (23) and (29), it is found that the maximum relative 
difference between the bulge profile and the circle is 


es 2 
Max. (- = =KH-1)5, when r=—— 


h /2° 

This vanishes when H=1, in agreement with the solution for the special 
stress-strain curve (12), since expe~l+e to this approximation. 
When H>1 (which is so for all save exceptional degrees of cold-work) 
the profile of the bulge lies above the circle (29), and below when H <1. 
Expressed otherwise, the polar circle of curvature cuts the original plane 
outside the perimeter of the die when H> 1 and inside when H<1. This 
was found experimentally by Brown and co-workers at all stages up to 
the pressure maximum (H~$ for the hard copper used by Brown and 
Sachs and ~4 for the hard aluminium alloy used by Brown and 
Thompson). 

From (27) the polar strain is 


rl 


E= 


[1-48-15]. wee (0) 


Combining this with (28) we obtain 


e-| 5 Saale eens e131) 


Ge We 16 
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Thus, the effect of the rate of work-hardening on the (p/a, <) relation 
is comparatively slight up to strains of order 0-3, because of the square 
root dependence and the small contribution made by the term H/16. 
(When H=1, equation (31) is the second approximation to the solution 
of (15), as is to be expected.) The relation between bulge height and 
polar strain, or between bulge height and polar curvature, is much more 
dependent on the value of H. We see from (30) that, for a given height, 
the most severely worked metals (small H) have the largest polar strains 
but the smallest polar curvatures; this is in agreement with the 
observations of Brown and Thompson. 

An inspection of equation (28) indicates that a rough criterion for the 
range of validity of this solution (say to an accuracy of 1 per cent) is 


(H+1)< <}. 


Thus, when H=1, the solution is probably reliable up to a polar height 
of ta. When H=3, a value typical of moderately cold-worked metals, 
the solution is probably reliable up to a polar height of a/2,/2~0-35a. 
It is only when H is less than about 4 that the range of validity probably 
extends past the pressure maximum which, according to (24), occurs at a 


polar height of 
9 1/2 
lan | "s 


Thus, the solution does not appear capable of explaining the secondary 
bulge observed by Brown and co-workers ; this may indeed be peculiar 
to annealed metals with a well-rounded stress-strain curve. It has not 
yet been found possible to obtain any simple approximate solution for 
such metals. 
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SUMMARY. 


The magnetic vector potential distribution in axially symmetrical 
systems containing current carrying conductors and iron of arbitrary 
shape can be obtained experimentally with good accuracy using resistor 
networks. Two networks are possible, one for plotting the flux function 
rA, and one for plotting the vector potential A, itself. The latter one is a 
modification of a field plotting network described previously. Some 
confirming measurements on a model of an infinitely long coil were 
carried out. | 


§1. Lyrropuction. 


THE determination of the magnetic field due to currents flowing in 
conductors requires the mapping of the magnetic vector potential. In 
the presence of iron a mathematical analysis becomes almost impossible, 
but it was shown by Pérés and Malavard (1938) and by Peierls (1946) that 
the problem can be solved experimentally with the help of the electrolytic 
tank analogue in those cases where all currents flow in straight parallel 
conductors. The problem is then a two-dimensional one, and only one 
component (¢.g. the z-component) of the vector potential exists, which 
satisfies a Poisson equation in z and y. In the electrolytic tank analogue 
the iron circuits are represented by wax models ; the current conductors, 
supposed to be perpendicular to the surface of the electrolyte, are repre- 
sented by a number of probes dipping into the electrolyte, and feeding 
the appropriate currents into the tank. 

This model representation cannot be applied to rotationally sym- 
metrical problems, where all currents flow in circles concentric with the 
axis of the system. The vector potential has then only an angular 
component A, which has to satisfy this equation : 


GC Aveo Agel CAs ae Ul 4 
Tea oe a Ate eee 
ae! OF ro Ay 10 9 o 
where J is the current density in the conductors, in amps/cm.?. 
Writing | raj, ead: Ole eer eer. 0 (2) 


* Communicated by the Author. 
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this equation can be put into the form 


Gla ee te Ge awd AAG Ss 4dr 
aii. aaa ap ae Oe rime GLO 6 (3) 
It was recently shown by Peierls and Skyrme (1949) that an experi- 
mental solution for this function Y can be obtained by using two tank 
measurements, first setting up the model in a tank with equal height of 
electrolyte, and then in a tilted or wedge tank. The difference between 
the two measurements can be used as a correction for the results obtained 
in the first tank, provided the correction is small, 7.e. measurements 
are taken not too near the axis. 


§2. APPLICATION OF Resistor NeTworK METHOD. 


A direct measurement either of the vector potential A,, satisfying 
equation (1), or of the function (rA,), satisfying equation (3), can be 
obtained with a resistor network similar to the type described recently 
by the present author. 

Consider the equation 


Gu oF 


ae + ay =G(F xX, y). a Sohite tiet etn (4) 


In the earlier paper (Liebmann 1950) it had been shown that the difference 
equation corresponding to equation (4) is satisfied by a network of inter- 
connected resistances of equal value Ry. The values of the function F 
are represented by the voltages at the resistor junctions, and the function 
G is related to the currents I, fed into the junctions from an external 
source, through the equation 


G=--IR (Agee ot ot ne 
Now let us apply the transformation 
di Vipins L—=2/r 
yi-k + kA, == 15 (hss. Pree See) 
UE res y=logr 


where k is any real constant. This transforms the functions F(a, y) and 
G(F, x, y) into f(r, z) and g(f, r, z) and equation (4) into 


arf aap ape a0 
dz? | Or 7eore Wf, 7; 2), 2 1G OS io oe 
or with 
Teh geags Pe a 
into a. 
jae of lcoy 
bat Bett oe 0 2 
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if we now choose the four resistances R,....R,, joined in the network 
junction P, (see fig. 1) such that 
R= Re=7- "Ro, | 
R= (ry! *#—1_!-*)Ry/(1—h), reer FT 10} 
R,=1,)_*—ro**)Ro/(1—k), 
Fig. 1 
P4 
R4 
R 
Po 
R3 
P2 
Fig. 2 


R, being an apparatus constant, we obtain a network which solves the 
differential equation (7) to a higher degree of approximation than the 
difference equation which would replace equation (7) in numerical methods 
of solution fF. 


+ Except for k 0, the difference equation as used in most numerical techniques 
is only correct to second order terms, whereas the specified network is correct to 
third order terms, whether k->0 or not. This point is mostly not very important 
in regions distant from the axis, but cannot be overlooked as the axis is 
approached. 
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In the previous paper (Liebmann 1950) it was shown, through the 
application of Kirchhoff’s law, that if the function F is represented by 
voltages, then the function G is represented by currents fed into the 
network points P, from an external source (see fig. 2). Similarly, we 
find that the function g* is related to the currents I fed in at Py by 


Tee —g7(Az)Ri) “- oss 2 
where R, is given by equation (10). 


§3. Resistor Network FOR FINDING (rA,). 


The two alternative equations (1) and (3) for the vector potential A, and 
for (rA,) can therefore be solved by constructing the appropriate networks. 
Taking first equation (3), we have k=—1, and 


R,=R,=rR», 
Ro=(79' —127)Ro/2, reheat 
R3=(r4’—197)Ry/2. J 


The currents to be fed in are then 
4 
lb=+7, J(Azy Res 2k ee 


‘Obviously, current would only be fed in at the positions occupied by the 
conductors. If iron is present, the network must be open-circuited along 
the boundaries of the iron, the network being correctly terminated there, as 


discussed in the earlier paper, in order to fulfil the boundary condition 
O(r Ag) 


on 
§4. Ruststor NETWORK FOR Frnpine Ay. 


It will mostly be more attractive to use a network which solves directly 
for Aj. Then k=-+-1, and the previously described network, representing 
the Laplacian operator for axially symmetrical fields, can be applied. The 
4 1 
Tou = y) Az?/R,; this 
would require the use of the iteration method described before (Liebmann 
1949). This need is eliminated by the following modification: Each 
network junction is connected through an external resistance of value 
R;=rR, to a point of zero potential (fig. 3). | 

We then have these values of resistances : 


R,=R,=R,/"9, 
R,=Ry log(79/73), 


currents to be fed-in would then have to be 


R,=Ry log(74/"9, CS 
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This network now represents the operator 


Cao LO 1 
ie ae =| 
which differs from zero only at the position of current carrying conductors. 
at these positions, the current 


ot 4dr 
T= 79 J (Az/Ri= 7p Ir(4z)/Ro - Par re ey eee ORs) 


is fed into the network. It is thus possible to measure the vector potential 
distribution inside current carrying regions within the limitations imposed 
by the finite size of the interval. Iron circuits are again represented by 
open circuiting the appropriate boundaries, as dA,/dn—0, for infinite 
permeability; the values of the resistances R; have to be doubled along the 
boundaries and quadrupled at interior corners to maintain correct relations 
at the boundaries. If it is desired to represent correctly a finite value 


of permeability, the network constants have to be modified, by increasing 
the value of the constant R, in the ratio of (firon /Hair), I all resistances 
within the area occupied by the iron, or by tapping down the voltage in the 
ratio of 1:~ at the boundary points, using the resistances R; as potential 
dividers, in order to satisfy the boundary condition Ag;,,, =1/mu Agair- AS 
previously described, curved boundaries can be simulated by shunting the 
network resistances ; the values of the resistors R, have then to be adjusted 
in a similar way, as mentioned before when referring to the boundaries of 
iron circuits, taking the change of size of the effective area of the network 
element (4zx Ar) into account. 

Ifa problem to be solved involves regions which are open in the direction 
away from the axis, then it will be necessary to follow the technique 
applied in similar cases of electrostatic fields. A small scale model is set 
up on the network, in which the distance of the current carrying conductors 
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from the axis is small compared with the radial extension of the network, 
and the approximate values of the vector potential at a suitable radial 
distance r, are determined. The model is then enlarged so that this. 
relative radial value rg now corresponds to the edge of the network, and 
the potential values required at this distance from the axis are forced on to. 
the model by feeding in additional currents at these points. Alternatively, 
a suitable flux line may be made an equipotential for the voltage. The 
axis is, of course, an equipotential line of zero potential. 

Instead of determining the values of A, at the mesh points, one can 
plot out directly the lines of equal A, using the ganged potentiometer 
arrangement described before (Liebmann 1950). 


§5. PracticaL REALIZATION AND TESTS OF ACCURACY. 


The correctness of the discussed method was confirmed by several 
experiments with a linear resistance strip. This strip (fig. 4) represented 
one section of a precision (r, z) resistance network being built in this 
laboratory. The gradation of the resistance values was on a logarithmic 


Fig. 4. 
416K 1396 Bil 1026 446 367 575 446 Ohms. 


basis (equation (14)), and the values were the same as those employed in the 
previously described resistance network (R j=2000 ohms, except 
that the resistance values were kept within about +0-2 per cent of 
their nominal values, and that the termination strip was extended to go — 
out from r=20, the edge of the network proper, to 7=50 by a few large 
steps. Each junction, at distance r from the axis, was connected to the 
axis (V0) by a resistor of value R;=rRp, according to equation (14). The 
finite size of the interval makes the nominal value rR, incorrect by a small 
factor, of the order of 1-8 per cent at r=1 and 1-2 per cent at r=2, the 
correction factor diminishing to a negligible amount for greater values of r. 
Expressing the length r pucugh its equivalent resistance value leads to the 


corrected value Ri=4 R22’ (1/R;); these corrected values were used in 


Mi 
the finalexperiments. The linear resistance strip can be used to investigate 
problems in which A, does not depend on z, e.g. the case of infinitely long 
coils of constant radii. In one such experiment a single layer coil of ro=aa 
units was assumed, the potential being set to the relative value A,=1 at 
ry=10. Table I. gives the experimentally found values and the theoretical 
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values A,—Ajr/ry. The agreement between experimental values and 
theoretical values is very good, particularly in view of the small number of 
components involved, which does not bring the averaging properties of the 
network fully into play. 


TABLE [, 
Aj(r) exp. 
} 
rir Before correction, After ce A(r) cale. 
R5=rRy. Ri=4R32(1/R,) 
1 
0 0 0 0 
0-1 0-1002 0-1000 0-:1000 
0-2 0-2012 0-2004 0-2000 
0-3 0-3012 0-3004 0-3000 
0-4 0-4011 0-4004 0-4000 
0:5 0:5009 0-5002 0-5000 
0-6 0-6007 0-6002 0-6000 
0-7 0-7005 0-7001 0-7000 
0:8 0-8003 0-8001 0-8000 
0-9 0-9002 0-9001 0-9000 
1 1 i} 1 


— exact 
© 7o= Junits 
x Fo=Gunits 
+ "o =(Zunits 


The vector potential A, falls away as 1/r outside an infinitely long single 
layer coil, end an incorrect field distribution will be obtained at some 
~ distance from the coil unless the model is made very small in relation to the 
network, giving low “ resolution ”’ near the coil itself. This difficulty can 
be eliminated by first measuring the field distribution of a small scale 
model in which the potential A,=0 at r=50 units, and by then setting up 
the full size model and impressing the potential values found in the first 
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test on the new boundary points. To study the quantitative side of this: 
procedure the field distributions outside single layer coils of radii ro=3, 
ro=6 and r)=12 mesh units were investigated. The results are plotted in 
fig. 5, where the full curve represents the theoretical curve Ag(r)=Apro/?- 
For the smallest coil, good agreement between theoretical and experimental 
values (with A,=0 at r=50) is found up to about 7/7922-5. The error 
is here —0-004 Ay. From graphs of the kind given in fig. 5 information 
about the appropriate size of the small scale model for the first test can be 


obtained. 
Fig. 6. 


lron shroud at r-12 


[ron shroud at r=13 


The effect of an iron shroud on a coil of inner radius r;—8 cm. and outer 
radius ry=12 cm. is shown in fig. 6. Curve 1 gives a plot of A,(r) for a coil 
without iron being present. The boundary at r—20 was set to the value 
obtained from a quarter scale model as described before (circled points). 
The total current fed into the network was 0-2500 V, ma, where V, is the 
feeding voltage of the network *), or J=6-25 x 10-5 V, amps/cm.2. 

Within the coil, the measurements give Aj,=0-02508 rV, gilberts, hence 
H,=1/r O(rA,)/0r=0-0501, Vo oersteds. Calculating H, direct from the 
formula H,=47/10 AT/cm. we obtain 0-0500 V, oersteds; this appears to be 
a quite satisfactory agreement with the experimental results +. 


i ee 
* All measurements of voltages and currents were referred to resistance 
measurements, as this gives a higher accuracy and greater convenience; the 
feeding voltage V, appears then as a common multiplication factor. 
+ The continuations of the inner and outer branches of the A, curve, which 
obey the (kr) and (k/r) relationships, would intersect at r— 10-2. the centre of 
gravity of the coil windings. f 
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Placing an iron shroud (u= oo) tightly around the coil (at r=12) changes 
the shape of the A, curve within the coil considerably (curve 2, fig. 6), and 
increases the value of A, inside the coil to A,=0-04243 rV>, i.c. the iron 
shroud has increased the slope of the A, curve, and therefore the magnetic 
field within the coil, by a factor 1-69,. If the iron shroud is placed a small 
distance away from the coil, at r=13 (fig. 6, curve 3) the A, curve becomes 
modified, and the factor of increase in field strength in only 1-60). A finite 
permeability (w1000) was then simulated, which produced very slight 
modifications of the A, curve near its maximum; the change is too small 
to be represented in the graph. 


§ 6. CONCLUSION. 


It is seen from these tests that a satisfactory accuracy can be obtained. 
As the modifications to existing field plotting networks, satisfying Laplace’s 
equation, can be easily effected, the measurement of vector potential 
distributions becomes a very practicable proposition, even within current 
carrying regions and in the presence of iron. The method should therefore 
have applications not only in the design of betatron and synchrotron 
magnets, and in similar problems, but also in many ordinary problems of 
electrical machine design. The same principle can also be applied to the 
mapping of other vortex fields of rotational symmetry, when the vorticities. 
and positions of the vortex rings are given. 
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CII. Double Stars due to Cosmic Rays. 


Bye ee aL 
Imperial College, London*. 


SUMMARY. 


Cosmic ray stars are distributed at random over a photographic plate 
when distances exceeding 1 cm. are considered. For small distances of 
the order of 1 mm. the number of close pairs exceeds that expected on a 
random distribution. The effect is not found on plates exposed at sea- 
level and there is some evidence that it is more marked if dense matter 
surrounds the plates. Several mechanisms are suggested to account for 
the effect but none is free from objection. 


§1. INTRODUCTION. 


MEASUREMENTS were made (Li and Perkins 1948) on the surface 
distribution of 1230 stars on 54in.2 of 100, C2 plates exposed for 
forty days at 3650 m. altitude. The plates were divided into (1 x) in. 
areas to facilitate searching. It was found that the distribution of star 
frequencies between these areas conformed to the Poisson distribution, 
showing that the stars were randomly distributed over areas of these 
dimensions. 

At smaller distances, however, the stars showed a tendency to be in 
close proximity of each other more frequently than would be expected 
if they were distributed at random. (See also Leprince-Ringuet and 
Heidmann 1948.) 

Thus, for a random distribution, the number of close pairs of stars, 
separated by a projected distance in the plane of the emulsion <r, 
should be a 

oe NGG) ere 
OO no aim U “ab ’ 
where 


Ni =number of stars in the 7th area, 


a, b=sides of the areas=1 in. and (4) in. respectively, and fr, 18 
small compared with 6. 

The actual observed number of close pairs, P(r,), was compared with 
the corresponding Q(r,). It was found that P(r,) was always greater 
than P(r,) up to r,=0-03in., which was then the largest distance 
measured. Moreover, the ratio (P—Q)//Q of the observed to the expected 
deviations was always greater than one. Different mechanisms were 


suggested to account for this excess of close pairs, which was termed 
“double stars ”’. 


* Communicated by Sir George Thomson, F.R.S. 
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In the present paper, the results are given of further measurements, 


extended to larger distances r,, and made on plates exposed under various 
conditions. 


§2. EXPERIMENTAL RESULTS. 
1. Measurements were made on the following sets of plates : 


(a) Plates exposed at 3650 m. altitude without any absorbing material 
immediately above the plates. (A total of 1230 stars on 54 in.?.) 

(6) Plates exposed at sea-level with no absorbing material immediately 
above the plates. (140 stars on 27 in.2.) 

(c) Five boxes of plates exposed in a narrow crevasse (~10 cm. wide) 
in a glacier at 3650 m. altitude, at the following depths below the surface 
of ice: 0m., 1m: 2m., 3:5 m., and 5-5m. (1126 stars on 77-5 in.’). 

(d) Four boxes of plates exposed at 3650 m. altitude with the following 


thicknesses of lead above them : 0 cm., 5 cm., 15 cm., and 35 cm. (374 stars 
on 16 in.’). 


The plates were all 100, thick of the Ilford C2 type. For the sake of 
convenience, they will be referred to as the “‘ 3650”’, ‘“ sea-level’, 
“ice ’, and “ Pb ” plates, respectively. 

2. Measurements were extended to include all close pairs with a 
projected separation of r,,<(0:1)in. The formula in (1) for the expected 
number of close pairs from a random distribution is no longer valid, 
since 7, is comparable with 6 now. The following exact formula can be 
derived: 


Ni(Ni—1) ar,” jpn se esi Re 
2 ab 37ab 2rab 


Qr)=z 


which holds true for 0<r,, <b. 
In the following table, the values of Q, P, (4P—4Q)?/4Q, and 
(P—Q)/(2 Ni) are tabulated, using the data of all four sets of plates : 


ry, (1/100) in. 2 3 4 5 6 7 8 9 10 
P(rp) 46 107-5 182-5 280 239 505 640 791 936 
Q(rp) 37-6 82-2 159-1 230-6 329-2 444:1 573-1 717-7 876-5 


(AP —AQ)?/4Q 1-3 O47 0-0 9:5 IS 55 O-Ose n0-d08 POSS 1-2 
103 (P—Q)/(2'N2) OS 82 17-2 21:2 21:2 23-3 25-6 20-8 


The last row shows the excess of close pairs per 1000 stars. There seems 
to be a saturation after r,~(6/100) in. Whether this is real or whether it 
simply arises from the fact that Q(r,) increases as ~r,,> while (P—Q) 
increases more slowly with r, so that the effect gets relatively small at 
large distances and is lost in the statistical fluctuations, it is not possible 
to say. (Also (P—Q) must be equal to zero at 7,—a, since, by definition, 
Q, <a 2 Ni(Ni—1)/2=P, a). 
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By applying the x?-test, it is found that the probability that this 
excess is only the result of statistical fluctuations is small, ~0-5 per cent. 
Hence we can conclude that this excess does exist. 

3. By comparing results of the different sets of plates, we found that 

(a) There is no excess of P over Q in the “ sea-level ”’ plates, or else 
the effect is so small that it is outweighed by the statistical fluctuations. 


Tp, (1/100) in, 2 Sree 5 6 ‘i 8 oF Ttto 
10° (P Q)(ZNi) —29 72 221 —50—-129 64 21 —86 —7-2 


Tt is seen from the above figures that the values of (P—Q)/(2Nz) fluctuate 
between positive and negative. And the x?-test gives a probability of 
10 per cent for these values of (P—Q) to be due to statistical fluctuations. 

(b) If we group together (i) all the plates which were exposed at 3650 m. 
with no absorbing material immediately above them (7. e. the * 3650 ” 
plates, the first box of the ‘“‘ ice” plates, and the first box of the “ Pb” 
plates), and (ii) those which were exposed under a thickness of ice or lead 
(i.e. boxes 2 to 5 of the ‘“‘ice”’ plates and boxes 2 to 4 of the “ Pb” 
plates), we notice that the values of (P—Q)/(2Nz) for the latter group are 
greater than the corresponding ones of the former group : 


1p, (1/100) in. 2 3 4 i) 6 7 8 9 10 
108; (P—Q)/ -(ii) 44 125 14:0 20:0 30-0 24:4 38-6 37-7 38-7 
(2N¢) 9) 7) 2:5 976°8) 11-0 17:3) 19-9» 14:9: 16-0 2r-2oeiie 


But the probability that the difference is due to statistical fluctuations 
is 25 per cent, which is not small enough to enable us to establish the 
existence of a difference. 


§3. DISCUSSIONS. 


1. Possible Errors: Let us examine the possible errors from which an 
apparent excess (P—Q) may arise. 

(a) Edge effect—The damage of the emulsion near the edge of the 
plates and at the ruled lines of the areas reduces the effective area of each 
rectangle to less than half sq.in. by an unknown amount. This might 
affect the calculation of Q to a significant extent. To eliminate this, let 
us study only the stars inside an imaginary rectangle of 1 cm. by 2 cm. 
at the centre of each area. The coordinates of the stars enable us to 
exclude all stars outside this central rectangle, thus giving a true area 
of 2 cm.* without any error due to incorrect ruling and damaged emulsion 
at the edge. 

We took, as a random sample, the plates exposed under 0 m., 1 m., 
and 2m. of ice. There are 52 areas altogether. The distribution of 
stars in these 52 1-cm. by 2-cm. rectangles is as follows : 


n OC 20 Bindi ob 62a 88 Galo 
number of 1x2 cm.? areas oa 


with n stars LL O:-l 4s 6 6ee8S oh aie ee oe 
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Thus, considering the central 1 cm. by 2 cm. areas only, we have 
ENi(Ni—L/2=(1X3X24+4X4%3-4 22. £2 13X12+1 17x 16)/2 
=1508; 
uid 4 4(a-+-b)r,, =| 


ab 37rab ce 
Here a=2 cm. and b=1cm. By comparing these expected numbers of 


close pairs with the observed number inside the central areas, we get 
the values of (P—Q)/(2N2) : 


yp; (1/100) in. 2 3 4 5 6 8 9 10 
103 . (P—Q)/(ZNi) 8-0 23:0 21:0 15-4 29-8 23:0 11-3 12:8—141 


The values of (P—Q)/(X Nz) for the same 52 areas but considering the 
whole (1 x4) in. area are given below : 


43 12:3 15:7 16-1 30-2 248 28-2 26-8 13-0 


A comparison of the two sets of values shows that the “ Edge effect ” 
does not contribute much, if any, to the excess (P—Q). 

(b) Uneven Thickness of Emulsion. Uneven thickness of the emulsion 
has an effect of concentrating stars and thus giving rise to a spurious 
excess, (P—Q). But it can be shown that even for a variation of thickness 
of up to 20 per cent, the effect is still very small. 

First consider a simple case. Suppose in the area (1 x4) in., half of 
the emulsion has a thickness d, and the other half has a thickness d’, 
djd'=1-2. Let n and n’ be the numbers of stars in the two halves 
respectively. Then n/n’=d/d'=1-2. 

If the stars were randomly distributed, then the number of close pairs 
that we should observe in the two halves is 

Tmt ee nat 


BE ES area ear rea 


for small values of r,, s being ($) x ($) in’. 
On the assumption that the emulsion thickness is uniform, we calculate 
the expected number of close pairs in the whole area as 


Hence the spurious excess of P over Q due to non-uniform thickness 
of the emulsion is 


2 ike : 
P—Q= —= [2n?+2n"—(n-+n')?]= . (n—n')° 


and 
(P—Q)/(2Ni)=(P—Q)/(n-+n')=a71,? . (n—n')?/4s(n-+7'). 
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Inserting n/n’=1-2, we get 
ary? 0-04’ 


(P—Q)(EN)= TE SS 


Assuming a value of n’=6, say, we get 


%p, (1/100) in. 2 3 4 etc. 

103. (P—Q)/SNi) 0-14 0-31 0:56 ete. 
These numbers are negligible compared with the observed values of 
(P—Q)|(2Ni). If the emulsion consists of two unequal areas of different 
thicknesses, the spurious excess can be shown in a similar way to be 
smaller still. Hence we conclude that the contribution due to non-uniform 
thickness of the emulsion is negligible. 

(c) Stars missed owing to Inefficient Searching.—Inefficiency in searching 
might also give rise to a spurious excess of P over Q. For instance, if 
most of the stars missed are those which are at a large distance to any 
other star, then a spurious (P—Q) can result. 

The efficiency of searching depends on the time spent on an area. 
At the rate (1 in.2/day) with which the plates were searched, we found 
that about 20 per cent of the stars were missed. The “ice ” plates have 
been searched twice, first at the normal speed and then at half that speed. 
The overall efficiency thus obtained was 99 per cent. 

Assume the extreme case in which all the stars missed were those that. 
were at a distance >(1/10) in. from any other star, so that they would 
increase the values of (Q(r,) but not P(r,). Assume a conservative 
efficiency of searching, in the case of the “‘ ice ”’ plates, of only 98 per cent. 


n=Average number of stars/area=1126/155=7-26. 
n' = Average number of stars/area, corrected for inefficient 
searching = 7-26/0-98=7-41. 
n'(n'—1)/n(m—1) =1-025. 
Q'(r,)=corrected values of Q(r,,)=1-025 (Qr,). 

We have 
trp. (1/100) in. 2 08 get ge! epee ee ee 
P(r) 17 43-5 82 127 182:5 229 294 352 416 

tp 17-3. 38-5 67-5 103-9 147-8 198-2 255-3 318-7 388-1 
10°. (P—Q)/(2N7z) —-3 4:3 12-6 20:0 30-1 26-7 33:6 28-9 24-1 

We see that after correction, assuming the extreme case, we are still left 
with large values of (P--Q)/(2N7z). Moreover, if the inefficiency of searching 
were the cause of the excess, (P—Q), then there should be a vast difference 
between the results of the “ice’’ plates and the other plates—the values 
of (P—Q)/(2 Nz) for the “ice ”’ plates would be much smaller than those 
for the rest of the plates. This difference was not found. This shows that 
the inefficiency in searching does not contribute much to the values of 
(P—Q). 

2. In this section we give some calculations of the possible causes of 
the * double stars ” and compare them with the experimental data. 
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(a) Regenerative process—If we assume that the double stars are pro- 
duced by penetrating radiation, emitted isotropically from the primary 
stars, producing secondary stars, then the number of double stars with 
separation <r is (Li and Perkins 1948). 


R=(2Nt) qo L [3d/4-+-d/2 log, (r/d)], 
for r>d, neglecting absorption. 


Here d =thickness of the emulsion=100p, 
L =number of nuclei (excluding hydrogen) per c.c. of emulsion 
=5-6 x 1022, 
q =number of penetrating particles from each star, 
o =cross-section of production of secondary stars by the penetrating 
particles. 


Since the values of 7 we are considering are >d, we can take r=r,. By 
comparing the calculated values of R/('N7) with the experimental values of 
(P—Q)/X Nz), we get 


go=(13°8 + 5-4) x 10-*4 cm?. 


A cross-section of the order of 10-24 cm. requires that g, the number of 
penetrating particles emitted per star, to be ~14. It seems almost impos- 
sible that there are so many particles emitted per star with energy large 
enough to produce another star. 

The values of g for known particles can be estimated as follows: In 
recent experiments with [ford G5 emulsion at 3560 m. altitude (Brown e¢ 
al. 1949), 1223 fast particles (with specificionization <1-5timesthe minimum 
ionization of a singly charged particle) were found to be emitted from 4923 
stars.. Allowing for an equal number of fast neutral particles, we get a 
value of 0-5 fast particle per star. Taking the number of grey tracks (ion- 
ization between 1-5 to 5 times the maximum ionization ofa singly charged 
particle) per star found for the star group with 7 and 8 heavily ionizing 
tracks as the average number of grey tracks per star for all stars, and multi- 
plying by 2-25 to allow for neutrons, we get a value of 3 per star. Assuming 
that half of these have energies large enough to produce another star, 
we get a total of 2 particles emitted per star, capable of producing a secon- 
dary star. Examples of secondary stars produced by fast particles emitted 
from primary stars have actually been observed with electron sensitive 
plates. Not all the fast particles end in this way. The contribution they 
offer, if assumed to be emitted isotropically, can thus account for less than 
2/13-8 or 15 per cent of the effect observed. 

(b) Single penetrating particles producing one or more stars in their 
passage through the emulsion. In this process, all stars are assumed to be 
produced by penetrating particles capable of producing more than one 
star before they lose all or a large part of their energies. A certain direc- 
tional distribution of the penetrating particles has to be assumed. An 
examination of a few particular cases will give us some idea of the cross- 
section required by such a process. 
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(i) If we assume that the penetrating particles are all coming in a 
direction perpendicular to the plane of the emulsion, it can easily be shown 
that the number of double stars is 

R=(Ni) d/2A, 
where \—mean free path between successive productions of stars in the 
emulsion. ; 

Taking 20/1000 as the value of (P—Q)/(2Nz), we get A=2°5 mm Or 
o=01 X10-** ema: ite tian 

(ii) Ifweassumethat the penetrating particles are all coming ina direction 
parallel to the plane of the emulsion, the number of double stars would be 


R=(2N2) 7/A, 
giving a value of A=9-3 cm. or o=2 x 104 cm’. 


Fig. la. Fig. 1. 
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(iii) If we assume that the penetrating particles are coming equally 
probably in all downward directions, then the expression for the number of 
‘double stars is simply the same as that for regenerative process except that 
q is replaced by unity, giving co=14 x 10-24 em?.. 

Case (iii) is presumably closest to the actual directional distribution, so 
that with an assumed value of c=10-*4 cm.2 this process could account 
about 7 per cent of the observed double stars. 

(c) Stars produced by sheafs of particles. If we assume that the stars 
are produced by sheafs of particles, which are generated in the material 
surrounding the plates and are emitted isotropically in all downward 
directions, we can calculate the number of double stars produced. 

Let K=number of sheafs generated per unit solid angle per unit 
volume per unit time, 
x —their average conical angle, 
S average number of particles per sheaf, 
T =total time of exposure, 
A =m. f. p. for star production in the emulsion, 
\’ =m. f. p. for absorption in the surrounding material, 
A =area of plates. 
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‘Then the number of stars produced in A is 


: 2 i aI Ua ced 7A cos 0 ; d sec 6 Am 
(EN)=| if iF iT (=) a (=) .p™ sin 6 dp db dd 


=2n KTAs 2 d/A. 


If two members of a sheaf simultaneously produce a star in the emulsion, 
we have a pair of stars with a projected separation <a«p sec 6. Let us take 
ap sec# < 27, as the condition defining a sheaf so that all stars produced 
in the emulsion by such a sheaf will have an average projected separation 
<r,. Then the number of double stars with projected separation Say his 


2% pn/2 72ry/xsec 0 Mi , 
R= i | | P KT. (AS) pifeelt AMS) (< =) _pisin 0 dp dé dd 
OF Os. 0 p 2 A 


‘2 (ie 
=%r KTA. as 5 as i hale —exp — (2 a =) | tan 6 dé 


ib sath (a) tan 6 dA 


Sp )=C ofp 


(s—1) 
2 


d 
ay 
(s-1) 
RAD 


R/(ZNi)= 


‘The integral can be evaluated graphically if we assign a value for «A’. 
Fig. 2 shows the curves of C.f(r,) against 7, for «A’=(1/2000) cm 
(1/200) cm., and (1/20) cm. respectively. They are normalized so that 
C.f(7,)=21/1000 at r,—(1/10) in. The factor C is then 0-0027, 0-0040, and 
0-0075 respectively. The small circles represent the experimental values of 
(P—Q)/(2Nz). It can be seen that the shape of the curve changes very 
slowly with the values of «\’._ It is not a sensitive method to fix the value 
of «X’ by finding the curve which fits the experimental data best. 

If we assume A=18 em. (corresponding to a cross-section of 10~*4 cm.? 
for the average nucleus in the emulsion), s can be calculated from the value 
of C. Thus we have 


xn’ (1/20) em. (1/200) em. ~—_ (1/2000) em 
fe 0075 -0040 0027 
8 28 15 ll 


We see that if we want to make s reasonably small, «A’ would have to be 
extremely small. It is very difficult, indeed, to explain the double stars by 
this process even if we assume that all stars are produced in this way. But 
suppose we accept that s=15 and A=18 cm. are reasonable, then we can 
explain «\’=(1/200)em. by assuming a small conical angle, say «=(1/1000) 
rad., giving \’=5 cm. which means that the particles are unstable, having a 
very short life-time. Thus we have tentatively the following set of values : 
A=18 em., ’=5em., «=(1/1000) rad., and s=15, which can roughly 
account for the observed number of double stars. 
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§4. CONCLUSION. 


It was seen in the last section that it is very difficult to account for the 
abserved number of double stars with our present knowledge. The 
isotropic regenerative process accounts for at most 15 per cent of the effect. 
Collimated penetrating showers, especially if in cascades, may also account 
for a fraction of the effect, though it could only be a small one owing to 
their large conical angles. It is, however, possible that these double stars. 
are due to something hitherto unknown. For instance, neutral particles. 
emitted in large number from primary stars producing secondary stars, or 
production of stars by sheafs of unstable neutral particles with small conical 
angles such as appear to produce the close pairs of electrons recently 
observed (Carlson et al. 1950). 
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SUMMARY. 


The electromagnetic signals used in communication are subject to the 
general laws of radiation. One obtains a complete representation of a 
signal by dividing the time-frequency plane into cells of unit area and 
associating with every cell a “ ladder ”’ of distinguishable steps in signal 
intensity. The steps are determined by Linstein’s law of energy 
fluctuation, involving both waves and photons. 

This representation, however, gives only one datum per cell, viz., 
the energy, while in the classical description one has two data: an 
amplitude and a phase. It is shown in the second part of the paper that 
both descriptions are practically equivalent in the long-wave region, or 
for strong signals, as they contain approximately the same number of 
independent, distinguishable data, but the classical description is always 
a little less complete than the quantum description. In the best possible 
experimental analysis by an electronic device the number of distinguishable 
steps in the measurement of amplitude and phase is only the fourth 
root of the number of photons. Thus it takes a hundred million photons 
per cell in order to define amplitude and phase to one per cent each. 


COMMUNICATION theory has up to now developed mainly on 
mathematical lines, taking for granted the physical significance of the 
quantities which figure in its formalism. But communication is the 
transmission of physical effects from one system to another, hence 
communication theory should be considered as a branch of physics. 
Thus, it is necessary to embody in its foundations such fundamental 
physical data as the quantum of action, and the discreteness of electric 
charges. This is not only of theoretical interest. With the progress of 
electrical communications towards higher and higher frequencies, we are 
approaching a region in which quantum effects become all-important. 
Nor must one forget that vision, one of the most important paths of 
communication, is based essentially on quantum phenomena. 

Some years ago, I have proposed a mathematical framework for the 
representation of signals (Gabor 1946). I have been rightly criticized 
for having left out noise, which is an essential feature of all communica- 
tions. This will be remedied here, and at the same time the description 
will be brought in line with modern physics. But as the mathematical 
frame will serve as a useful foundation, it will be necessary to give first 
a short review of it. 


* Communicated by the Author. 
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§1. CLASSICAL REPRESENTATION OF SIGNALS. 


The previous work started from the observation that the description 
of a signal in the conventional way, as a continuous function of time, is. 
redundant and non-physical. A continuous function contains in any 
interval, however small, an infinity of data, corresponding to an infinite 
range of frequencies. A similar objection can be raised against the 
Fourier representation, which involves infinite time. In the new 
description one considers the signal simultaneously as a function of 
frequency and time. It is convenient to use only positive frequencies in 
the description. This can be done by introducing a certain complex 
function, whose real part is the physical signal. (The theory of these 
complex or “ analytical ” signals has in the meantime received interesting 
additions by the work of J. A. Ville (1948)). If the time-frequency 
half-plane, fig. 1, is divided, by any network, into cells of unit area, 


frequency V —» 


Cime €— 


Information diagram. The time-frequency half-plane is divided up into cells. 
of unit area and an elementary signal is associated with each, with 
a coefficient cz. ; 


AiAv=1, one finds that the signal in any domain containing a sufficiently 
large number of cells is fully described by associating two real data, or 
one complex datum, with every cell. In other words, each cell has tw 
degrees of freedom. 

One can represent an arbitrary signal in an infinity of ways as a linear 
function of certain “ elementary ” signals, associated with the individual 
cells. There is, however, one description of particular interest, in which 
the elementary signals are harmonic functions, modulated with a 
“gaussian ”’ signal, 7.e., they have envelopes of probability shape. 
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(Fig. 2.) These share with other functions the property that their Fourier 
transforms have the same shape, but they are unique in that respect 
that the product of their “ effective ’ duration and of their effective 
frequency width is the smallest possible for any function. Thus, it can 
be said that these elementary functions overlap as little as possible. 
7 have also the advantage that the familiar concepts “ amplitude ”’ 
and “phase” can be used in connection with them in the same way as 
with infinite harmonic functions. The complex elementary function is 
Cos + ,Sin, if we denote by Cos the even, and by Sin the odd type of 
real elementary signal. Multiplying these with suitable complex 
coefficients ¢,,, as indicated in fig. 1, we can represent any arbitrary 
signal. Time-description and frequency-description (Fourier integral) 


Elementary signals of the “ cosine type ”’ (even), and of the “ sine type ”’ (odd). 


can be considered as extreme special cases of this representation. In the 
first case the even elementary signal degenerates to a delta-function, 
and the odd one to its derivative, in the second case both become ordinary 
harmonic functions. 

The “ matrix ” representation, illustrated in fig. 1, is proof against the 
objections raised against the pure “ time ” or “* frequency ” descriptions, 
but it does not go far enough. The infinity of data has been reduced to 
a finite number, but these data, 7.e., the coefficients c,,,, are still supposed 
to be exactly defined. But a single exact datum still contains an infinite 
amount of information, 7.e. an infinite number of “ yeses or noes” 
In reality, of course, these amplitudes, like every physical datum, have 
a certain amount of uncertainty or “noise.” This has been taken into 
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consideration in the mathematical theories of Shannon (1948) and 
Tuller (1949), where the noise amplitudes are assumed as known. But 
we cannot be satisfied with this in a physical theory. Even if all 
accidental imperfections of the instruments are eliminated, there remain 
certain basic uncertainties, which we are now going to investigate. 


§2. SraTisTicAL PROPERTIES OF THE INFORMATION CELL. 


In order to connect the mathematical scheme with physical reality, 
we must first observe that every physical signal has a certain energy 
associated with it. We can also associate a certain energy with every cell”, 
and, for simplicity, we will say that it is “‘ contained ” in the cell. 

We will limit the discussion to electric communications, though most 
of our results will apply equally well to sound, or, in short, to com- 
munication by any quantity which is considered as continuous in classical 
physics. We observe first that all electric signals are conveyed by 
radiation. Even if lines or cables are used in the transmission, by the 
Maxwell-Poynting theory, the energy can be located in empty space. 
Hence we can apply to our problem the well-known results of the theory 
of radiation. 

For simplicity, we consider our communication system as having the 
uniform temperature T. The uncertainty connected with the concept 
of temperature will produce certain fluctuations in the energy of the 
cells, which we can calculate by the rules of statistical thermodynamics 
once we know the law for the mean thermal energy €, of a cell, in function 
of the temperature. This we obtain at once, if we observe that every cell 
in the signal has two degrees of freedom, of which only one counts for the 
purpose of statistics, as the other is of the nature of a phase. Thus, b7 


Planck’s law 


2 hv 
at eo ee 


T exp (hv/kT)—V 
where v is the mean frequency of the cell. In other words, we identify 
every information cell with a “ Planck oscillator ’’. 

This requires perhaps a little more explanation, as physicists are less 
familiar with a discussion of radiation in terms of frequency and time 
than in terms of frequency and space. But the first case is immediately 
reduced to the second if we imagine the signal propagating with the 
velocity c, and plot the information diagram against the length ct instead 
of against the time t. The state of the field in such a linear system 


* The energies of the elementary signals of the type discussed will not add up 
exactly to the energy of the whole signal, because they are not quite orthogonal 
but the error will vanish in the mean over large number of cells. One can, 
however, to meet objections, consider instead an orthogonal set of elementary 
signals, such as the signals with “ limited spectrum ’’, introduced by Shannon 
(1949) and by Oswald (1949, 1950) which have uniform spectral density inside and 
ZerO outside a frequency strip. This makes no difference to the following 
discussion, as no reference will be made to any special type of elementary signal. 
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(in which we consider one state of polarization only) can be represented 
by two systems of progressive waves in opposite directions, only one of 
which represents the signal in which we are interested. We count the 
degrees of freedom in this wave system—which is what we have done in 
the last section—and give it the energy @, for each amplitude, disregarding 
the phases. This is the application to the linear case of von Laue’s well- 
known derivation of Planck’s law for the radiation density by super- 
position of plane waves (von Laue 1914). 

The thermal energy €, does not in itself represent ‘noise’, i.e., an 
uncertain disturbing factor. It becomes disturbing only by its 
fluctuations. In order to obtain the mean square fluctuations of the 
energy, we apply Einstein’s formula 

Sh—(eq— Pa kTe eT, SEO odin OF 
which gives peels ern. eee (3) 

Einstein’s interpretation of this equation is well known (cf. Born 1949). 
The second term has been identified by H. A. Lorentz with the fluctuations 
due to the interference of waves with random phases.* But the first 
term suggests that the energy is concentrated in light quanta or photons 
of energy hv, which fluctuate in any element as if they were independent 
particles. As both effects are present simultaneously, and are acting 
as if they were independent of one another, it is not possible to make 
a simple physical picture of the process. But fortunately this is not 
necessary. Applying Einstein’s argument to the case when a signal is 
present, so that the mean energy € in the cell exceeds the thermal mean 
energy €p, we obtain in Appendix I. 


S—=(e—ea ve Qeep—ep. ss se 4) 
Expressing the energy by the number of photons N in the cell, so that 
e=NAv, €>=N,hv, this becomes 

SN?=(N—N)?=N(1+2N,)—N2 . . . . . (5) 
with Nee Liexp iil) 1). . o. ac sue ee (6) 


* In technical theories of thermal noise, it is usually forgotten that it is not 
the noise power but its fluctuations which cause the disturbance. But if the 
quantum effect is small and the second term in equation (3) predominates, the 
r.m.s. value of the noise power fluctuation is equal to the noise power itself, 
hence this error is without consequences. 

From equation (3), neglecting the quantum term, one can easily derive 
Nyquists well-known rule (Nyquist 1928) that a resistance R can be considered 
as containing a “ noise generator ” with a mean square electromotive force 

H2=4kTRAp. 
The proof can be given in the same form as Nyquist has done, by substituting 
a cable with wave impedance R for the resistance. But one must add the 
condition that not only the noise power, but also its fluctuations follow the same 
rule in the resistance as in the cable. This is no arbitrary rule ; the necessity of 
a universal law of mean square fluctuations follows directly from the second 
principle, as Szilard (1925) has shown. 
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N,, the number of “ thermal ” photons per cell, is a large number for 
the frequencies used in electrical communications. One can use the 


approximate formula 
ut eed 
oe a A=0-7AT 
which gives Ny=210 for a wavelength of 1 cm. and a temperature of 
300°K. On the other hand for visible light, A=5-10-5 cm., the 


approximate formula 
N,=exp (—Av/kT)=exp (—1/0-7AT) 


gives N,=e-®4=10-*®. Hence for visible light, at ordinary temperatures, 
there is practically no thermal noise, and the fluctuation becomes pure 
“quantum noise ”, which follows the law 5N?=N. 

These results enable us to construct a complete physical representation 
of a signal. We see that the state of an information cell is completely 
determined by a stochastic number N, the number of photons. We can 
now construct a scale or “ladder” of distinguishable states, on which 
every step corresponds to a reasonably ascertainable difference. It is 
an evident suggestion to adopt the r.m.s. fluctuation of N as the unit 
step.* With this convention, the number of steps distinguishable below 
a maximum level N,, is, by equation (5), 


ae jae aN ={." aN 
— S Xr (5N2)12 J Ne [N(1+2N,)—Na}L? 


= —2_ (Ny, (1+2,) Nap? 10-4) 2} ye 
Teo N eek ie a 3 7 EN 
The last formula is valid for large signals. 

In the useful terminology introduced by D. M. MacKay (1950), S is 
the “* proper scale ’ of the photons. Fig. 3 illustrates the representation 
of a signal in three dimensions, with the photon scale at right angles to 
the time-frequency plane. If instead we plotted log, S, the ordinates 
would give directly the equivalent number of binary selections or “‘ bits ”’. 
This is the number of “yeses’’ or “‘noes”’ required to fix the position of 
the signal on the ladder. 

It is of interest to enquire about the minimum energy required for the 
transmission of the first “bit” of information. By our convention, this 
cannot be less than the r.m.s. value of the thermal energy fluctuations, 
though it can be more. Combining equations (1) and (3), one obtains 

a hv 
‘Cant = (deq) — 2 sinh (hv/2kT) . . . . . . (8) 


eg 

* By the theorem of Bienaymé and Tchebycheff, the probability of an error 
k times exceeding the r.m.s. error is smaller than 1/k?, whatever the law of the 
fluctuations may be. 


(7) 
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But, on the other hand, this energy cannot be less than one quantum, hv. 
Thus the energy required for the first ‘“‘ bit” is either <,,, a8 given by 
equation (8), or hv, whichever is the larger of the two. As shown in 
fig. 4, the two lines cross at hvy=0-96kT, which is only slightly less than kT. 


Fig. 3. 


Cine x 


Representation of an arbitrary signal. A ladder or “proper scale“ of 


distinguishable values is erected in every cell, on which the occupation 
is marked off. 


Fig. 4. 


0 %, 


The energy required for the transmission of the first “ bit ” of information. 


Thus up to a critical frequency v,,, an energy kT is sufficient for the first 
step, but no communication is possible with an energy of less than £T. 


AL2 
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The interesting feature of this result is its generality; it applies to the 
unknown processes in the nervous system as well as to electrical 
communications. 

§3. STATISTICS OF SIGNALS. 


Up to this point, we have directed our attention only to one cell. 
A short digression on signals and ensembles of signals may not be out of 
place before we return to the physical analysis of our results. 

A signal is a system of, say, n cells, preferably, but not necessarily, 
contiguous in the information plane. We consider a large number of 
such systems in a stationary transmission, in which they differ only by 
their position in time, not in frequency, and we speak of them as of an 
‘ensemble’. Evidently the analogy with statistical mechanics is not 
very perfect. In statistical mechanics, we can either follow a system in 
an ensemble over a long time, or we can look at all systems in the ensemble 
simultaneously, while here we can take only the first view. It is also 
somewhat questionable whether the often used expression “ ergodic ” 
is justified. In its original sense, due to Willard Gibbs, it means that 
each system (each signal) spends equal times in all states compatible 
with a given energy, which is not true for most stationary transmissions 
usually considered as “‘ ergodic’ in communication theory. Thus, we 
prefer to avoid this term. 

The most important mean value in such an ensemble is that of the 
entropy. In order to clarify the connexions between physics and com- 
munication theory, it may be useful to consider this problem in two stages. 
In the first, we consider all configurations of the system of 1 cells, 
compatible with the energetic conditions of the transmission as equally 
probable, and define the entropy as k log P, where P is the number of 
all these possible configurations. In the second stage, however, we give 
them different probabilities or “weights”. The first stage is in close 
connection with physics; actually it is the problem of calculating the 
entropy of an “ergodic ” system, in the original meaning of the word. 
According to quantum statistics, all simple, accessible states have equal 
probability, and the levels of Planck oscillators are, of course, simple. 
(Cf. Jordan 1933.) ; 

As an example, let us estimate the mean entropy of groups of n cells 
in a transmission in which the mean energy level is 8%, and the r.m.s. 
deviation from the mean is AS,,, for brevity. (The suffix n had to be added, 
as the standard deviation is dependent on the size of the group.) For 
simplicity, we assume (JS,,)?< 82, 7.e., a small degree of ‘‘modulation’’. 
This allows us to equate the number P of possible states to the number of 
points with positive, integer coordinates inside a spherical shell in 
n-dimensions, whose mean radius R is given by 

R2=n8?, 
while its thickness is 
24R=nA8,/R. 
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Well-known formule for the volume of a n-dimensional sphere give 
eclal4)i® 
~TGn+h 


(a7/4)a” eras Urls = AS 
= Perry Mi Se  endin VMAS, 


nR"-2AR_ 


We have used Stirling’s formula to approximate the I’-function for 
large arguments. For sufficiently large n, however complicated the 
signals, they must behave as if they were independent, and we must have 
asymptotically 


n.(AS,,)? > const. ; 


hence, apart from a constant, the entropy per cell goes to 


J lOe (se58) tee ae Anal) (10) 
a result also obtained by Shannon (1949). 

Apart from the factor k, the entropy (10) is also a measure of the 
information capacity per cell of the transmission, in accordance with its 
definition as the logarithm of the number of possible, equally probable, 
selections. Thus, the concepts of information capacity and of entropy 
are closely parallel, even identical. 

Consider now for a moment the central problem of statistical com- 
munication theory, which has been investigated by Shannon, the 
communication of symbols, as supplied by some source. If the source 
supplied them with equal probability, one could connect the possible 
and equally possible states of the signal with the symbols by an infinity 
of possible codes, and transmit information at the maximum rate, given 
by (10) multiplied by the number of cells per second. But if, as it is mostly 
the case, the symbols have different probabilities p,, it would be 
uneconomical to give them equal shares in the available signal entropy. 
Those which occur more often, and hence carry less information, should be 
given less than the more unexpected ones. In this way, though it is 
not, of course, possible to make the channel carry more information than 
given by (10), it is possible to make it carry more symbols, each symbol in 
the average being worth less. Shannon (1948) has shown that in the 
case of optimum coding each bit of channel capacity can carry 1/H 
symbols, where H is defined as 


H-=—2p los’p, bis/symbol, Sees. = .. (11) 


i.e., aS the mean value of a symbol in terms of binary selections. This 
“ selective entropy’ is defined by the source of symbols, and it is at an 
appreciable remove from the physical entropy of the signal previously 
discussed, as may be already seen from the fact that it appears as its 
divisor. The interesting properties of the expression (11), demonstrated 
by Shannon, follow rather from its mathematical form than from its 
intrinsic relation with Boltzmann's “H”’. 
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There exists, however, a general limiting relation between selective 
information and entropy in the thermodynamical entropy of the source 
considered as a physical system, which was revealed by Szilard (1929). 
This is that a signal corresponding in information value to an s-fold 
selection enables the receiving system to reduce the entropy of the 
transmitter by a maximum of k log s. It follows, as he proved in detail, 
that any mechanism acquiring this information must increase the entropy 
by a minimum of klogs, thus safeguarding the second principle of 
thermodynamics. 


§ 4. THE CONNECTION BETWEEN QUANTUM REPRESENTATION 
AND CLASSICAL DESCRIPTION. 


Comparing the representation of a signal, illustrated in fig. 3, with 
the classical, mathematical description, it appears that we have lost 
something. Previously we had two data per cell; we are left now with 
one only, a function of the quantum number N, which corresponds to the: 
energy. What has happened to the phase ? 

It is not difficult to give an answer to this question on general lines. 
N is not an exact datum but a stochastic number, which represents a. 
finite amount of information. But between one datum of finite accuracy 
and two, there is no transcendental abyss of the kind which exists between. 
a simple and a twofold infinity. One finite datum can very well contain 
two independent, finite data, provided that their aggregate information 
is not more than the original. It will be shown that this is indeed the case, 
and that the maximum amount of information on amplitude and phase: 
are both contained in the single photon scale. 

An electromagnetic signal can be physically analysed in various ways, 
but it will be instructive to consider first extreme cases only. One 
extreme is a counter, an instrument which records single photons. It 
follows immediately from Heisenberg’s uncertainty principle that the 
“time resolution ”’ of such an instrument cannot be better than a whole 
cycle, hence the phase remains entirely unobservable. The other extreme 
is any classical field-measuring instrument, capable of recording the 
electromagnetic field in the signal as a function of time. It may be 
called for brevity a proportional amplifier, as amplification of weak signals 
is one of its essential functions. There is no need to consider intermediate 
instruments, as it will be seen that in a certain range of very weak signals. 
every classical amplifier operates as a sort of ‘‘ proportional counter ’’. 

One limit for the operation of any such instrument is given by the well- 
known uncertainty relation of the quantum theory of radiation (ef. 
Heitler 1944, p. 68) 


ANAG Al 2 od Bee 
where AN is the uncertainty in the photon number N, and 4¢ the 


uncertainty in the phase ¢. But this gives merely an upper limit in our 
case, as we want to determine amplitude and phase simultaneously, 
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each as accurately as possible. It will be seen that, in fact, the limiting 
accuracy which is practically obtainable is much below what might be 
expected from the inequality (12). 

We will approach the problem by a detailed analysis of a particular 
type of proportional detector-amplifier. Subsequently, we will try to 
improve its performance to the extreme limit. It will then be found that 
all special features of the device vanish from the formule, thus the final 
result can be considered as of general validity. 

Assume that the electromagnetic signal is introduced into a rectangular 
wave guide, in the TE,, mode, 7.e., with an electric field in the x-direction, 
independent of the x-coordinate, which is at right angles to the direction 
of propagation z. An electron beam passes through two small holes 
along the z-axis, in the plane of maximum intensity, with velocity v 
and current intensity J. The alternating acceleration and retardation 


Fig. 5. 
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Analysis of an electromagnetic signal in a wave guide by an electron beam. 


of the electrons produces velocity modulation in the beam, which by 
well-known methods can be translated into current modulation, 7.e.,. 
producing an alternating component superimposed on the mean value J. 
One method utilizes the “ bunching ”’ of the electrons which takes place: 
at a certain distance from the guide, but one can also deflect the beam 
by a constant field and make it play between two collecting electrodes 
close together. In either case, one obtains an alternating current which 
in first approximation is proportional to the d.c. current and to the 
relative accelerations and retardations suffered by the electrons in the 
wave guide. 

Let us measure the energy exchange between the electrons and the 
field in quantum units, hv, where v is the mean frequency of the signal. 
(It will soon be seen that, in order to make the exchange intense, the 
waveband of the signal must be made so narrow that it is permissible 
to take the arithmetic mean.) Let N be the mean number of photons in 
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an information cell, which passes through a cross-section of a guide in 
a time dJt=1/v. During this time N=d At/e electrons pass through it, 
and exchange in the mean 7 quanta with the field, either by losing or 
by gaining energy. The positive number 7 is the essential parameter 
of the process. If % is a large number, which is possibie only if N is also 
large, the interchange will be essentially classic; if m is small, quantum 
phenomena will dominate.* 

The detailed calculations may be found in the Appendices II. and IILI., 
only the results will be discussed here. The first contains the classical 
theory of exchange, valid for large n and N, the second a wave-mechanical 
calculation valid for small exchange parameters. In the classical theory 
the result is 


E 32 /27e2\ sin? 0 v c \2142 Av = 
Be, rege Be) ee Ne oe les 
== (sr) Bb |S Gal i od 


Thus the exchange parameter is proportional to the square root of the 
photon number, as may be expected. Ofthe dimensions of the wave guide 
the width b appears explicitly in the factor of N, while the depth a is 
contained in the angle 

OS=nvalo OE ee 


which is one-half of the “transit angle”. c is the velocity of light. 
In addition, there appears a factor which is the reciprocal of 


the fundamental number which connects photons with electrons. 
Anticipating that the best measurement will require an intense 


interchange, we now try to increase the coefficient of N in equation (13) 
by all available means. First we make the factor sin? 6/9 a maximum. 
This is 0-723 and is obtained with 6—67°, i.e., a transit angle of 134°. 
This disposes of the depth a of the wave guide. The optimum width b} 
is determined by the condition that the group velocity 


Peeer aah 


must be as small as possible. But the smallest value is reached when U 
vanishes at the low-frequency limit, v—}Av, of the band. Substituting 
those values into equation (14), one obtains in the optimum case 


—, L:5(v\ fAv\t — 
ne = rav(s) (+) N. 5 ‘ E 3 5 @) 5) 


* Quantal energy exchange between electrons and the field in a wave guide at 
high quantum numbers has been previously discussed by Lloyd P. Smith (1946), 
but we cannot agree with most of his results. Monokinetic electrons and 
exchange of sharply defined quanta on the one hand, well-defined entrance phases 
and short transit times on the other, are mutually exclusive phenomena by the 
Uncertainty Principle; hence we believe that only certain averages over Smith’s 
detailed results have physical significance. 
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All special features of the device have vanished in this formula, apart, 
perhaps, from the unimportant factor, 1:5. But it is evident that the 
factor of N must always be much smaller than unity, while its best value 
as will be shown later, is just unity. There exists, however, a further 
possibility for improving the performance. Assume that we can make 
each electron perform repeated passages through the guide, each transit, 
in the opposite direction exactly half a cycle after the last. (This is 
possible in principle, as the optimum transit angle is about 134°.) If 
the frequency were known beforehand, the number of passages would 
be limited only by the consideration that, by repeated gains or losses, 
the electron would be bound to fall out of synchronism. But if the 
signal had a single frequency, known in advance, there would of course be 
no communication. However, even if the frequency is known beforehand 


Fig. 6. 
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Repeated passages. 


only within v--}Ayv, one.can make the number P of passages as great as 
v/ Av, without risking a phase error of more than +47, and it can be shown 
that these passages are almost of equal value, so that n is increased very 
nearly by a factor P. The number of passages required to make the 
coefficient of N in equation (15) unity is 


w\ U2 / p \1/4 
Pop =0°5 (2) Ga et enernes af (16) 


i.e., at least of the order ten. But this number must not exceed v/Av ; 
hence we obtain the condition that, in order to realize optimum conditions, 
the frequency band must be so restricted that 


(%) =202(5) fepiaiietets 0) %(17) 
Av Vv 


This is a somewhat surprising result. In the mathematical representation, 
it did not matter whether we divided up the frequency band into broad 
or narrow strips. But by the intervention of the number 137, 
(20-2=(137/1-5)%*), it turns out that only narrow frequency bands are 
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capable of accurate analysis by means of electrons. (If ions with charge Z 
were used, one would have to replace 137 by 137/Z?, and the condition 
would be less stringent.) 

Thus, at least in theory, the device could be perfected up to the 
optimum performance. The practical difficulties are of SOUnse evident. 
In practice one would rather replace the wave guide by a high-Q 
resonator, but this would somewhat complicate the theory. 

We thus find, assuming P passages, in the “ classical ’” case 


HOTS 09 


while the corresponding wave-mechanical formula, valid for very weak 
interchange, is (Appendix III.) 


é= ar(:) (2) ® eS 


Resi 


0 1 2 3 4 5 6 7 8 ) 10 


c 72 (?) Ay — 
187 \c) Ae 


Mean quantum exchange between electrons and photons. The initial tangent. 
is calculated from wave mechanics, the asymptotic parabola classically. 


Apart from the factor 2/3, the coefficient of N is the same in both cases, 
but this time the quantum exchange 7 is proportional to the photon 
number itself, not to its square root It can be said, therefore, that 
for small photon concentrations, the device acts as a counter, at large con- 
centrations, as a field-measuring instrument. The intermediate region is 
difficult to calculate, but, as shown in fig 7, the two branches can be 
connected by a plausible curve. 

One can also interpret the results in this way: If there are few photons 
present, there will be few collisions, and—as shown in Appendix III.— 
equal probabilities of gains and losses, at any instant. With increasing 
photon concentration, repeated collisions will increase in number, and the 
resulting loss or gain increases with the square root of the photon number 
only, but this resultant has now a prevailing direction, which changes. 
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its sign with the frequency of the signal. At this stage the “classical 
field’? has developed. _ 

Having ascertained that, with certain reservations, we can make the 
exchange as strong as we like, we ask the question: If we know the 
average photon number N, how must we adjust the electron beam 
current J, and the exchange parameter 7 in order to measure the field 
amplitude EK with maximum accuracy ? And having made these adjust- 
ments, how many steps shall we be able to distinguish in the scale of the 
sei amplitudes ¢ Evidently this question has a precise meaning only in 
the “ classical’ range of large n and N, and the following considerations 
relate only to this case. In order to simplify the problem, we neglect 
the thermal noise, 7.e., we put N,=0, so that the relative accuracy on 
the photon scale e would be 1/N1*, and the total number of steps in the 
photon ladder 2N12. The calculations are carried out in Appendix IV., 
here we give only the physical considerations 

The quantity to be measured is the electric amplitude in the information 
cell, which is proportional to the square root of the photon number. The 
measured quantity, on the other hand, is the alternating electron current, 
which, as mentioned above, is proportional to n.J or to n. M for not too 
strong signals, M being the mean number of exploring electrons per cell. 
For the optimum, we impose the condition that the relative mean square 
deviation of the quantity measured from the quantity to be measured 
shall be as small as possible, i.e., 


—— Smee ee, THe (20) 
(nM)? 

where the proportionality factor C is determined from the condition 
(nM—CN??)=0 


We have to choose n and M so as to satisfy the condition (20) for given N. 
There must be an optimum, for these reasons: A too weak interchange n 
would leave the cell unexplored. A too strong interchange, on the other 
hand, will interfere with the object of the measurement and spoil it. Though 
in the mean electrons are as often accelerated as retarded, fluctuations 
in the numbers M, and M, of electrons which pass through accelerating 
and retarding phases might produce extra photons, which could not be 
distinguished from those belonging to the signal, or annihilate some. 
The spurious photons are generated according to a law 


8,N=n(M,—M,)=n(SM,—8M,) . . . . . (21) 


as M,=M,=4M. We have written 5,N for this number, considering it 
as a fluctuation which must be eid to the natural fluctuation 6,N, 
whose law is 5,N?=N. The two fluctuations must be considered as 
independent. 

It is already evident from the above that the fluctuations in the number 
of beam electrons, i.e., the ‘shot effect’, play an important part in 
these phenomena. A too weak current has a high relative fluctuation. 
A too strong current, especially aided by a large exchange, will again spoil 


1176 D. Gabor on the 


the object. It may be noted that we have here a type of uncertainty whach 
springs directly from the fact that photons and electrons are discrete, without 
any reference to the physical values of h and of e.* 

The relative error according to equation (20) is calculated on the basis of 
equation (21), together with such evident assumptions as the independence 
of the fluctuations of n and M and of the ‘“‘natural” part of 3N. We 
assume also ‘“‘ normal’ shot effect, 5M’—M. The result is 


os 1/2\2 
GMS ONE ae (<5 + a): (22) 
(nM)? 4N %” mM AN 
This is a minimum for 
nM=nhyM=ON 2 72=N 0 x 2 


from which M=2N1?. This gives the simple rule that, for optimum analysis 
of the signal, one must take one electron for every step in the scale of 
the photons, and the interchange » must be itself equal to one 
distinguishable step at the level N. This again is a general rule, quite 
independent of the special model from which we started. 
Substituting these values, we obtain for the minimum of the mean 
square relative error in the measurement of amplitudes 
(oa) | pon aris 2 
See Se ee 
(nM)? min 4N n 4N Ni? 
As we are dealing with large photon numbers only, the first term at the 
right-hand side can be neglected with respect to the second. Thus we 
see, applying the same rule which we have used in constructing the 
photon ladder, that the smallest distinguishable relative step in the 
amplitude scale is 1/2 times larger than the square root of the corre- 
sponding quantity in the photon scale. In other words, the proper 
scale of the amplitudes will contain always less than the square root of the 
number of steps in the photon scale; apart from the factor \/2 for the 
reason that the optimum setting is, of course, possible for one level only. 
Having determined the proper scale of the amplitudes, a simple 
application of the Uncertainty Principle (Appendix IV.) shows that the 
proper scale of the phase must also contain less steps than the square 
root of the photon scale. Thus, summing up, we see that the classical 
description of the signal, by being too detailed, gives in fact a somewhat 
smaller total amount of information than the quantum description. 
The classical method of description, though theoretically inferior, may 
of course be still the best practically in the range of frequencies used for 
electrical communications, where efficient photon counters are’ not 


be Pa (24) 


‘ It may be noied, however, that the shot effect, that is to say the random 
arrival of the exploring electrons, is a practical rather than a fundamental 
limitation. As Prof. R. E, Peierls has pointed out to me, it must be ideally 
possible to measure amplitude and phase simultaneously up to the limit given 
by the relation 12. In fact if one could make the exploring electrons arrive at 
regular intervals, all terms but the first on the right-hand side of equation 22 
would vanish. But as it is practically impossible to suppress the shot effect to 
any significant degree, our results retain their practical importance. 


0 
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available. Conditions are different in the optical region, where detectors. 
of the counter type—such as the eye—are not far from perfection. In 
this region, analysis in terms of electromagnetic waves is as yet technically 
impossible, but it is interesting to note that, even if it were possible, it 
would not be very practical. Progress in the field of microwaves is now 
actually approaching a region where the two different methods of analysis 
may become competitive. We see from our results that it takes about 
a hundred million photons per information cell in order to define amplitude 
and phase of the signal to 1 per cent each. Remembering that at 1 cm. 
wavelength the number of thermal photons per cell is only about 200, 
it may be seen that the time may not be far off when the imperfections. 
of the classical method of description will manifest themselves even in 
electrical communications. 

Will it ever be possible to generate and to amplify light by an extension 
of microwave techniques? Generation may not be unthinkable, but 
we can show at once that a light-amplifier for the reception of weak-light 
signals is impossible. It follows from our results that an electromagnetic 
signal in which amplitude and phase are defined with a relative accuracy « 
must carry at least a power 

gear AV))) Maen eh eels A tek (20 
Applied to visible light with A=5000 A., a line width of 1 A. and «=0-01, 
we obtain 4:7 watts. Even if amplitude and phase were defined to 
10 per cent only, the power would still be 0-0005 watts, a very high intensity 
in terms of visible light. This example demonstrates that what we see 
are always photons, not waves. 

It may be hoped that these considerations have shown that the concepts 
of information theory may well prove their usefulness when applied to. 
problems of physics. 
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AGP PEON De 


Enxercy Friucruations In AN INFORMATION CELL IN THE PRESENCE 
oF A SIGNAL. 

As in the case of thermal equilibrium, the mean square energy 
fluctuation is the sum of a classical and of a quantum term. 

The classical term. Consider, for simplicity, a (complex) Fourrier- 
component KE, of the signal, and a component Ey of the thermal noise, 
corresponding to two different (circular) frequencies, say, a4 and wy. 
The instantaneous energy density, resulting from the interference, is 
proportional to 

E,Ef+E,E}+ {EgE¢ exp [j(w,—w»)t]+ ESE, exp [—j(w1—)t]}, 
where the asterisks stand for the conjugate complex values. The first 
two terms are the energy of the signal and the energy of the noise, the 
rest arises from interference. (Beats.) Let us write 

€=€5+ p+ egy. 
The mean value of the interference energy eg, is nil, but its mean square is 
ERE VADEI DES ASA ORR VE 
Using this and the relations 
esest— rest — 0 
which are evident, as there is no correlation between the signal and the 
noise, one obtains 


(<2)? [(eg Hep tegn) — (eg Hen P= G+ 2egey’ 
‘One knows the first two terms from Lorentz’s calculation, which gives, 
in the absence of a signal 


Diep 
ép— €p— eq, 


, hence (e—€)?=2ege, +e, 
or, as Eg =€— Eq, 
(c—€)?, = (2e—Ep)eq. . . . . . . . (1.1) 


This is the classical part of the fluctuation, arising from the interference 
of waves. The quantum term is always given by the “‘ law of rare events ” 
(e-e)2, Shve” oy % ieee (1.2) 


as if the energy were present in the form of particles (photons), with 
energy hv. Adding (1.1) and (1.2), one obtains equation (4). 


AUP Po RNeD- Lax. Ti 
‘CLASSICAL ENERGY EXCHANGE BETWEEN A WAVE AND AN ELECTRON 
Bram. 


In a rectangular wave guide with cross-section ab (fig. 5), the com- 
ponents of the electric field are, in the TE,, mode, 


H,=E, cos @ - COS 2zry (- 7) E,=E,—0, 
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where V is the phase velocity 


vf)" 


Thus the mean electric energy is, per unit length of wave guide, 


| 

gg Ho ab, 
and the mean total energy, electric and magnetic, is twice as much. 
This energy moves with the group velocity U=c?/V. The mean number 
N of photons in the e information cell is obtained by equating the energy 
flux expressed by N and expressed by E,. As the cell occupies a time 
At=1/Av, the flux per second, i.e., the power, is 


— ce \2 [1/2 
hvNAv = s— 3H jabU = st Te Bjabe| 1 = eee Ln 1) 
2bv 

Consider an electron which traverses the guide on the x-axis. Assume 

that its velocity v is large enough to neglect its change during the transit. 

Thus the number of quanta exchanged—gained or lost—will be, in the 
mean, 

\[*" wa 


n= 


b) 


[[" wal 
all ie|=— 


where the mean is to be taken with ache to the entrance time {). 
A simple calculation gives 


Pisce 


Fe mee eee ge =. (IL2) 
a 0 hp 


where @ is half the “ transit angle,” 7.e., 
0=n7valv. 


‘The condition 0<7 for the validity of (11.2) is always satisfied in an 
efficient arrangement, in fact, one must take 6 <j7 

Eliminating the field amplitude H, between (II.1) and (II.2), one 
obtains a direct relation between the photon number N and the exchange 
parameter 7 : 


eee oo, Bi c \2]~-12 Ap — 
i ae a | 1-(s5) | apy 


1 32 sin?6@ v c \2]712 Ap — 
Eo, Sg a ad | eas aan aos IN ee LE 
137 0 E (5) | ; ee) 


As explained in the text, the coefficient of N must be made as large as 
possible. The depth a appears only in the form sin? 6/0, whose maximum 
value 0-723 is reached at 067°, i.e., at a transit angle of 134°. The 
width 6b, on the other hand, appears in the factor 


yp. [1—(g. "= anu. 
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The smallest admissible value of U is obtained if the group velocity is. 
zero at the lower end of the frequency band, 2.e., if frequencies below 
=4A)y are cut off. In this case 


=}e|(v—}4r), 


hence the best value of the factor in question is the reciprocal of 


ce \2 41/2 C Av in 1 Av “as (Z)- 
st ()i], —aeeaan er alle ee 


assuming Av/v <1, i.e., the band to be narrow. Substituting these values, 


one obtains 
a 1-50 v /Av\12 
2 = (— mee etieate eae 
- Tal) a ue 


which is equation (15) of the text. 

If repeated passages are used, it is evidently permissible to consider 
them of equal value, so long as their number Pissmall. It may be asked 
when this assumption will lead to an appreciable error. An estimate can 
be made by introducing the concept of the “instantaneous frequency ’’* 
which will vary slowly between the limits v+3}4v. A somewhat long 
calculation, which may be omitted, gives the result that for large P’s the 
factor of n is not exactly P but approximately 


EES) o[-yee@- J 


where dy is the deviation of the instantaneous from the mean frequency. 
Assuming that the instantaneous frequency is uniformly distributed in the 
available band Av, one has 6v?=1/12(Adv)?. In the text, we have assumed. 
that the maximum number of admissible passages is P,,,,=v/4v. Even 


this gives an error of less than 14 per cent, which iaeen — approximate 
treatment. 


AS Ee be eNeD) es eee 


ENERGY EXCHANGE OF ELECTRONS AND WEAK FIELDS ACCORDING TO 
QuANTUM MEcHANICS. 


We apply to the problem the standard perturbation method of wave 
mechanics. As the beam width can be considered as small in comparison 
to the spatial periods of the field, it is sufficient to start from Schrédinger’s. 
one-dimensional equation with vector potential A,, 


oh OT Re OT eh Gh es 
ami Ot Sam Oa  Qame* on’ °° ° (HLL) 
We put 
Ap Ao COR Gg, Ay == A= 0, Wy=27, 
a 


* This concept is due to Helmholtz and has been first applied in communication 
theory by Balth. v.d. Pol. (1930). Cf. also J. R. Carson and T. ©. Fry (1937). 
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which corresponds to an electric field 
E,=Ey, sin wo, K,=E,—0, Hp=Ap/a¢; 


vq is here the frequency of the field, which previously we have called v, 
but in this calculation we will reserve v for the auteeucyy of the oscillators 
by which we will later replace the electron. 


Assume that the solutions for zero field are written in the form 
Pr=ty exp (—127e,t/h), 


where «, is the energy of the unperturbed electron in some state “n’”’. 
The general solution of equation (III.1) can be written in the form 


Pea eh) 
With this substitution, equation (III.1) becomes 
oa PS — oy Apa. se CSC) CTT.) 
mec ie 


Let us assume that the functions ¥,,(v) are orthonormal 


Multiplying equation (III.2) by Y*, integrating over the whole domain and 
making use of the orthonomality conditions, we obtain ordinary differential 
equations for the coefficients c,(¢) in the form 


d eA ‘ : 
Fe A! aeahe 2 mPnr 2SP [i(eo + Om) + exp [—i(%9—pm)tD} 
(IIT.3) 
h 0 
jee wah he = 
where Pur= 5 [OF zg hm Oa 
is the matrix element of the momentum p of the electron, and 


is the circular frequency which corresponds to the transition r— m. 

Assume now that at the instant ¢, the electron is in the state m=0. 
Thus at the start the coefficient c, is unity, all the others are zero. Let us 
also assume that during the time of interaction, 7.e., during the passage of 
the electron, the other states c,, increase slowly enough to neglect tran- 
sitions other than 0 +r. Thus we obtain, integrating equation (III.3) for 
the coefficient c, at the time ¢, 


neAy  fexp [i(wpte,)t]—exp [i(ot+oy)tol 
ee ae acta) Te 
exp [—17(wy—w,)t]—exp ee, 
“ (W9—,) ; 
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The absolute square of c, is the probability that after the time t—tf the: 
electron will be in the state “ r ”’, having absorbed an energy corresponding 
to the frequency w,»=o,. Substituting E, instead of Ay, we find finally 
for the probability distribution of the electron at time ¢ the expression. 


CN gee ea 


mhvo (wo +w,) (wo—o,,)” 
sin? $w,(t—tp) 
+2 cos Zwot oe . 5 (IIT.4) 


There are several interesting features in which this differs from the classical 

‘result. Note first, that the expression is even in w,, hence an electron 
has the same probability to gain or to lose a certain amount of energy. 
Note also that only the last term depends on the phase of the electric field 
at the instant t), the others depend only on the interaction time t—f). But 
even this last term has a frequency twice that of the electric field. It is 
rather doubtful whether this corresponds to an experimentally observable 
effect. Our calculation relates to weak fields, where n will be smaller than ~ 
unity; but in order to observe even n=1, the electron beam must be mono- 
chromatized to such an extent that its entrance phase becomes indefinite 
within a whole cycle. This is of no importance for our subsequent 
calculations, as we will consider only the mean value of (III.4), averaged 
over the entrance times ¢), in which the last term vanishes. 

There remains now the problem to apply equation (III.4), which is a 
well-known result of wave mechanics, to our problem of an electron 
traversing the oscillating field between two conducting plates. This means. 
calculating the matrix of the (mechanical) momentum p. We will evaluate 
this by a classical method, making use of the correspondence principle. 
This is justified in our case, as the de Broglie wavelength of the electrons. 
is always very small with respect to the distance a between the plates. 

In order to simplify the calculations, and to avoid singularities, we 
replace the electron by a plane density wave of total charge e (per unit 
cross-section) 


p(x, t)= Tear? [—(t—t)—2/v)2/7] D(x). 2. (IIL5). 
D(x) is an “ annihilation factor” which is unity inside the guide, i.e., for 
— $a <u“ <a, and zero outside. Apart from this, 7.e., if the electron is at 
more than about v7 from the walls, the distribution is gaussian. 7 is a 
small time, with which we will go ultimately to the limit zero. It will 
make the calculations easier if for the start we ignore the factor D(z), 
and make use of it only towards the end. 

We make use of the Fourier formula 


exp —[(t—a/v)/7 P= /(m)r { : exp [—(mrv)2—2riv(t—a/v)] dv. 
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This is to say that we decompose the electron into travelling waves of | 
density 


p, v=e/v exp —(mtv)? . exp [—2aiv(t—a/v)] dv. 


We have put the entrance time ty=0, as it plays no part in the mean values 
which we want to calculate. In the following formulz we will drop also 
the factor exp —(m7v)?, which cuts off the spectrum at very high fre- 
quencies, but which is practically unity in the frequency region to be 
investigated. 

The electric field produced by the electron follows from Poisson’s 
equation 


0K 
Let us decompose this too into harmonic components by putting 
K= [ K, dv. 
‘ .2€ : a 
One obtains Ke=3 = exp | —2niv G =) | 


for the progressive field waves.which accompany the free electron. But in 
a plane condenser with boundaries at «=-+4a, this oo induces surface 
charges 


M } : ‘ 1 
o,(+4a)=F Z, EAt2a) =F sa exP | —20iv G 5 =) | ; 


which cut off the field beyond x=-+-ta. 

We cannow calculate the electric moment of the charges in the condenser, 
i.e., of the electron charge and the surface charges, whose sum total is zero. 
The Fourier component of the moment in the frequency interval v . dv is. 


4a 
dv {| ape de talo,(}a)—o,(—4a))} =i Oe pom 5 Sin (mva/v) . exp (—2mivt) dv. 


a3 9? 


We now equate this to the “ standard oscillator ’’ of electromagnetic 
theory, 7.e., to an electron of charge e which oscillates with a frequency v 
around an equal and opposite charge, with an instantaneous amplitude X,. 
This gives 


X,dv=t iia —_ 573,2 8in (xva/v) . exp (—2rivt) dv. 


The mechanical momentum of this equivalent oscillator is mX, dv. Thus 
the momentum is, for the frequency interval v , dv, | 


P, dv= sin (mva/v) . exp (—2rivt) dv... .- (111.6) 
TV 


We cannot introduce this directly into equation (III.4), because that 
equation relates to discrete, steady states. In that case one has only to. 


4M2 
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take the absolute square of the matrix elements to form the quantities 
| Po, 2 which figure in equation (III.4). But here we have to do with 
continuous states, and we must also remember that the oscillators, 
according to equation (III.6), replace the electron during the finite transit 
time only. In order to obtain the equivalents of the “oscillator strengths ” 
| Mo, ?, we must now multiply the momentum amplitudes P, dv with the 
conjugate amplitudes P* du which belong to another interval » ,du; we 
must integrate over the whole domain of the ., and finally average over the 
transit time a/v. It is to say that we must calculate, instead of | 9%, [°, 


ga/v ioe) 
[PoP dv=25 dv [dt] PAPE du. 2 ae RL 


We had to add the factor two, because in wave mechanics v is considered as 
a positive quantity, while in the classical theory it runs from — oo to + oo. 
Substituting from (III.6), 


; mv\2 vy sin rva/v pa? © gin mpa/v 
pane (2) gnome ggj> eh 
— a/v —o 


a yp 
After integration with respect to t, the double integral becomes . 


exp [277(u—v)t] du. 


7 


1° sin apa i = 
=| sin 7pa/v eee sin (pu tae 
9 be pov 


With the substitutions v= (2u/v)—1 and 6=7va/v, this is transformed into 


7 


2 ,° sin $0(1+2) . sin $0(1—2) 2 ¢” cos @x—cos 0 i ha 
ai apes ete gee 
> (mv\2v sin? 6 
thus finally lo, az = 2(™*) sree ie iu gc te eee ee 


This means simply that during the transit time a/v the coherence region of 
any frequency v is v--v/2a, with a bandwidth v/a, which is the reciprocal of 
the transit time. A frequency v cooperates with any frequency , inside 
this band in producing the oscillator strength. 

Substituting (III.8) instead of | yo, |? into equation (III.4), where we 
drop the last term whose mean value is zero, and replacing t—t, by the 
transit time a/v, we obtain ; 


1 /27e?\ H2acv dv 
| ¢, F dy = il ) : 
h\ he 


Lee Py sin’ (vp+v)a/v sin? r(vy—v)a/v 
rae a ajo) (CSIs ee \ 
(1.9) 


This is the probability for an electron to be, after the passage, in a state 

v , dv, that is to say, to have absorbed the energy hv. Thus the energy 

spectrum of the originally monokinetic electrons after the passage through 
the field is of the form 

8(0) = sin? {5 (0-+-6,) ; sin? Ge 


a | O+0,)2 sie mo 0\ua (IIT.10) 
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where the frequencies (and energies) are expressed by the half transit- 
angles 


6=nva/v, Oy=7v,a/v. 


This function is shown in fig. 8 for a few values of 6). Experimental 
checking, though difficult, may not be impossible (cf. Macdonald and 
Kompfner 1949). 

The mean energy exchanged with the field, lost or gained, is by (III.10) 


1 EKacv 
137 dete * o> 


Tihvy =| hv |c, 2 dv = (IIT.11) 
0 


-180° 150 120 90 60 30 0 30 60 90 120, 150 180 


The function S(9), equation (III.10), which gives the energy distribution of 
electrons after the transit. Gains and losses have equal probabilities. 
The “resonance point’ at which the electron has gained or lost one 
quantum hy, of the oscillating field is at 6=6). 


where F(@,) is the integral 


© d@[ sin? (0+6,) | sin? (@—6)) 
= |r apart eae 
For 6=67°, which was the most favourable case in the classical theory, the 
value of the integral is found to be 1-22. Substituting this into equation 
(111.11) and expressing, as before the field E) by the number of photons, N, 
by means of equation (II.1), and finally giving the wave guide the same 
width 6 which was found the best in the classical case, one finds 


. 12 
= led (=) N, 


Vv 


where we have again written v for the field frequency, instead of vy). This 
is the energy exchange in weak fields, discussed in the text. 
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CONDITIONS FOR OptimuM EXPLORATION OF ELECTRIC FIELDS BY 
ELECTRONS. 
As explained in the text, we adopt as the criterion of optimum measure- 
ment the condition that 
(nM)? 


must be as small as possible. The proportionality factor C is given by the 
condition 


(nM—CN?2)=0, 


“ 


or, as we are dealing here only with the “ classical”’ case, in which the 
fluctuations are relatively small, C=nM .(N)-1?. As the fluctuations of n, 
the energy exchange parameter, and of M, the electron number, can be 
supposed as independent of one another, we can write nM=7n . M. 


Let us now first calculate 


n?M2—[(n-+-8n)(M+8M)2—(nM-+-nSM+M8n)?=72M2-+7128M?2-+ M25n? 
—WM?4+77°M+M 7. 
Here we have neglected the flunctuation term of the highest order, and we 


have assumed 
ondM=—0, on2=n, 6M?2=—M. 


The first is the assumption of the relative independence of the fluctuations 
of m and M. The second is the assumption that the absorption of a 
quantum is independent of those previously absorbed or lost, which is 
certainly admissible so long as the electron velocity is not changed 
appreciably. The third assumption means shot effect without space charge 
smoothing. This is a logical assumption in a theory in which each electron 
is considered to interact singly with the field. With these formule, we 
can now calculate 


(nM—CN1?)?—n?M?— 20 (7+ 8n)(M+8M)(N¥2-+ 3SN(N)-22) 


Oo 


+ C*(N224 38N(N)-12)?—= (7. M—CN?2)24 7 M-+M'n +402 = 


(IV.1) 
where the first term is zero. 


In order to calculate the last term, we consider the fluctuation SN of 
the photons as the sum of two independent components 8N, and 8N, 


SN=8N,+8N,, 


of which the first represents the ‘‘ natural ’’ fluctuations, subject to the law 
dN{=N, and the second is due to the exploring electrons, i.¢., it is the excess 
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of emitted over absorbed photons. Let M, be the number of beam elec- 
trons which pass through the field in a retarding phase, M, the number of 
accelerated electrons. As in the mean M,=M,—4M, and as in the mean 
each electron produces or absorbs 7 photons, we can write 


5,N=n(M,—M,)=n(8M,—8M,), 
and 5,N? =n (SME-+3M3)=7 (M,+M,)=n M, 


where we have assumed 6M,8M,=0, in accordance with our assumption of a 
normal, entirely random shot effect. 

We now substitute these results into equation (IV.1), and replacing C? 
by its value 72M’/N, we obtain for the mean square relative error in the 
measurement of the amplitudes 


(mM—CN?)? 1 1 4. n( i =) 


pars (IV.2) 


mM 4N?2 


me aN 
This is equation (24), discussed in the text. 
It is not so easy to give an exact analysis of the errors in the measure- 

ment of phase, because devices for the running measurement of phase, with 
uniform error, are by no means simple. But the order of magnitude can be 
seen immediately if we apply the uncertainty principle to the exploring 
electrons, in the form 

dedta 1. 
6t is the uncertainty in the time at which the electron passes some fixed 
plane, e.g., a fluorescent screen on which the waveform is recorded. The 
uncertainty in the energy « is, assuming that the electrons were mono- 
kinetic at the start, 

de=hvbn, 
n being the number of quanta which it may have absorbed from or lost to 
the field. But as 5n is of the order n1?, we obtain 


vdt=(n)-1?, 
In words, the oscillogram of the signal is traced with a spot whose mean 


width in the direction of time is a fraction 1/\/7 of a cycle. Thus one 
can determine the phase with the same order of relative error as the 


amplitude. 


ie Sses| 
CIV. Density Variations in Aluminium. - 


By J. L. SNonK, 


Philips Research Laboratories, 
Eindhoven *. 


[Received August 25, 1950.] 


SUMMARY. 
From careful determinations of the density of samples of aluminium 
in various conditions the conclusion is drawn that grain boundaries in 
aluminium may have a more spacious structure than is generally assumed. 


CoLtp WorK AND DENSITY. 


A priori considerations lead one to assume that the density of any metal, 
particularly however, that of metals with a close packed lattice structure 
will diminish as a result of cold work. 

Various authors have found such a decrease (Burgers), but very careful 
work by Smart, Smith and Phillips on copper has shown that the decrease 
observed becomes less and less as the amount of oxygen in the lattice 
diminishes so that in all probability in extremely pure copper there is 
only a very small effect of cold work on the density (<0-01 per cent). 

The above serves to show that one has to be careful in drawing con- 
clusions from results obtained on so-called pure metals. The only 
satisfactory procedure seems to be to vary the amount of a certain 
impurity and extrapolate the observed data towards concentration zero. 
In the majority of cases this has not been done and the results obtained 
remain open to serious doubt. 


Density AnD NUMBER OF VACANCIES. 


Our measurements on the density of aluminium in various conditions 
were started with the object of obtaining some direct data on the concen- 
tration of vacancies in this metal, the underlying idea being that in thermo- 
dynamical equilibrium this concentration would depend very much on 
the temperature and that it would be possible to ‘‘ freeze in” a certain 
number of vacancies by cooling rapidly from a high temperature, or, if 
the process of annihilation were very slow to eliminate them by a process. 
of cold work and annealing. 

Aluminium as an object for research was preferred to gold and silver 
or copper on account of its light weight, which makes the measurement 
of the density roughly four times more accurate than is possible for the 
heavier elements. Magnesium, which is even more favourable in that 
respect, was discarded because of its non-cubic structure and the difficulty 
in cold-working it satisfactorily. 


* Communicated by the Author. 
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Before describing the results, which were very unexpected, we shall 
give some details on the technique of the measurement and on the chemical 
composition of the samples. 


TECHNIQUE OF THE MEASUREMENT. 


Determination of the density on aluminium foil of half a millimeter 
thickness showed us—as was mentioned to us by Dr. W. Ch. van Geel of 
this laboratory before—that the surface absorbs variable amounts of 
water vapour. It is desirable, therefore, to make the surface play an 
unimportant part in the determination. 

The analytical balance used—a new one, capable of an accuracy of 
0-1 milligram—allowed of a maximum weight of 200 gram. As the 
standard shape of the object of which the density had to be determined 
we choose accordingly a cylinder with a diameter of 19 mm. and a length 
of 85 mm., having at one end a butt carrying a hole through which a. 
molybdenum wire of 50 micron thickness was passed, serving to suspend. 
the sample on one side of the balance. 

The determination of the “dry” weight of the sample to within 
0-1 milligram presented no undue difficulties, the weights used having been 
calibrated and the balance being placed in a room kept at constant. 
temperature. 

The determination of the weight in water presents more difficulties.. 
Doubly distilled water was placed in a massive block made of aluminium 
containing four channels, two vertical ones stopping above the bottom 
of the block and two horizontal ones connecting the vertical holes. The 
vertical holes had diameters of 50 mm. subs 20 mm. each. 

The wider vertical hole contained the sample, the smaller hole a Beckmann. 
thermometer allowing the temperature of the water to be determined 
to within a few thousandths of a degree centigrade. The temperature 
never deviated more than 2 degrees from the temperature to which all 
density determinations were reduced, viz. 20° C. The usual corrections. 
were of course applied. 

The sample as a rule was turned on the lathe only a few hours before a 
determination was made, but measurements taken some days later showed. 
no systematic deviations from this first value. Sharp corners were 
always avoided and the surface turned to a maximum of smoothness. 
The main difficulty lies in the capillary depression around the molybdenum 
wire. Olewm pedis bovis was finally chosen as a suitable means of reducing 
the uncertainties in the capillary behaviour to a minimum. 

The question whether a given determination is satisfactory from this 
point of view is answered by a determination of the apparent sensitivity 
of the balance under the conditions of the experiment. The adding or 
taking off of one milligram, should always cause the same displacement of 
the needle. If it does not, something is wrong with the wetting of the 


molybdenum wire. 
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Heat-treatments as a rule were carried out on pieces of aluminium 
bar diameter 20 mm. as taken from stock. After the treatment the sample 
was given a suitable shape for a determination of the density. Such 
determinations are called a-measurements. For comparison an additional 
heat-treatment was sometimes applied to a cylinder which had already 
served for a measurement of the density. On account of the oxide-layer 
formed it was then first turned off on the latheasecondtime. Measurements 
made on such a smaller specimen are called b-measurements. 


IMPURITIES IN THE ALUMINIUM USED. 

Aluminium was obtained in the form of long bars having a diameter 
of 20 mm. and a Vickers hardness of 29-2. Chemical analysis revealed 
the presence of the following impurity contents : 

Impurities in aluminium used 
Element Per cent by weight 


Si 0-006 
Cu 0-003 
Ga 0-002 
Fe 0-004 


‘The main impurity, however, is Mg, the amount of which was estimated 
from spectrochemical data to be at least 0-025 per cent by weight. 
No decrease in the Mg-content was observed after heat-treating the 
aluminium for many hours at 600° C. 


TaBLeE I, 
Number of Mean 
specimen Heat-treatment (in air) Vaph P29 grain-size 
(mm.) 
l(a) | As received 29-2 . | 2-69885 o- _- 
23) Wet We 35 
3 (a) | 50h. 400°C. 13-7 57 — 
FA} (ol einai, 59 °8 
5 (a) a »  +2h. 600°C. 14-4 81 — 
é(alu) Oa ‘ s6 54 
7 (a) Ps »  +20h. 600°C. | 13-3 78 96 -— 
8 (a) »? ” oe) 95 
9 (a) 2s »  +200h. 600°C. | 13-6 89) a5 a 
WA Sede. ce u 96 
1(b) | 50h. 400°C. 40 
OO 4502 ie 
a ae. 60 -- : 
12 (a) ” ” 50 uo 0 2 
9 (b) quenched from 600° C. in water 90 
10 (b) quenched from 600° C. in water 100 - : 
3(b) | 50h. 400° C.+200 h. 600° C. 78 
4 (b) 29 9? > be) 78 78 4 
5 (b) ? >? > 29 92 
6 (b) 99 29 29 99 95 94 4 
7 (b) » 0 99 10 ae 
ME Yh ‘ . 908 : 
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EXPERIMENTAL RESULTS. 


The outcome of our experiments was at first confusing and at any rate 
not in accordance with what we expected to find. It was found that 
heating at increasing temperatures caused the density to increase 
(Table I.) whereas quenching on the whole had no effect upon the measured 
density. The idea occurred to us that this might be an effect of the 
grain size on the mean density. 

A second series of experiments was then started in which the temperature 
of the heating was kept constant at 600° C., the times of heating being 
varied between wide limits. This led to essentially the same result 
(Table TI.). 


TaB_eE II, 
xnber me Heat-treatment Poo Mean 2091 
specimen S size (mm.) 
13 (a) 2 min. 600° C. (salt bath) 269841 <1 
14 (a) 5 76 1 
15 (a) 10 68 2 
16 (a) 25 76 2 
17 (a) 50 80 3 
18 (a) 100 81 3 
TaBLE III *, 
(Fine grained as cast.) 
. H Mean grain- 
Specimen eat-treatment Pao oe 
A As received 2209919) yl oy el 
A 2h. 600° C. BS TOTCd wee OLS 6 
B As received 2 G9938 a: <1 
B Zh. 600° C. D002Guaay dvi 4 


As a third approach we had a cast made from a fairly cold melt, chilling 
it rapidly so that a fine grained structure was obtained. This material 
also became very coarse after a prolonged heating at 600° C. and by this 
treatment again became denser (Table III.). The aluminium used for 
this experiment came from another stock and contained less magnesium 
(0-005 per cent). 


DISCUSSION OF RESULTS. 


The above-mentioned results are hard to explain. One may first of 
all think of a spurious effect caused by impurities present in the material. 
The pursuit of this thought leads, however, to no tangible results, in view of 


*The circumstance that the density of the cast samples on the whole is 
somewhat larger than for the extended bars may be due to the fact that it 


contained less magnesium. 
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the three very different series of heat-treatments applied to the aluminium,. 
which all lead to essentially the same result. The cast sample besides 
contained considerably less magnesium than the wrought samples. 

A second possibility is the appearance and gradual disappearance again 
of mysterious holes or gaps. We did not find a trace of such gaps on 
micro-etched samples. 

Somewhat reluctantly we arrive therefore at the conclusion that grain 
boundaries in aluminium have a much more spacious structure than is 
generally supposed. They might, for instance, act as “ retainers ” for a 
number of vacancies which are formed at high temperature. No specific 
picture can be offered at the present time. | 

It is hoped that our experiments will be repeated by other physicists. 
with a view to obtaining an independent check of the result presented. 
here. 


ADDENDUM. 


These results have been discussed with Professor N. F. Mott and 
Dr. F. C. Frank who point out that if the density changes are attributed 
to an open structure at the grain boundary they imply that each boundary 
is equivalent to a void layer 4(5 x 10-? x 2 x 10-4) em=330 A thick, which 
is hard to believe. It is almost as difficult to attribute the changes to 
precipitation of solutes, for if the whole of the analytically ascertained 
impurity were going into and coming out of solution, it would have to be 
changing volume by about 50 per cent. They have suggested that the 
aluminium specimens possibly contained gas (most probably hydrogen). 
which goes into solution at 600° C., stays in under ordinary cooling and 
particularly with water quench, but comes out again within 50 hours at 
400° C. in the form of gas bubbles which lower the density. 

This explanation cannot be entirely excluded, especially with regard 
to those samples (the majority) which came from a stock not prepared. 
in the laboratory. The third series of measurements were carried out on 
a cast which had been made in the laboratory under the exclusion of free 
hydrogen, and in the presence of only very little water vapour. It seems 
unlikely that this cast would contain hydrogen in appreciable quantity. 
No holes were seen in microscopic examination of specimens after electro- 
polishing, though one would expect to see them if gas bubbles were 
present to the extent required (one bubble of 30u, 100 of 3u or 10° of 
0-3 diamenter per grain). 

Whatever the cause, the inherent difficulties connected with experiments 
of this kind on aluminium, and the desirability of an independent check, 
are emphasized. 
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CV. CORRESPONDENCE. 


Radiative Corrections in the Non-Relativistic Limit. 


By W. Turrrine, D.Phil., 
Dublin Institute for Advanced Studies *. 


[Received September 10, 1950.] 


Iw calculating radiative corrections to the Klein—Nishina formula it has 
been assumed by Corinaldesi and Jost (1948) and by Heitler and Ma (1949) 
that the corrections vanish in the non-relativistic limit. All terms 
proportional to the second-order matrix element were interpreted as 
renormalization, which means defining the electron-charge by the 
Thomson formula. The question is now, whether Dyson’s (1949) 
renormalization procedure, the only one of which the consistenéy has 
been proved, leads to the same result, or whether there are finite 
‘corrections proportional to the Thomson cross-section analogous to the 
‘case of the magnetic moment. It has been shown by Schafroth (1949), 
by direct calculation, that after renormalization the fourth-order matrix 
element vanishes for zero momentum of the incident photon. We will 
give a simple proof that this holds in all approximations. 

One obtains the graph of the nth order to the Compton effect by 
inserting successively the two external photon lines in any of the 
electron lines of the electron self-energy graphs of order n—2. An 
electron self-energy graph consists of the main electron line and some 
closed electron lines and internal photon lines. If one inserts the two 
external photon lines into the main electron line (energy-momentum py) 
for the limit of zero momentum of the photon, using 


d L 
ap (p—m)-!=—(p—m)“1y,(p—m)* ; p=?,7' so gt (1) 
Le 
and summing over all insertions and all self-energy graphs, one obtains 
7) 
— Hoe : — Bete eth aA gly a2 
(p—m)| ae a8 (») | (em) (2) 


for the matrix element. After renormalization S’(p), as defined by Dyson, 
has the form 
S'(p)=(p—m)-“[1+ B(p?—m?)]+m*C(p2—m*) fF, « - (3) 


* Communicated by the Author. ; 
+ Dyson’s renormalization casts contributions from proper self-energy parts 


into the form (3) whilst on multiplying the proper self-energy parts the 
contribution from the general self-energy part will be 

— (p—m)4{1—B—Cm™(p—m)]™ 
but one readily sees that this is again of the form (3). 
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where B(0)=0. Inserting into (2) we get 
(2)=[e'(p—m)-te?-++ e°(p—m)-1](1-+ B(0)) —2[e(pe*) + e*(pe?) |B '(0) 
+4(e1p)(ep)[(p—m)B”’(0)+ (p—m)*m*0"(0)]. 
Everything but the uncorrected matrix element vanishes since one can 
choose for the potentials of the external photons a gauge system in which 


(pe')=(pe*)=0. 
Contributions from the closed loops will have the form 


| | dq trIT(q-+k,—m)-*;, 


where h, are some energy-momentum fourvectors, over some of which 
one might have to integrate. Inserting one external photon line 
successively in all electron lines one gets for zero momentum 


—e, { dq = trIT(q+-k;—m)~*e;, 
qu ry 
which vanishes on partial integration, provided that the integral 
converges. Otherwise one has to apply regularization in order to avoid 
contributions from the boundary. 

It is evident that the same argument applies to a case where the electron 
interacts any number of times with the electromagnetic field, and that 
therefore in all those cases the radiative corrections vanish in the 
non-relativistic limit. 
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Multiple Processes in the Production of Pairs of Electrons by y-Radiation. 


By J. E. Hooprr and D. T. Kine, 
H. H. Wills Physical Laboratory*. 


[Received September 28, 1950.] 
[Plate XXXV.] 


Ir is shown by Heitler (1944) that the cross-section for the simultaneous 
emission of » quanta by a multiple Bremsstrahlung process is, at high 
energies, of the order Z?r?/137". On this basis one would expect the 
cross-section for the analogous case of the production of pairs of electrons 
by a y-ray in the electrostatic field of a nucleus to be of the same order 
of magnitude. 


* Communicated by Professor C. F. Powell. 
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We have recently made an investigation of the pairs of electrons 
created in photographic emulsions exposed at an altitude of 68,000 ft., and 
resulting from the passage of y-radiation arising in the decay of 7-particles. 
Among 1380 pairs we have found two examples which appear to corre- 
spond to the simultaneous production of two pairs of electrons. A mosaic: 
of photo-micrographs of one of these events is shown in Pl. XX XV. In both 
cases the associated tracks are of sufficient length to allow us to determine 
the grain-density in the tracks and the momenta of the particles by the 
scattering method. All the particles have a specific ionization equal to 
the minimum value for a particle of charge |e}, and the total energies of 
the particles in the two events are 220 and 2000 MeV. respectively. 

The energies of the individual particles are shown in the Table. 


Event (A) Event (B) 
oe 40+5,  95-+10, 8+3, 960-4150, 
(MeV.) 60-8, 22+ 4 20+4, 980+150 
Ey (MeV.) 220+25 2000-4300 


These measurements show that nearly all the particles are less massive 
than protons. Further, in one case two and in the other, three of the 
particles are less massive than pz.-mesons. These facts, and the observed 
even number of secondary particles, allow us to exclude the possibility 
that either of the events is due to the collision of a fast neutron with a 
proton in the emulsion in which a shower of z-mesons is created. We 
therefore attribute them to the production of pairs of electrons and 
positrons. 

The possibility remains that the events are due to a process first observed. 
by Bradt and Peters (1950) and by Occhialini (1949), in which one of the 
electrons of a pair, created in the emulsion by a y-ray, subsequently 
generates a second pair in the field of a nucleus. We have observed 52 
examples of pair creation by a fast particle—the so-called “ tridents ”— 
and it would be possible to explain the present observations in terms of a 
similar occurrence in which the production of the second pair has occurred, 
by chance, very close to the point of origin of the “ primary ” pair. Such 
an event would be difficult to distinguish from one involving the simul- 
taneous production of four electrons in the field of a nucleus, especially 
in the case of pairs of small angular divergence such as are commonly 
produced by y-rays of great quantum energy. We believe, however, that 
the following argument allows us to exclude this possibility :— 

The number of observed events in which a “trident” is found to 
be associated with a pair of electrons enables us to estimate the average 
‘mean free path” in the emulsion for trident production. The mean 
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value is about 100em. The probability of the trident occurring within 
10 of the point of origin of the pair is therefore 1 in 104; and the 
probability of observing two such events among 1300 pairs is about 
1in 105, These values are so small that we can disregard this explanation 
of the origin of the two events under discussion. Further, the probability 
that, by chance, two unrelated pairs should be superposed is less than 
1 in 102 per cm.2 of the emulsion in the conditions of the experiment. 
We therefore regard our observations as providing strong evidence for 
the multiple production of pairs of electrons by a y-ray quantum in the 
field of a nucleus, the probability of the process being of the order of 
magnitude suggested by quantum-electrodynamics. 

The observations on the “ tridents ”’ associated with pairs of electrons 
will enable us to determine the variation with energy of the cross-section 
for the production of “ pairs” by fast electrons. Preliminary results 
suggest that there is a rapid increase in the value of the cross-section for 
the process as the energy of the incident electron increases. We may 
direct attention to the fact that some caution should be observed in 
observations of this type because of the possibility of confusion due to 
‘““Bremstrahlung’’. A fast electron loses energy in inelastic collisions 
which result in the creation of quanta. Especially in the case of electrons 
of high energy, these quanta may be emitted in directions nearly 
coincident with the line of motion of the recoiling electron. If, further, 
such a quantum produces a “ pair ’’, the resulting electrons may appear 
to arise at a point nearly coincident with the track of the original electron. 
Such events can be distinguished only with difficulty from true trident 
production. The possibility of this source of confusion was first pointed 
out to us by Professor Peierls at the Edinburgh Conference in 1949. 
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CVI. Notices of New Books and Periodicals received. 


Practical Applications of Spectrum Analysis. By H. Dryeue. [Pp. 245.] 
14 plates. (Chapman and Hall, 1950.) Price 40s. net. 


Ix a general conclusion to this useful book the author states that in his 
experience more than half of the spectrographic analyses which have been 
published during the last thirty years are untrustworthy. Used in conjunction 
with a standard set of wave-length tables, this book should do much to remedy 
this state of affairs. 

No attempt is made to describe the theoretical laws of spectroscopy, but 
the relevant aspects of the mechanism of radiation, of sources of luminosity 
and of the optical principles of the spectrograph are dealt with in so con- 
structive a fashion as to lead to a confident and critical choice and use of the 
most appropriate apparatus for a very wide range of spectrographic analyses. 


- Vol, 41. Pl XXXV. 
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The chapters on the measurement and interpretation of the spectrograms 
obtained in these analyses, together with the comprehensive tables of Raies 
Ultimes and the wealth of critical comment on special points of interpretation, 
taken element by element, embody great experience and form a valuable guide 
to the less experienced and an indispensable reference work to the more 
experienced workers in this field. 

The writing throughout is clear in style, and the general level of production 
excellent, except for the reproduction of the plates, the quality of which is 
somewhat disappointing for a book at this price. D. McG. 


Progress in Coal Science. Edited by D. H. Bancuam. [Pp. 456.] (Butter- 
worths Scientific Publications, London, 1950.) 


Tats book comprises twenty-six review articles, many of which have appeared 
in the Monthly Bulletin of the British Coal Utilisation Research Association 
during the years 1943-1947. They are divided into the following sections :— 
Modern Experimental Techniques, Fine Particles, Constituents of Coal, Organic 
Chemistry of Coal Products, and Chemical Aspects of Combustion and Gasifica- 
tion. 

The first section contains articles on crystallography, absorption spectro- 
graphy, mass spectrometry, and chromatography, all as applied to chemical 
analysis. These should be of use to chemists generally, inside and outside the 
Coal industry, as also should the articles on the production and measurement of 
fine particles. The third section contains among others an account by 
Prof. V. M. Goldschmidt of his work on the rare elements in coal ash, and a 
review by Dr. C. A. Seyler of his work on coal petrology. In the fourth section 
is a translation of a German article by O. Gruber on “ The Progress in Coal Tar 
Research in the period 1910 to 1940,” and a review of the Fischer-Tropsch 
reaction. The last section contains an article by Prof. D. T. A. Townsend, 
which gives a general account of producer gas research, and more specialised 
articles on fuel beds, combustion kinetics, and other topics. 

The book should be most valuable to all scientists in the coal and allied 
industries and, as indicated, some of the articles should find a wide pitta 


Astronomy of Stellar Energy and Decay. By Martin Jounson. [Pp. 216, 
Plates 4.] (Faber and Faber, 1950.) Price 18s. Od. 


Tuts book, which aims at bridging the gap between popular astronomy and 
the research literature of stellar physics, can be strongly recommended. It is 
deliberately set in two parts. The first, though framed for the general intelligent 
reader, can also be read with interest by students of wider scientific training. 
Part II. is labelled ‘‘ student’s introduction to analysis of the facts and argument 
of the theories.” As its title suggests it is more qualitative, though never 
heavily mathematical, in its treatment. It supplements Part I. and serves as a 
very useful introduction to more advanced. treatises. The two parts together 
can be specially recommended for private reading at the university degree 
stage. ee Mo: 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 
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